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.94.2% 384.9 / 17.6 87coding (413/1533 nt)pntApyridine nucleotide transhydrogenase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base . (39/43);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAA  >  NC_000913_3_pae_pgi/1677342‑1677631
                                                                                                                                                     |                                                                                                                                            
cGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAAc                                                                                                                                              >  2:389927/1‑150 (MQ=255)
   gaCTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAAcggc                                                                                                                                           <  4:563866/150‑4 (MQ=255)
   gaCTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCgg                                                                                                                                           <  3:252708/150‑1 (MQ=255)
   gaCTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCgg                                                                                                                                           <  1:446847/150‑1 (MQ=255)
      tGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGAt                                                                                                                                       >  4:300764/1‑150 (MQ=255)
           gctgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCg                                                                                                                                  <  1:389927/150‑1 (MQ=255)
             tgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGcc                                                                                                                                >  3:241308/1‑150 (MQ=255)
             tgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGcc                                                                                                                                >  3:521204/1‑150 (MQ=255)
             tgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGcc                                                                                                                                >  4:118224/1‑150 (MQ=255)
             tgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGAGATAA‑CGGCGATGTTCGcc                                                                                                                                >  1:65483/1‑150 (MQ=255)
                  cAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGa                                                                                                                           >  3:745703/1‑150 (MQ=255)
                   aaTGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAg                                                                                                                          <  4:618430/150‑1 (MQ=255)
                    aTGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGc                                                                                                                         <  2:298473/150‑1 (MQ=255)
                       gcggcCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTa                                                                                                                      >  1:534398/1‑150 (MQ=255)
                        cggcCAGACCTGCAACACCCGCACCAATTGCTATCACTCTAGCCGGTGACACTTTCACGGACGCAGTAATTTGCCCGGTAAAAAAGCGCCCAATAAAATGTTCCTCTTCAACAATAGCGCGATAA‑CGGCGATGTTCGCCATCGAgataaa                                                                                                                     >  4:212957/1‑145 (MQ=255)
                         ggcCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                    <  2:1196786/150‑1 (MQ=255)
                            cAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCg                                                                                                                 >  3:473759/1‑150 (MQ=255)
                             aGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                    <  2:431176/146‑1 (MQ=255)
                             aGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                    >  1:431176/1‑146 (MQ=255)
                              gACCTGCACCACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                               <  2:93758/150‑1 (MQ=255)
                              gACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                               <  3:731718/150‑1 (MQ=255)
                               aCCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTcc                                                                                                              <  1:533759/150‑1 (MQ=255)
                                 cTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAg                                                                                                            <  3:141125/150‑1 (MQ=255)
                                   gCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCg                                                                                                          >  2:1058158/1‑150 (MQ=255)
                                     aaCACCCGCACAAATCACCTGCGCTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGAt                                                                                                        <  2:443403/150‑1 (MQ=255)
                                          ccGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTTCGTCCAGCGATTGTgc                                                                                                   >  2:333368/1‑150 (MQ=255)
                                          ccGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgc                                                                                                   >  2:928936/1‑150 (MQ=255)
                                            gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                 <  1:928936/150‑1 (MQ=255)
                                            gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                 <  4:246010/150‑1 (MQ=255)
                                            gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                 >  3:246010/1‑150 (MQ=255)
                                            gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCACAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgtgc                                                                                                 <  3:452550/150‑1 (MQ=255)
                                              cccAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                               <  1:1021668/149‑1 (MQ=255)
                                              aCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                               <  1:1150363/150‑1 (MQ=255)
                                              aCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                               <  4:521204/150‑1 (MQ=255)
                                              aCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                               <  1:25883/150‑1 (MQ=255)
                                                cAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTga                                                                                             <  1:136854/150‑1 (MQ=255)
                                                cAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTga                                                                                             <  2:85434/150‑1 (MQ=255)
                                                cAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTga                                                                                             <  1:137051/150‑1 (MQ=255)
                                                 aaTCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgag                                                                                            <  4:30442/150‑1 (MQ=255)
                                                  aTCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga                                                                                           <  4:500562/150‑1 (MQ=255)
                                                     aCCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAc                                                                                        >  2:947586/1‑150 (MQ=255)
                                                      ccATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg                                                                                       >  2:243831/1‑150 (MQ=255)
                                                        aTCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcgcg                                                                                     <  2:78549/150‑1 (MQ=255)
                                                        aTCACTGTTGCCGGTGGCACTTTCCCGGCCTCAGGAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGCTAA‑CGGCGACGTTCGCCATCGAGCTAAGTGCGTCCAGCGATTGCGCGCGTGAGATAcgcg                                                                                     >  2:998225/1‑150 (MQ=255)
                                                           aCTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGca                                                                                  <  1:606198/150‑1 (MQ=255)
                                                             ttttGCCGGTGGCACTATCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGAGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCAAGCAATTGTGCGCGTGAGATACGCGGcaca                                                                                >  2:1144021/1‑150 (MQ=255)
                                                              tttGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACag                                                                               >  1:939271/1‑150 (MQ=255)
                                                                tGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACagag                                                                             <  1:745127/150‑1 (MQ=255)
                                                                 gCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGt                                                                            <  2:1006640/150‑1 (MQ=255)
                                                                  ccGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTc                                                                           >  2:1096126/1‑150 (MQ=255)
                                                                     gTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCAt                                                                        >  4:318845/1‑150 (MQ=255)
                                                                        gCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGc                                                                     >  2:598508/1‑150 (MQ=255)
                                                                           cTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCAt                                                                  >  1:1074943/1‑150 (MQ=255)
                                                                           cTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCAt                                                                  >  2:337272/1‑150 (MQ=255)
                                                                           cTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTAATGTGCCGCTTCCCCCATGCCGCGGTAA‑CGGCAATGTTCCCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCAt                                                                  <  3:735208/150‑1 (MQ=255)
                                                                            tttCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCGGCTTCAACAATAGCGCGATAACCGGCGATGTTGGCCATCGAGCCTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCAt                                                                  <  4:436975/150‑1 (MQ=255)
                                                                                ccGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACgg                                                             >  1:1090845/1‑150 (MQ=255)
                                                                                 cGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCagtgacac                                                            >  2:34106/1‑143 (MQ=255)
                                                                                 cGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGt                                                            >  4:172082/1‑150 (MQ=255)
                                                                                         gTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgttac                                                    <  2:1074943/150‑1 (MQ=255)
                                                                                          tAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgttacg                                                   >  1:1133973/1‑150 (MQ=255)
                                                                                          tAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgttacg                                                   <  1:511414/150‑1 (MQ=255)
                                                                                                        aagaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGttttt                                     <  2:65483/150‑1 (MQ=255)
                                                                                                         agaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTg                                    <  3:300764/150‑1 (MQ=255)
                                                                                                         agaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTg                                    <  2:215331/150‑1 (MQ=255)
                                                                                                          gaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGc                                   <  3:331103/150‑1 (MQ=255)
                                                                                                           aagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCa                                  >  3:450563/1‑150 (MQ=255)
                                                                                                           aagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCa                                  >  4:527410/1‑150 (MQ=255)
                                                                                                                 ccAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGTGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAAtt                            >  4:402871/1‑150 (MQ=255)
                                                                                                                  cAAATTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTc                           >  3:158519/1‑150 (MQ=255)
                                                                                                                   aagttCATGTGCCGCTACAACAATGGCGCGATAA‑CGGCGATGTCCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACTTTCCGCAAGTTTTTGCATTAATTcc                          <  4:407793/147‑1 (MQ=255)
                                                                                                                   aaaTTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTcc                          >  2:487756/1‑150 (MQ=255)
                                                                                                                   aaaTTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTcc                          <  1:509081/150‑1 (MQ=255)
                                                                                                                   aaaTTCATGTGCCGCTTCAACAACTGCGCGATAA‑CGGCGAGGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAAATCCGCGGCTCAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTcc                          <  2:939271/150‑1 (MQ=255)
                                                                                                                    aaTTCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCg                         >  1:1006086/1‑150 (MQ=255)
                                                                                                                      ttCATGTGCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGa                       <  2:751528/150‑1 (MQ=255)
                                                                                                                          tgtgCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCt                   >  2:692632/1‑150 (MQ=255)
                                                                                                                           gtgCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGa                       <  1:247797/145‑1 (MQ=255)
                                                                                                                           gtgCCGCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGa                       >  2:247797/1‑145 (MQ=255)
                                                                                                                            tgCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTCTTTGCATTAATTCCGGAGTCTg                  >  2:1217202/1‑150 (MQ=255)
                                                                                                                                gCTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAg             >  4:80400/1‑150 (MQ=255)
                                                                                                                                 cTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCaa                                           >  1:73308/1‑119 (MQ=255)
                                                                                                                                 cTTCAACAATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCaa                                           <  2:73308/119‑1 (MQ=255)
                                                                                                                                   tCAACGATGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGcc          <  2:197677/150‑1 (MQ=255)
                                                                                                                                     aacaaTGGCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAg        >  3:624388/1‑150 (MQ=255)
                                                                                                                                           ggCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa  <  2:259133/150‑1 (MQ=255)
                                                                                                                                           ggCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa  <  1:781808/150‑1 (MQ=255)
                                                                                                                                           ggCGCGATAA‑CGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa  <  3:318845/150‑1 (MQ=255)
                                                                                                                                            gcgcGATAACCGGCGATGGTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa  <  2:230055/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
CGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAA  >  NC_000913_3_pae_pgi/1677342‑1677631

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