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)
      tGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGAt                                                                                                                                       >  4:300764/1‑150 (MQ=255)
           gctgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCg                                                                                                                                  <  1:389927/150‑1 (MQ=255)
             tgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGcc                                                                                                                                >  3:241308/1‑150 (MQ=255)
             tgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGcc                                                                                                                                >  3:521204/1‑150 (MQ=255)
             tgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGcc                                                                                                                                >  4:118224/1‑150 (MQ=255)
             tgcGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGAGATAACGGCGATGTTCGcc                                                                                                                                >  1:65483/1‑150 (MQ=255)
                  cAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGa                                                                                                                           >  3:745703/1‑150 (MQ=255)
                   aaTGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAg                                                                                                                          <  4:618430/150‑1 (MQ=255)
                    aTGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGc                                                                                                                         <  2:298473/150‑1 (MQ=255)
                       gcggcCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTa                                                                                                                      >  1:534398/1‑150 (MQ=255)
                        cggcCAGACCTGCAACACCCGCACCAATTGCTATCACTCTAGCCGGTGACACTTTCACGGACGCAGTAATTTGCCCGGTAAAAAAGCGCCCAATAAAATGTTCCTCTTCAACAATAGCGCGATAACGGCGATGTTCGCCATCGAgataaa                                                                                                                     >  4:212957/1‑145 (MQ=255)
                         ggcCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                    <  2:1196786/150‑1 (MQ=255)
                            cAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCg                                                                                                                 >  3:473759/1‑150 (MQ=255)
                             aGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                    <  2:431176/146‑1 (MQ=255)
                             aGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGt                                                                                                                    >  1:431176/1‑146 (MQ=255)
                              gACCTGCACCACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                               <  2:93758/150‑1 (MQ=255)
                              gACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTc                                                                                                               <  3:731718/150‑1 (MQ=255)
                               aCCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTcc                                                                                                              <  1:533759/150‑1 (MQ=255)
                                 cTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAg                                                                                                            <  3:141125/150‑1 (MQ=255)
                                   gCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCg                                                                                                          >  2:1058158/1‑150 (MQ=255)
                                     aaCACCCGCACAAATCACCTGCGCTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGAt                                                                                                        <  2:443403/150‑1 (MQ=255)
                                          ccGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTTCGTCCAGCGATTGTgc                                                                                                   >  2:333368/1‑150 (MQ=255)
                                          ccGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgc                                                                                                   >  2:928936/1‑150 (MQ=255)
                                            gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                 <  1:928936/150‑1 (MQ=255)
                                            gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                 <  4:246010/150‑1 (MQ=255)
                                            gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgcgc                                                                                                 >  3:246010/1‑150 (MQ=255)
                                            gCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCACAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTgtgc                                                                                                 <  3:452550/150‑1 (MQ=255)
                                              cccAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                               <  1:1021668/149‑1 (MQ=255)
                                              aCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                               <  1:1150363/150‑1 (MQ=255)
                                              aCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                               <  4:521204/150‑1 (MQ=255)
                                              aCCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGt                                                                                               <  1:25883/150‑1 (MQ=255)
                                                cAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTga                                                                                             <  1:136854/150‑1 (MQ=255)
                                                cAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTga                                                                                             <  2:85434/150‑1 (MQ=255)
                                                cAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTga                                                                                             <  1:137051/150‑1 (MQ=255)
                                                 aaTCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgag                                                                                            <  4:30442/150‑1 (MQ=255)
                                                  aTCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTgaga                                                                                           <  4:500562/150‑1 (MQ=255)
                                                     aCCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAc                                                                                        >  2:947586/1‑150 (MQ=255)
                                                      ccATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcg                                                                                       >  2:243831/1‑150 (MQ=255)
                                                        aTCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATAcgcg                                                                                     <  2:78549/150‑1 (MQ=255)
                                                        aTCACTGTTGCCGGTGGCACTTTCCCGGCCTCAGGAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGCTAACGGCGACGTTCGCCATCGAGCTAAGTGCGTCCAGCGATTGCGCGCGTGAGATAcgcg                                                                                     >  2:998225/1‑150 (MQ=255)
                                                           aCTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGca                                                                                  <  1:606198/150‑1 (MQ=255)
                                                             ttttGCCGGTGGCACTATCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGAGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCAAGCAATTGTGCGCGTGAGATACGCGGcaca                                                                                >  2:1144021/1‑150 (MQ=255)
                                                              tttGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACag                                                                               >  1:939271/1‑150 (MQ=255)
                                                                tGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACagag                                                                             <  1:745127/150‑1 (MQ=255)
                                                                 gCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGt                                                                            <  2:1006640/150‑1 (MQ=255)
                                                                  ccGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTc                                                                           >  2:1096126/1‑150 (MQ=255)
                                                                     gTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCAt                                                                        >  4:318845/1‑150 (MQ=255)
                                                                        gCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGc                                                                     >  2:598508/1‑150 (MQ=255)
                                                                           cTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCAt                                                                  >  1:1074943/1‑150 (MQ=255)
                                                                           cTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCAt                                                                  >  2:337272/1‑150 (MQ=255)
                                                                           cTTTCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTAATGTGCCGCTTCCCCCATGCCGCGGTAACGGCAATGTTCCCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCAt                                                                  <  3:735208/150‑1 (MQ=255)
                                                                            tttCCCGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCGGCTTCAACAATAGCGCGATAACCGGCGATGTTGGCCATCGAGCCTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCAt                                                                  <  4:436975/150‑1 (MQ=255)
                                                                                ccGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACgg                                                             >  1:1090845/1‑150 (MQ=255)
                                                                                 cGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCagtgacac                                                            >  2:34106/1‑143 (MQ=255)
                                                                                 cGGCCGCAGTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGt                                                            >  4:172082/1‑150 (MQ=255)
                                                                                         gTAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgttac                                                    <  2:1074943/150‑1 (MQ=255)
                                                                                          tAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgttacg                                                   >  1:1133973/1‑150 (MQ=255)
                                                                                          tAATTTGCCCGGTAAAGAAGCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCacgttacg                                                   <  1:511414/150‑1 (MQ=255)
                                                                                                        aagaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGttttt                                     <  2:65483/150‑1 (MQ=255)
                                                                                                         agaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTg                                    <  3:300764/150‑1 (MQ=255)
                                                                                                         agaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTg                                    <  2:215331/150‑1 (MQ=255)
                                                                                                          gaagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGc                                   <  3:331103/150‑1 (MQ=255)
                                                                                                           aagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCa                                  >  3:450563/1‑150 (MQ=255)
                                                                                                           aagCGCCCAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCa                                  >  4:527410/1‑150 (MQ=255)
                                                                                                                 ccAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGTGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAAtt                            >  4:402871/1‑150 (MQ=255)
                                                                                                                  cAAATTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTc                           >  3:158519/1‑150 (MQ=255)
                                                                                                                   aagttCATGTGCCGCTACAACAATGGCGCGATAACGGCGATGTCCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACTTTCCGCAAGTTTTTGCATTAATTcc                          <  4:407793/147‑1 (MQ=255)
                                                                                                                   aaaTTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTcc                          >  2:487756/1‑150 (MQ=255)
                                                                                                                   aaaTTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTcc                          <  1:509081/150‑1 (MQ=255)
                                                                                                                   aaaTTCATGTGCCGCTTCAACAACTGCGCGATAACGGCGAGGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAAATCCGCGGCTCAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTcc                          <  2:939271/150‑1 (MQ=255)
                                                                                                                    aaTTCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCg                         >  1:1006086/1‑150 (MQ=255)
                                                                                                                      ttCATGTGCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGa                       <  2:751528/150‑1 (MQ=255)
                                                                                                                          tgtgCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCt                   >  2:692632/1‑150 (MQ=255)
                                                                                                                           gtgCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGa                       <  1:247797/145‑1 (MQ=255)
                                                                                                                           gtgCCGCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGa                       >  2:247797/1‑145 (MQ=255)
                                                                                                                            tgCCGCTTCAACAATGGCGCGATAACCGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTCTTTGCATTAATTCCGGAGTCTg                  >  2:1217202/1‑150 (MQ=255)
                                                                                                                                gCTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAg             >  4:80400/1‑150 (MQ=255)
                                                                                                                                 cTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCaa                                           >  1:73308/1‑119 (MQ=255)
                                                                                                                                 cTTCAACAATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCaa                                           <  2:73308/119‑1 (MQ=255)
                                                                                                                                   tCAACGATGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGcc          <  2:197677/150‑1 (MQ=255)
                                                                                                                                     aacaaTGGCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAg        >  3:624388/1‑150 (MQ=255)
                                                                                                                                           ggCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa  <  2:259133/150‑1 (MQ=255)
                                                                                                                                           ggCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa  <  1:781808/150‑1 (MQ=255)
                                                                                                                                           ggCGCGATAACGGCGATGTTCGCCATCGAGCTTAGTGCGTCCAGCGATTGTGCGCGTGAGATACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATaaaa  <  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: