Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,666,978 A→G 5.5% L290S (TTG→TCG) ‡ hcaT ← putative 3‑phenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,666,9780AG 5.5% 536.6 / 15.6 164L290S (TTG→TCG) ‡hcaTputative 3‑phenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base A (73/81);  new base G (5/4);  total (78/86)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 171 total aligned reads displayed.

CGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATAC  >  NC_000913/2666829‑2667127
                                                                                                                                                     |                                                                                                                                                     
   tACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACg                                                                                                                                                    <  1:1565640/150‑1 (MQ=255)
     cTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATTTAACGCTTGGCGGCCAGGTGGCTGACCGTGATGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGcc                                                                                                                                                  <  1:1966989/150‑1 (MQ=255)
      tGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCg                                                                                                                                                 >  1:1431495/1‑150 (MQ=255)
        ccgccCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTa                                                                                                                                               >  1:1151034/1‑150 (MQ=255)
         cgccCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAg                                                                                                                                              >  2:2342955/1‑150 (MQ=255)
         cgccCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAg                                                                                                                                              <  1:475126/150‑1 (MQ=255)
         cgccCAGCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCACTATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAg                                                                                                                                              <  2:2322736/150‑1 (MQ=255)
            ccATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTg                                                                                                                                           <  2:1595106/150‑1 (MQ=255)
              aTCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTGCCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACTACGGCGTTGCCGTGGTTgat                                                                                                                                         <  1:126632/150‑3 (MQ=255)
              aTCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACTACGGCGTTGCCGTGGTTgat                                                                                                                                         <  2:1093677/150‑3 (MQ=255)
              aTCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCt                                                                                                                                         <  1:1059075/150‑1 (MQ=255)
                  cgACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCa                                                                                                                                     >  2:1982833/1‑150 (MQ=255)
                      ggCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAAt                                                                                                                                 >  2:1737925/1‑150 (MQ=255)
                      ggCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAAt                                                                                                                                 >  2:445392/1‑150 (MQ=255)
                        cAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGc                                                                                                                               <  2:123327/150‑1 (MQ=255)
                             tAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGACCCAg                                                                                                                          >  1:233718/1‑150 (MQ=255)
                                 ccGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCa                                                                                                                      <  1:497596/150‑1 (MQ=255)
                                      tGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCacgacg                                                                                                                 >  1:958998/1‑150 (MQ=255)
                                      tGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTAcg                                                                                                                 >  1:2292168/1‑150 (MQ=255)
                                      tGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTAcg                                                                                                                 >  1:2292245/1‑150 (MQ=255)
                                         aaaCGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTAcgccg                                                                                                              <  1:19603/150‑1 (MQ=255)
                                          aaCGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCACCGCCGTAGTTGCTCCCATCATGCCCCAGCGCACTAcgccgc                                                                                                             >  2:2243019/1‑150 (MQ=255)
                                           aCGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATGTGCACCACTATCAACCACGGCACCGCCGTAGGTGCTCCCATAATCCCCCAGCGCACTACGCCGCa                                                                                                            >  2:356949/1‑150 (MQ=255)
                                            cGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCaa                                                                                                           <  1:2115629/150‑1 (MQ=255)
                                               aTGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATc                                                                                                        >  2:1711124/1‑150 (MQ=255)
                                               aTGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATc                                                                                                        >  1:1212468/1‑150 (MQ=255)
                                                    ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGTCCCAGCGCACTACGCCGCAAATCGCCga                                                                                                   >  1:715361/1‑150 (MQ=255)
                                                    ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCga                                                                                                   >  1:36829/1‑150 (MQ=255)
                                                    ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCga                                                                                                   >  2:1961499/1‑150 (MQ=255)
                                                    ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCga                                                                                                   >  2:1373043/1‑150 (MQ=255)
                                                    ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCga                                                                                                   >  1:2100611/1‑150 (MQ=255)
                                                    ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCga                                                                                                   >  1:1952911/1‑150 (MQ=255)
                                                    ttCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCga                                                                                                   >  1:2289482/1‑150 (MQ=255)
                                                      cGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCgaga                                                                                                 <  1:1209770/150‑1 (MQ=255)
                                                      cGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCgaga                                                                                                 <  2:705424/150‑1 (MQ=255)
                                                       gCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGAt                                                                                                >  2:1156205/1‑150 (MQ=255)
                                                        cTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATc                                                                                               <  2:310372/150‑1 (MQ=255)
                                                           cccTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAAc                                                                                            >  1:2297522/1‑150 (MQ=255)
                                                           cccTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAAc                                                                                            >  1:972858/1‑150 (MQ=255)
                                                           cccTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAAc                                                                                            >  1:972828/1‑150 (MQ=255)
                                                                 cGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCata                                                                                      >  2:4887/1‑150 (MQ=255)
                                                                 cGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCata                                                                                      >  1:1518675/1‑150 (MQ=255)
                                                                     cagcaATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATcgc                                                                                  <  2:517172/150‑1 (MQ=255)
                                                                      agcaATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCa                                                                                                                           >  2:2113008/1‑108 (MQ=255)
                                                                      agcaATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATcgcg                                                                                 >  1:1933077/1‑150 (MQ=255)
                                                                       gcaATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGt                                                                                <  1:1864505/150‑1 (MQ=255)
                                                                          atatAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCa                                                                             <  2:1419452/150‑1 (MQ=255)
                                                                             taACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTa                                                                          <  1:1536969/150‑1 (MQ=255)
                                                                              aaCGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTAc                                                                         >  2:307011/1‑150 (MQ=255)
                                                                                 gCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAc                                                                      >  2:1126961/1‑150 (MQ=255)
                                                                                     ggcggcCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggcg                                                                  >  1:1840651/1‑150 (MQ=255)
                                                                                          ccAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGaaaa                                                             <  2:313423/150‑1 (MQ=255)
                                                                                             ggTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGtt                                                          >  1:1041394/1‑150 (MQ=255)
                                                                                             ggTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGtt                                                          <  2:218180/150‑1 (MQ=255)
                                                                                             ggTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGAAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGtt                                                          >  1:304845/1‑150 (MQ=255)
                                                                                              gTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggc                                                                   >  1:51619/1‑140 (MQ=255)
                                                                                               tGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTtta                                                        <  2:1884243/150‑1 (MQ=255)
                                                                                               tGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTtta                                                        <  2:1048933/150‑1 (MQ=255)
                                                                                                 gCAGACCGTGAAGGTTCCGCTATGCAGAATTTGCACCACTATCAACCACGGCAACGCAGTAGTTGCTCCCATAATGCCCCAGCGCACTAAGCCGCAAATCTCCGAGATCAACAGCATATCGCGTGCTCTACAACGGCGGAAAAGTttatt                                                      <  2:703674/150‑1 (MQ=255)
                                                                                                 gCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatt                                                      >  2:1269191/1‑150 (MQ=255)
                                                                                                 gCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatt                                                      <  2:650632/150‑1 (MQ=255)
                                                                                                 gCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatt                                                      >  2:346601/1‑150 (MQ=255)
                                                                                                  cAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatta                                                     <  2:1431495/150‑1 (MQ=255)
                                                                                                    gACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACt                                                   <  2:865605/150‑1 (MQ=255)
                                                                                                      ccGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCa                                                 <  2:2280669/150‑1 (MQ=255)
                                                                                                         tGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAgcg                                              >  1:2185162/1‑150 (MQ=255)
                                                                                                                  cGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa                                     <  2:1840651/150‑1 (MQ=255)
                                                                                                                  cGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa                                     <  2:1695452/150‑1 (MQ=255)
                                                                                                                  cGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa                                     <  2:1695497/150‑1 (MQ=255)
                                                                                                                  cGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa                                     <  1:2001170/150‑1 (MQ=255)
                                                                                                                  cGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa                                     >  1:23759/1‑150 (MQ=255)
                                                                                                                  cGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa                                     <  1:2131891/150‑1 (MQ=255)
                                                                                                                   gCAATGCAGAATTTGCGCCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAAt                                    <  1:815484/150‑1 (MQ=255)
                                                                                                                   gCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAAt                                    <  1:1283061/150‑1 (MQ=255)
                                                                                                                   gCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAAt                                    <  2:682767/150‑1 (MQ=255)
                                                                                                                   gCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAAt                                    <  1:1283072/150‑1 (MQ=255)
                                                                                                                   gCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAAt                                    >  1:2074938/1‑150 (MQ=255)
                                                                                                                        gCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACtt                               <  1:1625348/150‑1 (MQ=255)
                                                                                                                            aaTTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGc                           >  1:1343128/1‑150 (MQ=255)
                                                                                                                               ttGCACCACTATCAACCACGGCACCGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGccac                        <  2:2400672/150‑1 (MQ=255)
                                                                                                                               ttGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGccac                        <  2:2335826/150‑1 (MQ=255)
                                                                                                                               ttGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGccac                        <  2:262077/150‑1 (MQ=255)
                                                                                                                                     cacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGcc                  >  1:2051678/1‑150 (MQ=255)
                                                                                                                                        tATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAg               >  1:431479/1‑150 (MQ=255)
                                                                                                                                        tATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTAGTCCTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAg               >  2:784245/1‑150 (MQ=255)
                                                                                                                                        gagcAACTACGGCGTTGCCGTGGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCCCTCCAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAg               >  1:538174/4‑150 (MQ=255)
                                                                                                                                           cAACTACGGCGTTGCCGTGGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGa            >  2:2081986/1‑150 (MQ=255)
                                                                                                                                           cAACTACGGCGTTGCCGTGGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGa            >  2:510672/1‑150 (MQ=255)
                                                                                                                                           cAACTACGGCGTTGCCGTGGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGa            >  2:605891/1‑150 (MQ=255)
                                                                                                                                           cAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGa            <  1:1737925/150‑1 (MQ=255)
                                                                                                                                           cAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGa            <  1:445392/150‑1 (MQ=255)
                                                                                                                                                aCGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACcaca       <  2:2289482/150‑1 (MQ=255)
                                                                                                                                                aCGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACcaca       >  1:4753/1‑150 (MQ=255)
                                                                                                                                                aCGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACcaca       >  2:780109/1‑150 (MQ=255)
                                                                                                                                                  ggCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACaaa     <  2:1635153/150‑1 (MQ=255)
                                                                                                                                                   gCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCTCTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAAt    <  2:1827975/150‑1 (MQ=255)
                                                                                                                                                   gCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAAt    <  2:1261844/150‑1 (MQ=255)
                                                                                                                                                     aaCGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATAc  <  1:1602547/150‑1 (MQ=255)
                                                                                                                                                     aaCGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATAc  <  2:2295506/150‑1 (MQ=255)
                                                                                                                                                     aaCGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATAc  <  1:661505/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGATACTGCCGCCCATCGCGACGGCAGAGTAAACCGCCTGTAAACGGATGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATAC  >  NC_000913/2666829‑2667127

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