Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,738,088 T→A 26.6% F134Y (TTC→TAC)  yiaM → 2,3‑diketo‑L‑gulonate TRAP transporter small permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,738,0880TA26.6% 621.6 / 270.9 357F134Y (TTC→TAC) yiaM2,3‑diketo‑L‑gulonate TRAP transporter small permease protein
Reads supporting (aligned to +/- strand):  ref base T (116/146);  new base A (52/43);  total (168/189)
Fisher's exact test for biased strand distribution p-value = 9.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 366 total aligned reads displayed.

CCATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTT  >  CP009273/3737939‑3738237
                                                                                                                                                     |                                                                                                                                                     
ccATTCCTTAATCTTATTTATTTGTTGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCTGCCTGATGTATGCCGCTTGTCTGCCCACCAGCATCGTTATCTCAtt                                                                                                                                                       >  1:724414/1‑150 (MQ=255)
ccATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCAtt                                                                                                                                                       >  1:4277161/1‑150 (MQ=255)
ccATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCAtt                                                                                                                                                       >  1:585516/1‑150 (MQ=255)
 cATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTc                                                                                                                                                      >  2:4939763/1‑150 (MQ=255)
  aTTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCt                                                                                                                                                     >  1:1061945/1‑150 (MQ=255)
   ttCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCtt                                                                                                                                                    <  2:484662/150‑1 (MQ=255)
   ttCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCACTTGTCTGCCCACCAGCCTCGTTATCGCATTCtt                                                                                                                                                    <  1:169183/150‑1 (MQ=255)
   ttCCTGAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATAATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTGATCGCATTCtt                                                                                                                                                    >  2:2132129/1‑150 (MQ=255)
    tCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACttt                                                                                                                                                   >  2:389566/1‑150 (MQ=255)
     ccTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTg                                                                                                                                                  >  1:2861232/1‑150 (MQ=255)
      cTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGa                                                                                                                                                 >  2:3483276/1‑150 (MQ=255)
        tAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAAc                                                                                                                                               <  1:1473122/150‑1 (MQ=255)
          ctCTTATTTATGTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAAGTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTg                                                                                                                                             <  2:1056828/149‑1 (MQ=255)
          aTCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTg                                                                                                                                             >  1:1820140/1‑150 (MQ=255)
                tttatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCAt                                                                                                                                       >  1:2096244/1‑150 (MQ=255)
                 ttatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCAtt                                                                                                                                      >  2:4436531/1‑150 (MQ=255)
                   atttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTta                                                                                                                                    >  1:5107183/1‑150 (MQ=255)
                   atttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTta                                                                                                                                    >  1:4202701/1‑150 (MQ=255)
                      tgtgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTatat                                                                                                                                 <  1:1132921/150‑1 (MQ=255)
                      tgtgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTatat                                                                                                                                 <  1:5203714/150‑1 (MQ=255)
                        tgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCa                                                                                                                               >  2:3591215/1‑150 (MQ=255)
                        tgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCa                                                                                                                               <  2:409125/150‑1 (MQ=255)
                          gagcACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAAc                                                                                                                             <  1:445054/148‑1 (MQ=255)
                           cgcACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACt                                                                                                                            <  1:3024206/150‑1 (MQ=255)
                            gcACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtc                                                                                                                           <  1:2574740/150‑1 (MQ=255)
                             cACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtca                                                                                                                          <  1:3134083/150‑1 (MQ=255)
                             cACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtca                                                                                                                          <  2:2022570/150‑1 (MQ=255)
                                tGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTTTGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtcatca                                                                                                                       <  1:2821280/150‑1 (MQ=255)
                                tGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtcatca                                                                                                                       <  2:4251389/150‑1 (MQ=255)
                                  gCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCAcg                                                                                                                     <  1:4666078/150‑1 (MQ=255)
                                   ccTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCAcgc                                                                                                                    >  2:3055465/1‑150 (MQ=255)
                                    cTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGAGTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCAcgcg                                                                                                                   >  2:273319/1‑150 (MQ=255)
                                     tGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGt                                                                                                                  >  2:3228428/1‑150 (MQ=255)
                                             aaCGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtc                                                                                                          >  2:4568447/1‑150 (MQ=255)
                                             aaCGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtc                                                                                                          >  1:2065396/1‑150 (MQ=255)
                                              aCGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtct                                                                                                         <  1:3783262/150‑1 (MQ=255)
                                              aCGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtct                                                                                                         <  1:687985/150‑1 (MQ=255)
                                               cGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtctc                                                                                                        <  1:422526/150‑1 (MQ=255)
                                                  tAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTaa                                                                                                     <  1:1781992/150‑1 (MQ=255)
                                                  tAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTaa                                                                                                     >  2:838615/1‑150 (MQ=255)
                                                   aaaaaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAAc                                                                                                    >  1:4773582/1‑150 (MQ=255)
                                                   aaaaaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAAc                                                                                                    <  1:4854174/150‑1 (MQ=255)
                                                     aaaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGt                                                                                                  >  1:3579121/1‑150 (MQ=255)
                                                        cAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGc                                                                                               >  1:454299/1‑150 (MQ=255)
                                                            tCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCAcc                                                                                           >  2:3809678/1‑150 (MQ=255)
                                                                 aCTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAg                                                                                      <  1:906654/150‑1 (MQ=255)
                                                                   tGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGa                                                                                    >  2:589570/1‑150 (MQ=255)
                                                                    ggAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGTTGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAg                                                                                   <  2:1487102/150‑1 (MQ=255)
                                                                    ggAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAg                                                                                   <  2:1718745/150‑1 (MQ=255)
                                                                        cGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTa                                                                               >  2:4150555/1‑150 (MQ=255)
                                                                         gATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTaa                                                                              >  1:2917997/1‑150 (MQ=255)
                                                                         gATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTaa                                                                              >  2:3393125/1‑150 (MQ=255)
                                                                         gATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTaa                                                                              >  1:4076003/1‑150 (MQ=255)
                                                                          attattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAAt                                                                             <  2:2981586/150‑1 (MQ=255)
                                                                             attCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCAt                                                                          <  2:1309104/150‑1 (MQ=255)
                                                                             attCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCAt                                                                          >  2:4815709/1‑150 (MQ=255)
                                                                               tCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATgg                                                                        >  1:2533922/1‑150 (MQ=255)
                                                                                cACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGc                                                                       >  2:4043236/1‑150 (MQ=255)
                                                                                cACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAAAGTCGCCACCGCAAGGAGCTTAATCATGGc                                                                       >  2:3329181/1‑150 (MQ=255)
                                                                                   cGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgt                                                                    >  1:473963/1‑150 (MQ=255)
                                                                                    gATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgtg                                                                   >  2:5315487/1‑150 (MQ=255)
                                                                                       ccTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTg                                                                >  1:839416/1‑150 (MQ=255)
                                                                                       ccTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTg                                                                >  2:1847545/1‑150 (MQ=255)
                                                                                       ccTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTg                                                                <  1:815685/150‑1 (MQ=255)
                                                                                        cTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCACGGAGCTTAATCATGGCTGTGCTGa                                                               >  1:4606774/1‑150 (MQ=255)
                                                                                         tCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAt                                                              <  1:3968889/150‑1 (MQ=255)
                                                                                         tCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAt                                                              <  2:4467017/150‑1 (MQ=255)
                                                                                          cGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAtt                                                             >  2:2321203/1‑150 (MQ=255)
                                                                                          cGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAtt                                                             <  2:3118414/150‑1 (MQ=255)
                                                                                             tttACCCATCGGCCTGATTTATGCAGCTTGTCTGCCCAACAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAttttt                                                          <  1:3692884/150‑1 (MQ=255)
                                                                                                aCCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTg                                                       >  2:609474/1‑150 (MQ=255)
                                                                                                 cccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTgg                                                      <  2:5090243/150‑1 (MQ=255)
                                                                                                   cATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGc                                                    >  1:1500331/1‑150 (MQ=255)
                                                                                                      cGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgt                                                 <  2:1344106/150‑1 (MQ=255)
                                                                                                       ggCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATATTTCTGGGctgtc                                                >  1:1924222/1‑150 (MQ=255)
                                                                                                        gCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtct                                               <  2:1818938/150‑1 (MQ=255)
                                                                                                         ccTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtctg                                              <  2:888707/150‑1 (MQ=255)
                                                                                                         ccTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtctg                                              <  2:5202417/150‑1 (MQ=255)
                                                                                                                tATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGggtg                                       >  1:944745/1‑150 (MQ=255)
                                                                                                                tATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCAATTATATCAACTCATCACGCGTAGCAATTCACTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTTggtg                                       >  2:5146183/1‑150 (MQ=255)
                                                                                                                  tGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGggtggt                                     >  1:2689013/1‑150 (MQ=255)
                                                                                                                      gCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTg                                 <  2:687811/150‑1 (MQ=255)
                                                                                                                      gCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTg                                 >  2:1777609/1‑150 (MQ=255)
                                                                                                                      gCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTg                                 <  1:2134803/150‑1 (MQ=255)
                                                                                                                        ttGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCt                               <  1:389566/150‑1 (MQ=255)
                                                                                                                               cccACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGtt                        >  2:805983/1‑150 (MQ=255)
                                                                                                                                     aGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTa                  >  1:320569/1‑150 (MQ=255)
                                                                                                                                     aGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTa                  >  1:1633192/1‑150 (MQ=255)
                                                                                                                                     aGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTa                  >  2:4287342/1‑150 (MQ=255)
                                                                                                                                       ccTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTAtt                <  2:2446501/150‑1 (MQ=255)
                                                                                                                                          cGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGcc             >  1:2216452/1‑150 (MQ=255)
                                                                                                                                          cGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGcc             <  1:2283067/150‑1 (MQ=255)
                                                                                                                                           gTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCt            <  2:4391860/150‑1 (MQ=255)
                                                                                                                                           gTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCt            >  1:1221035/1‑150 (MQ=255)
                                                                                                                                            ttATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTg           <  1:3992622/150‑1 (MQ=255)
                                                                                                                                                cGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCa       >  1:3241746/1‑150 (MQ=255)
                                                                                                                                                cGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCa       >  2:4032118/1‑150 (MQ=255)
                                                                                                                                                  cATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACt     >  1:3496083/1‑150 (MQ=255)
                                                                                                                                                  cATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACt     >  2:4168956/1‑150 (MQ=255)
                                                                                                                                                     tCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACtgtt  >  1:1342606/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTT  >  CP009273/3737939‑3738237

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: