Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,933,680 G→T 7.2% S98* (TCG→TAG)  hsrA ← putative multidrug or homocysteine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,933,6800GT 7.2% 964.5 / 30.9 250S98* (TCG→TAG) hsrAputative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base G (124/108);  new base T (9/9);  total (133/117)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 253 total aligned reads displayed.

GCCCCACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGC  >  CP009273/3933531‑3933829
                                                                                                                                                     |                                                                                                                                                     
   ccaccaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAAt                                                                                                                                                    >  2:2072686/1‑150 (MQ=255)
     accaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTa                                                                                                                                                  <  2:2784505/150‑1 (MQ=255)
      ccaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAg                                                                                                                                                 <  2:3895957/150‑1 (MQ=255)
       caGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACGACCAGCTGTGGTAGCGAATTAGa                                                                                                                                                >  2:3465067/1‑150 (MQ=255)
          aCCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAg                                                                                                                                             <  1:2608994/150‑1 (MQ=255)
           cccGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGt                                                                                                                                            >  2:3487019/1‑150 (MQ=255)
           cccGGCATGGCGACAAAAGTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGt                                                                                                                                            <  2:853862/150‑1 (MQ=255)
                 aTGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAg                                                                                                                                      >  2:3424555/1‑150 (MQ=255)
                 aTGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAg                                                                                                                                      >  1:1828032/1‑150 (MQ=255)
                     cGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCa                                                                                                                                  <  2:3199302/150‑1 (MQ=255)
                      gACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCag                                                                                                                                 >  2:1488413/1‑150 (MQ=255)
                      gACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCag                                                                                                                                 >  1:1733332/1‑150 (MQ=255)
                       aCAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCaga                                                                                                                                >  2:3276745/1‑150 (MQ=255)
                        cAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCagag                                                                                                                               >  2:2355961/1‑150 (MQ=255)
                         aaaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAgaga                                                                                                                              >  1:974760/1‑150 (MQ=255)
                           aaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAc                                                                                                                            <  1:2607531/150‑1 (MQ=255)
                            aTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAAcc                                                                                                                           >  1:1178563/1‑150 (MQ=255)
                             ttCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAccc                                                                                                                          <  1:2874274/150‑1 (MQ=255)
                                aaTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt                                                                                                                       >  1:1965746/1‑150 (MQ=255)
                                aaTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt                                                                                                                       <  1:2165936/150‑1 (MQ=255)
                                  tACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtgt                                                                                                                     <  1:234968/150‑1 (MQ=255)
                                  tACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtgt                                                                                                                     >  1:3716734/1‑150 (MQ=255)
                                    cTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGa                                                                                                                   <  2:632521/150‑1 (MQ=255)
                                     tGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGaa                                                                                                                  <  2:444001/150‑1 (MQ=255)
                                      ggAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAAc                                                                                                                 >  1:16337/1‑150 (MQ=255)
                                        aagaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAg                                                                                                               >  1:4009728/1‑150 (MQ=255)
                                        aagaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAg                                                                                                               >  1:14273/1‑150 (MQ=255)
                                           aagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTc                                                                                                           <  2:473058/150‑1 (MQ=255)
                                            agTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTc                                                                                                           <  2:1176424/150‑1 (MQ=255)
                                             gTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCa                                                                                                          <  2:2488325/150‑1 (MQ=255)
                                             gTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCa                                                                                                          <  2:3525205/150‑1 (MQ=255)
                                               tCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACg                                                                                                        >  1:3237476/1‑150 (MQ=255)
                                                   tACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCac                                                                                                    >  2:56184/1‑149 (MQ=255)
                                                     cGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAgg                                                                                                  <  2:3696605/150‑1 (MQ=255)
                                                      gAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAggg                                                                                                 >  2:312161/1‑150 (MQ=255)
                                                       aGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGt                                                                                                <  1:3819501/150‑1 (MQ=255)
                                                        ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTa                                                                                               >  1:2201393/1‑150 (MQ=255)
                                                                   cAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACTACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAcg                                                                                    <  2:1515625/150‑1 (MQ=255)
                                                                   cAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAcg                                                                                    >  1:2613784/1‑150 (MQ=255)
                                                                       aaGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTa                                                                                <  2:2952035/150‑1 (MQ=255)
                                                                       aaGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTa                                                                                >  1:1353225/1‑150 (MQ=255)
                                                                       aaGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTa                                                                                <  2:4041139/150‑1 (MQ=255)
                                                                        aGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCACGGGTAAAAATGCGACGCGTAc                                                                               >  1:256871/1‑150 (MQ=255)
                                                                        aGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTAc                                                                               <  2:1943851/150‑1 (MQ=255)
                                                                        aGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGGAc                                                                               <  2:1148162/150‑1 (MQ=255)
                                                                          gCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCg                                                                             <  2:3743665/150‑1 (MQ=255)
                                                                           ccAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGa                                                                            >  2:61275/1‑150 (MQ=255)
                                                                            cAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGaa                                                                           >  1:3600647/1‑150 (MQ=255)
                                                                             aGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAg                                                                          >  2:1375323/1‑150 (MQ=255)
                                                                             aGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGACTTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAg                                                                          >  2:2891233/1‑150 (MQ=255)
                                                                             aGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAg                                                                          >  1:2280525/1‑150 (MQ=255)
                                                                                 gAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGAt                                                                      >  2:2649348/1‑150 (MQ=255)
                                                                                 gAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGAt                                                                      >  1:2425244/1‑150 (MQ=255)
                                                                                  aGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATc                                                                     <  1:4107732/150‑1 (MQ=255)
                                                                                  aGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATc                                                                     >  1:1671172/1‑150 (MQ=255)
                                                                                   gCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCg                                                                    <  2:70571/150‑1 (MQ=255)
                                                                                   gCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCg                                                                    <  2:2224151/150‑1 (MQ=255)
                                                                                     aaCAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGc                                                                  >  1:2250540/1‑150 (MQ=255)
                                                                                       cAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTa                                                                <  1:749193/150‑1 (MQ=255)
                                                                                              catcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATcc                                                         <  1:637032/150‑1 (MQ=255)
                                                                                                 catTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCt                                                      <  1:2459936/150‑1 (MQ=255)
                                                                                                  atTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCtt                                                     <  1:892837/150‑1 (MQ=255)
                                                                                                      cgccgcCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCg                                                 <  2:397425/150‑1 (MQ=255)
                                                                                                          gccTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAAt                                             >  1:878923/1‑150 (MQ=255)
                                                                                                             tATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAg                                          <  2:1931587/150‑1 (MQ=255)
                                                                                                              aTCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGc                                         >  1:1655689/1‑150 (MQ=255)
                                                                                                                  ccTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCg                                     >  2:4146182/1‑150 (MQ=255)
                                                                                                                   cTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGc                                    >  2:3717241/1‑150 (MQ=255)
                                                                                                                   cTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGc                                    >  2:3717296/1‑150 (MQ=255)
                                                                                                                     gAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCa                                  <  1:1479874/150‑1 (MQ=255)
                                                                                                                       aTAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCAcc                                <  1:3676770/150‑1 (MQ=255)
                                                                                                                       aTAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCAcc                                >  2:2991326/1‑150 (MQ=255)
                                                                                                                        tAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCg                               >  2:3095105/1‑150 (MQ=255)
                                                                                                                          aCCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTc                             >  2:702853/1‑150 (MQ=255)
                                                                                                                          aCCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTc                             <  1:1048013/150‑1 (MQ=255)
                                                                                                                          aCCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTc                             >  1:3037860/1‑150 (MQ=255)
                                                                                                                          aCCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTc                             >  1:537247/1‑150 (MQ=255)
                                                                                                                            ccGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAg                           <  1:3359028/150‑1 (MQ=255)
                                                                                                                             cGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGc                          >  2:449617/1‑150 (MQ=255)
                                                                                                                               gAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGt                        <  1:1000249/150‑1 (MQ=255)
                                                                                                                                aaGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGta                       >  1:3459281/1‑150 (MQ=255)
                                                                                                                                 aGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGtat                      >  1:2024107/1‑150 (MQ=255)
                                                                                                                                   aCAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATaa                    >  1:3098905/1‑150 (MQ=255)
                                                                                                                                      aCCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACtg                 >  2:2427133/1‑150 (MQ=255)
                                                                                                                                         aGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACtgatg              <  2:3632249/150‑1 (MQ=255)
                                                                                                                                          gCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACtgatga             <  1:682407/150‑1 (MQ=255)
                                                                                                                                          gCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACtgatga             >  2:1165665/1‑150 (MQ=255)
                                                                                                                                           cTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTgatgat            <  1:1503007/150‑1 (MQ=255)
                                                                                                                                           cTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTgatgat            >  1:1689444/1‑150 (MQ=255)
                                                                                                                                           cTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTgatgat            >  1:1730231/1‑150 (MQ=255)
                                                                                                                                           cTGTGGTAGCGAATGAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTgatgat            <  1:4147906/150‑1 (MQ=255)
                                                                                                                                             gtgGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATgg          >  1:3651115/1‑150 (MQ=255)
                                                                                                                                              tgGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGc         <  1:1301049/150‑1 (MQ=255)
                                                                                                                                              tgGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGc         <  1:2911053/150‑1 (MQ=255)
                                                                                                                                               ggTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGTAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCt        <  2:2330655/150‑1 (MQ=255)
                                                                                                                                               ggTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCt        <  2:2022144/150‑1 (MQ=255)
                                                                                                                                                  aGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGAt     <  2:3346023/150‑1 (MQ=255)
                                                                                                                                                  aGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGAt     >  1:2713435/1‑150 (MQ=255)
                                                                                                                                                   gCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGAtt    <  2:3138205/150‑1 (MQ=255)
                                                                                                                                                     gAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGc  >  2:3299422/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCCACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGC  >  CP009273/3933531‑3933829

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: