Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,933,680 G→T 25.8% 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,6800GT25.8% 771.4 / 233.4 339S98* (TCG→TAG) hsrAputative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base G (132/119);  new base T (41/47);  total (173/166)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 345 total aligned reads displayed.

CCCACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGC  >  CP009273/3933533‑3933829
                                                                                                                                                   |                                                                                                                                                     
 ccaccaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAAt                                                                                                                                                    >  1:4698389/1‑150 (MQ=255)
  caccaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAAtt                                                                                                                                                   <  2:4423307/150‑1 (MQ=255)
  caccaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAAtt                                                                                                                                                   <  2:591698/150‑1 (MQ=255)
   accaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTa                                                                                                                                                  <  1:1408341/150‑1 (MQ=255)
   accaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTa                                                                                                                                                  <  2:320374/150‑1 (MQ=255)
     caGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGa                                                                                                                                                >  1:2444428/1‑150 (MQ=255)
      aGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGaa                                                                                                                                               >  1:3884688/1‑150 (MQ=255)
      aGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGaa                                                                                                                                               >  1:89196/1‑150 (MQ=255)
        aCCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAg                                                                                                                                             <  2:167152/150‑1 (MQ=255)
          ccGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTg                                                                                                                                           <  1:346542/150‑1 (MQ=255)
             gCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTgcgc                                                                                                                                        >  1:2107790/1‑150 (MQ=255)
             gCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcgc                                                                                                                                        >  1:2256610/1‑150 (MQ=255)
                tGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAgg                                                                                                                                     <  1:2715319/150‑1 (MQ=255)
                   cGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCa                                                                                                                                  <  1:1551662/150‑1 (MQ=255)
                     aCAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCaga                                                                                                                                >  1:3248532/1‑150 (MQ=255)
                        aaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGaa                                                                                                                             >  2:3466453/1‑150 (MQ=255)
                           ttCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAccc                                                                                                                          >  2:2779100/1‑150 (MQ=255)
                            tCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCa                                                                                                                         >  2:1869929/1‑150 (MQ=255)
                               aTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtg                                                                                                                      >  1:1670349/1‑150 (MQ=255)
                               aTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtg                                                                                                                      <  1:3214942/150‑1 (MQ=255)
                                 aCTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtgtg                                                                                                                    >  1:2874094/1‑150 (MQ=255)
                                   tGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGaa                                                                                                                  <  2:2615203/150‑1 (MQ=255)
                                   tGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGaa                                                                                                                  >  2:1319566/1‑150 (MQ=255)
                                   tGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGaa                                                                                                                  <  2:224445/150‑1 (MQ=255)
                                      aagaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAg                                                                                                               >  1:4489563/1‑150 (MQ=255)
                                        gaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGAc                                                                                                             >  1:207787/1‑150 (MQ=255)
                                        gaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGAGCGCAGGCCAGAGAACCCAATGTGAACAGAc                                                                                                             >  1:3798771/1‑150 (MQ=255)
                                         aagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAACAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACt                                                                                                            <  2:2011814/150‑1 (MQ=255)
                                            ttCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCAc                                                                                                         >  2:1703196/1‑150 (MQ=255)
                                               aTTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGc                                                                                                      >  1:3921480/1‑150 (MQ=255)
                                                 tACGAGGATAATCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAAACAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCaa                                                                                                    >  2:2365560/1‑150 (MQ=255)
                                                 tACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATACCCCGGAAGACAACCAGCTGTGGTAGCGACTTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGcca                                                                                                    >  2:860888/1‑148 (MQ=255)
                                                     aGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGt                                                                                                >  2:656290/1‑150 (MQ=255)
                                                       gATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTaa                                                                                              >  1:832885/1‑150 (MQ=255)
                                                        aTAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTaaa                                                                                             >  1:3624427/1‑150 (MQ=255)
                                                          aaGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTaaaaa                                                                                           <  2:322474/150‑1 (MQ=255)
                                                           aGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAAt                                                                                          <  2:1150088/150‑1 (MQ=255)
                                                           aGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAAt                                                                                          <  1:668054/150‑1 (MQ=255)
                                                             cgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAAGTAGAAAGGGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCACGGGTCCAAATGc                                                                                        >  2:772757/1‑150 (MQ=255)
                                                                gcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAc                                                                                     >  2:4512819/1‑150 (MQ=255)
                                                                gcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAc                                                                                     <  1:2143339/150‑1 (MQ=255)
                                                                     aaGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTCATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCACGGGTAAAAATGCGACGCGTa                                                                                >  1:262454/1‑150 (MQ=255)
                                                                      aGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTAc                                                                               <  2:2941078/150‑1 (MQ=255)
                                                                           aGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAg                                                                          <  2:3917189/150‑1 (MQ=255)
                                                                             ccGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCg                                                                        <  1:985046/150‑1 (MQ=255)
                                                                              cGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGa                                                                       >  1:1894657/1‑150 (MQ=255)
                                                                              cGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGa                                                                       >  1:5141/1‑150 (MQ=255)
                                                                              cGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGa                                                                       >  1:3307824/1‑150 (MQ=255)
                                                                                aGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATc                                                                     <  2:2444428/150‑1 (MQ=255)
                                                                                aGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATc                                                                     <  2:3998428/150‑1 (MQ=255)
                                                                                  cAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCgg                                                                   >  2:1500862/1‑150 (MQ=255)
                                                                                   aaCAGGCATCATCATTGCTCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGc                                                                  >  2:2098253/1‑150 (MQ=255)
                                                                                    aCAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGCAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCt                                                                 <  1:403884/150‑1 (MQ=255)
                                                                                    aCAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCt                                                                 <  2:2415282/150‑1 (MQ=255)
                                                                                     cAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTa                                                                <  2:347322/150‑1 (MQ=255)
                                                                                      aGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAg                                                               <  1:3639881/150‑1 (MQ=255)
                                                                                        gCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGTAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGcc                                                             <  1:2608534/150‑1 (MQ=255)
                                                                                          atcatcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCAt                                                           <  1:4714733/150‑1 (MQ=255)
                                                                                          atcatcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCAt                                                           <  1:1211602/150‑1 (MQ=255)
                                                                                               catTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCt                                                      <  1:2038605/150‑1 (MQ=255)
                                                                                                atTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCACTGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCtt                                                     <  2:4231141/150‑1 (MQ=255)
                                                                                                    cgccgcCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCg                                                 <  1:1824566/150‑1 (MQ=255)
                                                                                                    cgccgcCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCg                                                 <  2:875601/150‑1 (MQ=255)
                                                                                                     gccgccTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCACGGGTAAAATTGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCgg                                                >  1:2991896/1‑150 (MQ=255)
                                                                                                        gccTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGTCATCCGCTTACCGGAAt                                             <  1:4592652/150‑1 (MQ=255)
                                                                                                        gccTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAAt                                             <  1:1978731/150‑1 (MQ=255)
                                                                                                               cccTGAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATc                                      <  1:3453956/150‑1 (MQ=255)
                                                                                                               cccTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATc                                      <  2:1416319/150‑1 (MQ=255)
                                                                                                                   gAATAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCa                                  <  1:3567593/150‑1 (MQ=255)
                                                                                                                    aaTAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCAc                                 <  1:4192515/150‑1 (MQ=255)
                                                                                                                    aaTAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCAc                                 >  2:2505600/1‑150 (MQ=255)
                                                                                                                     aTAACCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCAcc                                >  2:4776758/1‑150 (MQ=255)
                                                                                                                     aTAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCAcc                                >  1:1424706/1‑150 (MQ=255)
                                                                                                                       aaCCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGt                              >  2:2710264/1‑150 (MQ=255)
                                                                                                                       aaCCCGGAAGACAACCAGCTGTGGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGt                              >  2:314840/1‑150 (MQ=255)
                                                                                                                         cccGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCa                            >  1:3077333/1‑150 (MQ=255)
                                                                                                                           cGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGc                          >  1:4259859/1‑150 (MQ=255)
                                                                                                                            ggAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCg                         >  2:389165/1‑150 (MQ=255)
                                                                                                                              aaGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGta                       <  1:2104215/150‑1 (MQ=255)
                                                                                                                              aaGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGta                       <  1:1319566/150‑1 (MQ=255)
                                                                                                                                 aCAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATaa                    <  2:1838954/150‑1 (MQ=255)
                                                                                                                                 aCAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATaa                    >  1:1093336/1‑150 (MQ=255)
                                                                                                                                  cAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAAc                   <  1:3520484/150‑1 (MQ=255)
                                                                                                                                   aaCCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACt                  <  2:2569216/150‑1 (MQ=255)
                                                                                                                                     ccAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGACTAAGCATCGCCACCGTCAGCGTATAACtga                >  2:4332643/1‑150 (MQ=255)
                                                                                                                                      cAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCACTGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACtgat               <  2:3000790/150‑1 (MQ=255)
                                                                                                                                       aGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACtgatg              >  1:2554181/1‑150 (MQ=255)
                                                                                                                                           gtgGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATgg          <  1:3774494/150‑1 (MQ=255)
                                                                                                                                           gtgGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATgg          >  2:213765/1‑150 (MQ=255)
                                                                                                                                           gtgGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATgg          <  2:2672918/150‑1 (MQ=255)
                                                                                                                                            tgGTAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGc         <  2:354738/150‑1 (MQ=255)
                                                                                                                                            tgGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGc         >  2:408560/1‑150 (MQ=255)
                                                                                                                                            tgGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGc         >  2:3096275/1‑150 (MQ=255)
                                                                                                                                             ggTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCACTGGGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCt        <  2:4541903/150‑1 (MQ=255)
                                                                                                                                               tAGCTAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGa      <  1:760758/150‑1 (MQ=255)
                                                                                                                                                aGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGTCACCGTCAGCGTATAACTGATGATGGCTGAt     <  2:3033501/150‑1 (MQ=255)
                                                                                                                                                aGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGAt     <  1:1777012/150‑1 (MQ=255)
                                                                                                                                                  cGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTg   <  1:2779112/150‑1 (MQ=255)
                                                                                                                                                   gAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATCACTGATGATGGCTGATTGc  >  1:451822/1‑150 (MQ=255)
                                                                                                                                                   gAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGc  >  1:1471535/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCCACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTGATGATGGCTGATTGC  >  CP009273/3933533‑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: