Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,600,721 C→A 8.5% S149* (TCG→TAG)  hyfF → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,600,7210CA 8.5% 1136.0 / 53.5 301S149* (TCG→TAG) hyfFhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (156/117);  new base A (16/11);  total (172/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 305 total aligned reads displayed.

CATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGTATACGCCAACGCCGCCAGCGTTATGCCGCAGGCAGAA  >  CP009273/2600573‑2600869
                                                                                                                                                    |                                                                                                                                                    
 aTGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCg                                                                                                                                                     <  1:2932472/149‑1 (MQ=255)
 aTGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCg                                                                                                                                                     >  1:64801/1‑149 (MQ=255)
  tGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCGGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCgg                                                                                                                                                    >  2:905917/1‑149 (MQ=255)
   gAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGc                                                                                                                                                   >  2:2981105/1‑149 (MQ=255)
      gTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGtt                                                                                                                                                >  1:3722089/1‑149 (MQ=255)
      gTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGtt                                                                                                                                                <  1:1584812/149‑1 (MQ=255)
      gTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGtt                                                                                                                                                >  2:625219/1‑149 (MQ=255)
        ggCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTc                                                                                                                                              >  1:376940/1‑149 (MQ=255)
          cacaCGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTAGGCGTTTCTg                                                                                                                                            >  1:1220348/1‑149 (MQ=255)
           acacGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTgg                                                                                                                                           >  2:3700051/1‑149 (MQ=255)
                gCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGc                                                                                                                                      >  1:5027406/1‑149 (MQ=255)
                gCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGc                                                                                                                                      >  2:1479247/1‑149 (MQ=255)
                 cGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCa                                                                                                                                     <  1:468524/149‑1 (MQ=255)
                   aGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCAtt                                                                                                                                   >  2:1077967/1‑149 (MQ=255)
                   aGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCAtt                                                                                                                                   >  1:2955581/1‑149 (MQ=255)
                    gCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGAGTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCAttt                                                                                                                                  >  1:526080/1‑149 (MQ=255)
                     cTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTa                                                                                                                                 <  1:215743/149‑1 (MQ=255)
                     cTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTa                                                                                                                                 >  2:5201889/1‑149 (MQ=255)
                      tttCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTAGGCGTTTCTGGTAGGCATTTAc                                                                                                                                >  2:2288222/1‑149 (MQ=255)
                       ttCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACg                                                                                                                               >  2:4823216/1‑149 (MQ=255)
                        tCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACgg                                                                                                                              <  1:1200943/149‑1 (MQ=255)
                         cACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGt                                                                                                                             >  2:1501643/1‑149 (MQ=255)
                           cccGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCa                                                                                                                           >  1:1709370/1‑149 (MQ=255)
                           cccGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCa                                                                                                                           >  1:4229340/1‑149 (MQ=255)
                            ccGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAg                                                                                                                          >  1:114162/1‑149 (MQ=255)
                               ttACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGt                                                                                                                       >  2:3257297/1‑149 (MQ=255)
                               ttACGCTGTGCGAGTACTACGGTTTCTTCCATCTGTTTTTGGTCACCATGCTGCTGGTTGGGACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTAGGCGTTTCTGGTAGGCATTTACGGTCAGCGt                                                                                                                       <  2:7068/149‑1 (MQ=255)
                                tACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGtt                                                                                                                      >  2:130254/1‑149 (MQ=255)
                                 aCGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTc                                                                                                                     >  1:1653638/1‑149 (MQ=255)
                                 aCGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTGTCTGGTAGGCATTTACGGTCAGCGTTc                                                                                                                     >  2:2503676/1‑149 (MQ=255)
                                   gCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCAt                                                                                                                   <  1:4823216/149‑1 (MQ=255)
                                       tgCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCt                                                                                                               <  2:736544/149‑1 (MQ=255)
                                        gCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTg                                                                                                              <  1:720949/149‑1 (MQ=255)
                                           aTTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGaa                                                                                                           >  1:3898884/1‑149 (MQ=255)
                                             tactacGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCCTTTACGGTCAGCGTTCATCGCTGGAAGc                                                                                                         >  2:2299879/1‑149 (MQ=255)
                                                tacGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGc                                                                                                      >  1:2597613/1‑149 (MQ=255)
                                                tacGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGc                                                                                                      >  1:151804/1‑149 (MQ=255)
                                                  cGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTAGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCAt                                                                                                    <  1:3159101/149‑1 (MQ=255)
                                                   ggTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCTTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATg                                                                                                   >  2:3999047/1‑149 (MQ=255)
                                                   ggTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATg                                                                                                   <  1:1498171/149‑1 (MQ=255)
                                                    gTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATgg                                                                                                  <  2:1288569/149‑1 (MQ=255)
                                                    gTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTAAGGTCAGCGTTCATCGCTGGAAGCTGCATgg                                                                                                  >  1:1317038/1‑149 (MQ=255)
                                                     tttCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGa                                                                                                 <  1:2470098/149‑1 (MQ=255)
                                                      ttcttcCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGaa                                                                                                >  1:2376605/1‑149 (MQ=255)
                                                      ttcttcCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGaa                                                                                                >  1:4745164/1‑149 (MQ=255)
                                                      ttcttcCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGaa                                                                                                >  1:4115706/1‑149 (MQ=255)
                                                       tcttcCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAg                                                                                               <  1:2253333/149‑1 (MQ=255)
                                                            cATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTAGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTAcat                                                                                          >  2:1982280/1‑149 (MQ=255)
                                                             aTCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTAcatc                                                                                         >  2:175897/1‑149 (MQ=255)
                                                               cTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTAcatcat                                                                                       <  2:793818/149‑1 (MQ=255)
                                                                  tttttGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCattat                                                                                    <  2:4099041/149‑1 (MQ=255)
                                                                     ttGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTg                                                                                 >  2:3636264/1‑149 (MQ=255)
                                                                     ttGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTg                                                                                 <  2:4135188/149‑1 (MQ=255)
                                                                       gTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTa                                                                               <  2:5127169/149‑1 (MQ=255)
                                                                        ttCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTAc                                                                              <  2:314991/149‑1 (MQ=255)
                                                                         tCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACt                                                                             >  2:1840709/1‑149 (MQ=255)
                                                                           aCCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACCTCATTATTTGTACTGt                                                                           >  2:1179848/1‑149 (MQ=255)
                                                                           aCCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGt                                                                           >  2:4181817/1‑149 (MQ=255)
                                                                           aCCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTGTGTGGTAGGCATTTACGGTCAGCGGTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGt                                                                           >  1:2191935/1‑149 (MQ=255)
                                                                           aCCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTAGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGt                                                                           >  1:3958271/1‑149 (MQ=255)
                                                                              aTGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTgg                                                                        <  2:3718257/149‑1 (MQ=255)
                                                                              aTGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTgg                                                                        <  2:3063480/149‑1 (MQ=255)
                                                                               tgctgcTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGgt                                                                       <  2:4163635/149‑1 (MQ=255)
                                                                               tgctgcTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTAGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGgt                                                                       >  1:3928480/1‑149 (MQ=255)
                                                                                gctgctGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGgtg                                                                      <  1:1363098/149‑1 (MQ=255)
                                                                                  tgctgGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTc                                                                    >  2:3525118/1‑149 (MQ=255)
                                                                                      ggTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCtt                                                                <  2:1418169/149‑1 (MQ=255)
                                                                                      ggTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCtt                                                                <  1:457682/149‑1 (MQ=255)
                                                                                      ggTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTGTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGcgt                                                                >  1:4958603/1‑147 (MQ=255)
                                                                                      ggTTGTGACCATCAATAAGCTGATTGTGATGTGGGCGGCGATCGATGCCTACTCCTTATGCTGGGAGTTTCTGGTTGGCTTTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTTTAATGTTGGTGTCGCtt                                                                <  2:3707568/149‑1 (MQ=255)
                                                                                       gttgttACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCttt                                                               <  2:2619922/149‑1 (MQ=255)
                                                                                        ttgttACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCtttt                                                              >  2:3395819/1‑149 (MQ=255)
                                                                                        ttgttACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCgttt                                                              >  2:1900448/1‑149 (MQ=255)
                                                                                          gttACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTAGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTgg                                                            <  1:682588/149‑1 (MQ=255)
                                                                                           ttACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGt                                                           <  1:2071694/149‑1 (MQ=255)
                                                                                            tACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCGTTTTGTc                                                          >  1:1927545/1‑149 (MQ=255)
                                                                                               cAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGt                                                       <  1:2542948/149‑1 (MQ=255)
                                                                                                aGCAATAACCTGAGTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGtt                                                      <  1:3504244/149‑1 (MQ=255)
                                                                                                 gCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTc                                                     >  2:847460/1‑149 (MQ=255)
                                                                                                  cAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCg                                                    >  2:1455106/1‑149 (MQ=255)
                                                                                                    aTAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGt                                                  <  2:4244095/149‑1 (MQ=255)
                                                                                                       aCCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTAcc                                               <  2:4624656/149‑1 (MQ=255)
                                                                                                         cTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTAATGTTGGTGTCGCTTTTTGTCTGTTAGGTAACGt                                             >  2:4414120/1‑149 (MQ=255)
                                                                                                           gATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGc                                           >  2:101663/1‑149 (MQ=255)
                                                                                                             ttGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTg                                         >  2:2429008/1‑149 (MQ=255)
                                                                                                               gtgATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTACATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGt                                       <  2:1600763/149‑1 (MQ=255)
                                                                                                                tgATGTGGGCGGCGATCGTAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGta                                      <  2:3539613/149‑1 (MQ=255)
                                                                                                                   tgtgGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCGTTTGGTCGGTTCGGTACCGTGCTGGTATAc                                   >  2:2112268/1‑149 (MQ=255)
                                                                                                                       ggcggcGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGTATACGCCa                               <  2:4937865/149‑1 (MQ=255)
                                                                                                                        gcggcgATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGTATACGCCaa                              <  2:2843735/149‑1 (MQ=255)
                                                                                                                        gcggcgATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGTATACGCCaa                              <  2:2399848/149‑1 (MQ=255)
                                                                                                                                    gacaccagCTTATGCTCGGCGTTTCTTGTAGTCATTTACGGTCAGGGTTGATCGCTGGTAGCTGCTTGGTAGTACTTAATTTTTTGTACTGTTGGTGTTGCTTTTGGTCTGTTCGGTACCGTGCTGGTATTAGCCTAAGCCGCCATCGt                  <  2:2587655/147‑1 (MQ=255)
                                                                                                                                      caccacCTTAAGCTAGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGTATACGCCAACGCCGCCAGCGTTa                <  2:1220348/149‑1 (MQ=255)
                                                                                                                                       accaccTTAAGCTAGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTGTGTTCGGTAACGGGCTGGTATAAGCCAACGCCGCCACCGGTAt               >  1:101269/1‑149 (MQ=255)
                                                                                                                                        ccaccTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGTATACGCCAACGCCGCCAGCGTTATg              <  2:4757077/149‑1 (MQ=255)
                                                                                                                                         caccTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGTATACGCCAACGCCGCCAGCGTTATGc             <  2:2194183/149‑1 (MQ=255)
                                                                                                                                          accTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCATTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTTGTCTGTTCGGTACCGTGCTGGTATACGCCAACGCCCCCACCGTTATTcc            >  1:855260/1‑149 (MQ=255)
                                                                                                                                            cTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCCTGCGGGTATAAGCCAACCCCCCCACCGGTAttaccc          >  1:5120826/1‑144 (MQ=255)
                                                                                                                                              tAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGTATACGCCAACGCCGCCAGCGTTATGCCgcag        <  1:412595/149‑1 (MQ=255)
                                                                                                                                                 gCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGTATACGCCAACGCCGCCAGCGTTATGCCgcaggca     <  1:925739/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGCTGCATGGAAGTACATCATTATTTGTACTGTTGGTGTCGCTTTTGGTCTGTTCGGTACCGTGCTGGTATACGCCAACGCCGCCAGCGTTATGCCGCAGGCAGAA  >  CP009273/2600573‑2600869

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