Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 1,677,492 Δ1 bp coding (412/1533 nt) pntA ← pyridine nucleotide transhydrogenase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,677,4910C.96.2% 362.1 / 7.5 79coding (413/1533 nt)pntApyridine nucleotide transhydrogenase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (37/39);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCA  >  NC_000913_3_pae_pgi/1677345‑1677635
                                                                                                                                                  |                                                                                                                                                
gaCTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAAcggc                                                                                                                                               <  1:923211/150‑4 (MQ=255)
         ctgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGc                                                                                                                                     <  2:336750/150‑1 (MQ=255)
            cgcCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCGAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCAt                                                                                                                                  >  3:256499/1‑150 (MQ=255)
             gcCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATc                                                                                                                                 >  1:6199/1‑150 (MQ=255)
              gcAATGGCTGCTCGACTTGCAACACCCGCACCAATCGCCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCg                                                                                                                                <  4:406068/149‑1 (MQ=255)
              ccAATGGCGGCCAGACCTGCAACACCCGCGCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTGAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCg                                                                                                                                <  2:746301/150‑1 (MQ=255)
              ccAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCg                                                                                                                                >  2:32711/1‑150 (MQ=255)
               cAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGa                                                                                                                               >  1:1105664/1‑150 (MQ=255)
               cAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGa                                                                                                                               >  4:9056/1‑150 (MQ=255)
                aaTGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAg                                                                                                                              >  3:495061/1‑150 (MQ=255)
                  tGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCt                                                                                                                            <  2:1105664/150‑1 (MQ=255)
                    gcggcCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTa                                                                                                                          >  3:434328/1‑150 (MQ=255)
                      ggcCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        >  1:1097752/1‑150 (MQ=255)
                      ggcCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        >  1:506198/1‑150 (MQ=255)
                          aGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGt                                                                                                                    >  1:248225/1‑150 (MQ=255)
                           ggccTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                                   <  2:346241/148‑1 (MQ=255)
                           gACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                                   <  1:399712/150‑1 (MQ=255)
                           gACCTGCAACACCCGCACCAATCACCAGCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                                   <  1:1143923/150‑1 (MQ=255)
                              cTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAg                                                                                                                <  2:546202/150‑1 (MQ=255)
                                         gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCGAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                     <  4:256499/150‑1 (MQ=255)
                                         gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                     >  2:987201/1‑150 (MQ=255)
                                         gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                     <  2:271244/150‑1 (MQ=255)
                                          cACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGcgcg                                                                                                    >  3:176587/1‑150 (MQ=255)
                                           aCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                                   <  2:1097752/150‑1 (MQ=255)
                                           aCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                                   <  2:506198/150‑1 (MQ=255)
                                            ccAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                                   <  1:718210/150‑1 (MQ=255)
                                               aTCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGAt                                                                                                                                           >  3:637553/1‑106 (MQ=255)
                                               aTCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGAt                                                                                                                                           <  4:637553/106‑1 (MQ=255)
                                                tCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGAt                                                                                              >  3:652884/1‑150 (MQ=255)
                                                  aCCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAc                                                                                            >  3:454059/1‑150 (MQ=255)
                                                   ccATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg                                                                                           >  4:599343/1‑150 (MQ=255)
                                                    cATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcgc                                                                                          <  2:248225/150‑1 (MQ=255)
                                                     aTCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        >  4:340515/1‑119 (MQ=255)
                                                     aTCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        <  3:340515/119‑1 (MQ=255)
                                                     aTCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcgcg                                                                                         >  4:487423/1‑150 (MQ=255)
                                                     aTCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcgcg                                                                                         <  3:9056/150‑1 (MQ=255)
                                                      tCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCgg                                                                                        >  4:654308/1‑150 (MQ=255)
                                                          ttttGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGtt                                                                                                                                        <  4:34927/98‑1 (MQ=255)
                                                          ttttGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGtt                                                                                                                                        >  3:34927/1‑98 (MQ=255)
                                                           tttGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACag                                                                                   <  1:192035/150‑1 (MQ=255)
                                                               ccGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTc                                                                               <  2:1126105/150‑1 (MQ=255)
                                                                       aCTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCa                                                                       <  4:728/150‑1 (MQ=255)
                                                                          ttCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCa                                                                    <  4:319373/150‑1 (MQ=255)
                                                                             ccGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACgg                                                                 >  4:386760/1‑150 (MQ=255)
                                                                                 ccGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCac                                                             <  2:103645/150‑1 (MQ=255)
                                                                                   gCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgt                                                           >  1:1113708/1‑150 (MQ=255)
                                                                                      gTAATTTGCCCGGTAAAGAAGCGCCCAAATGCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgttac                                                        <  2:288752/150‑1 (MQ=255)
                                                                                       tAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                                   >  4:269442/1‑90 (MQ=255)
                                                                                       tAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                                   <  3:269442/90‑1 (MQ=255)
                                                                                        aaTTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgttacgt                                                      >  4:153588/1‑150 (MQ=255)
                                                                                            tGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCg                                                  <  2:79132/150‑1 (MQ=255)
                                                                                               ccGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCaa                                               <  1:76614/150‑1 (MQ=255)
                                                                                                cGGTAACGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTcc                                                   >  4:558030/1‑146 (MQ=255)
                                                                                                cGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTcc                                                   <  3:558030/146‑1 (MQ=255)
                                                                                                       gttacgcCCAAATCCAAGTGCAGCTTCAACAATCGCGCGATAACGGCGATGTCGCCCGGCGAGCTTCGTGCGGCCAGCGATTGTCCGCGTGAGATACGCGGCACAGAGTACATCGCCTTCACGGTCACGTTACGTTCCCCCAGATGTTGc                                       >  2:467953/5‑150 (MQ=255)
                                                                                                       gaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        >  2:713643/1‑69 (MQ=255)
                                                                                                       gaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                        <  1:713643/69‑1 (MQ=255)
                                                                                                       gaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGc                                       <  4:309001/150‑1 (MQ=255)
                                                                                                          gcgcCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCAtt                                    <  4:82208/150‑1 (MQ=255)
                                                                                                           accccAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTa                                   <  4:410079/148‑1 (MQ=255)
                                                                                                             cccAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAAt                                 <  1:482928/150‑1 (MQ=255)
                                                                                                             cccAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAAt                                 >  4:731412/1‑150 (MQ=255)
                                                                                                              ccAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAAtt                                >  3:536783/1‑150 (MQ=255)
                                                                                                               cAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTc                               >  1:487246/1‑150 (MQ=255)
                                                                                                               cAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTc                               >  2:800693/1‑150 (MQ=255)
                                                                                                                   ttCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGa                           >  3:512260/1‑150 (MQ=255)
                                                                                                                     cATGTGCCGCTTCAACAATGGGGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGAtt                         <  4:25332/150‑1 (MQ=255)
                                                                                                                      aTGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTc                        >  1:831417/1‑150 (MQ=255)
                                                                                                                        gtgCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTg                      <  2:1100342/150‑1 (MQ=255)
                                                                                                                         tgCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTgc                     >  1:539372/1‑150 (MQ=255)
                                                                                                                             gCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAg                 >  4:2531/1‑150 (MQ=255)
                                                                                                                                 caacaaTGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCa             <  2:487246/150‑1 (MQ=255)
                                                                                                                                  aacaaTGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCAGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAg            >  3:623454/1‑150 (MQ=255)
                                                                                                                                     aaTGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgc                                                                                                       <  2:513076/56‑1 (MQ=255)
                                                                                                                                     aaTGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgc                                                                                                       >  1:513076/1‑56 (MQ=255)
                                                                                                                                        ggCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa      <  3:113123/150‑1 (MQ=255)
                                                                                                                                        ggCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa      <  4:454059/150‑1 (MQ=255)
                                                                                                                                        ggCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa      <  3:654308/150‑1 (MQ=255)
                                                                                                                                         gcgcGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAAc     >  4:600907/1‑150 (MQ=255)
                                                                                                                                            cgATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCa  <  4:176587/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCA  >  NC_000913_3_pae_pgi/1677345‑1677635

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: