Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,984,987 C→A 34.8% T38T (ACG→ACT)  yecA ← UPF0149 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,984,9870CA34.8% 480.0 / 293.3 296T38T (ACG→ACT) yecAUPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base C (78/115);  new base A (50/53);  total (128/168)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 303 total aligned reads displayed.

TCTGCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCT  >  CP009273/1984838‑1985136
                                                                                                                                                     |                                                                                                                                                     
      gTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAAt                                                                                                                                                 <  2:863502/150‑1 (MQ=255)
       tATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATa                                                                                                                                                >  1:3824367/1‑150 (MQ=255)
       tATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATa                                                                                                                                                >  1:2796309/1‑150 (MQ=255)
        aTCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATaa                                                                                                                                               <  2:1741458/150‑1 (MQ=255)
         tCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATaaa                                                                                                                                              <  1:3109537/150‑1 (MQ=255)
          cGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAAc                                                                                                                                             >  1:3888385/1‑150 (MQ=255)
             ccATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGt                                                                                                                                          >  1:4421226/1‑150 (MQ=255)
             ccATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGt                                                                                                                                          >  2:2527025/1‑150 (MQ=255)
             ccATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGt                                                                                                                                          >  1:4525666/1‑150 (MQ=255)
               atatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTcc                                                                                                                                        >  1:1963475/1‑150 (MQ=255)
               atatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTcc                                                                                                                                        <  2:2024239/150‑1 (MQ=255)
                tatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCa                                                                                                                                       >  2:3716114/1‑150 (MQ=255)
                tatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCa                                                                                                                                       >  2:1145596/1‑150 (MQ=255)
                tatGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCa                                                                                                                                       >  1:4842497/1‑150 (MQ=255)
                   gTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCt                                                                                                                                    <  1:5101446/150‑1 (MQ=255)
                   gTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCt                                                                                                                                    <  2:3652204/150‑1 (MQ=255)
                       aaaaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGc                                                                                                                                >  2:2611316/1‑150 (MQ=255)
                           gCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCcaca                                                                                                                            >  2:2398170/1‑150 (MQ=255)
                              aGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCa                                                                                                                         >  2:4252954/1‑150 (MQ=255)
                               gATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATACACCGTCCAGCTCCGCCACATCaa                                                                                                                        <  2:4978235/150‑1 (MQ=255)
                                    aTAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATg                                                                                                                   <  1:645515/150‑1 (MQ=255)
                                      aaaGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                                 >  1:1456772/1‑150 (MQ=255)
                                      aaaGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGc                                                                                                                 >  1:2284816/1‑150 (MQ=255)
                                        aGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGt                                                                                                               >  1:2915111/1‑150 (MQ=255)
                                         gcgcGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTg                                                                                                              <  2:3162664/150‑1 (MQ=255)
                                         gcgcGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTg                                                                                                              <  2:70992/150‑1 (MQ=255)
                                          cgcgTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTgg                                                                                                             >  1:2029257/1‑150 (MQ=255)
                                           gcgTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGt                                                                                                            <  2:1444254/150‑1 (MQ=255)
                                           gcgTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGt                                                                                                            <  2:2337722/150‑1 (MQ=255)
                                             gTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCa                                                                                                          >  1:3016299/1‑150 (MQ=255)
                                             gTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCa                                                                                                          >  2:3326260/1‑150 (MQ=255)
                                              tCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAg                                                                                                         >  1:261541/1‑150 (MQ=255)
                                              tCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAg                                                                                                         <  1:1663660/150‑1 (MQ=255)
                                               cATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgt                                                                                                        <  2:2736857/150‑1 (MQ=255)
                                                aTCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAgtg                                                                                                       <  2:4527454/150‑1 (MQ=255)
                                                  ctctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgtt                                                                                                     <  1:4084073/150‑1 (MQ=255)
                                                  ctctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgtt                                                                                                     <  1:3023533/150‑1 (MQ=255)
                                                  ctctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgtt                                                                                                     <  1:2895861/150‑1 (MQ=255)
                                                   tctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgttg                                                                                                    <  2:2830685/150‑1 (MQ=255)
                                                    ctTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGtgttgt                                                                                                   >  2:4480997/1‑150 (MQ=255)
                                                     tttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTa                                                                                                  <  2:1614860/150‑1 (MQ=255)
                                                      ttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAt                                                                                                 >  1:3799214/1‑150 (MQ=255)
                                                      ttCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAt                                                                                                 <  2:5060986/150‑1 (MQ=255)
                                                        ctctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAttt                                                                                               <  1:3589681/150‑1 (MQ=255)
                                                        ctctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAttt                                                                                               <  1:2791602/150‑1 (MQ=255)
                                                        ctctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTAttt                                                                                               <  1:1422154/150‑1 (MQ=255)
                                                         tctGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTg                                                                                              <  1:4860292/150‑1 (MQ=255)
                                                            gACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTc                                                                                           >  2:1639534/1‑150 (MQ=255)
                                                                      gcgGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGt                                                                                 >  2:4845260/1‑150 (MQ=255)
                                                                            caTAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGcc                                                                           >  1:156525/1‑150 (MQ=255)
                                                                             aTAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCa                                                                          <  1:3104253/150‑1 (MQ=255)
                                                                              tAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCAt                                                                         >  1:4578366/1‑150 (MQ=255)
                                                                                gTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTc                                                                       <  1:701877/150‑1 (MQ=255)
                                                                                gTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTc                                                                       >  1:2934290/1‑150 (MQ=255)
                                                                                   aGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCaa                                                                    >  1:1177561/1‑150 (MQ=255)
                                                                                          ccccacACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACtt                                                             >  2:3551381/1‑150 (MQ=255)
                                                                                           cccacACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACttt                                                            >  2:4670562/1‑150 (MQ=255)
                                                                                              acacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGt                                                         <  2:855089/150‑1 (MQ=255)
                                                                                               cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGtt                                                        <  2:1547291/150‑1 (MQ=255)
                                                                                                 cGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTa                                                      >  2:1436146/1‑150 (MQ=255)
                                                                                                  ggCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTaa                                                     >  2:5024445/1‑150 (MQ=255)
                                                                                                    cAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACg                                                   <  1:5080529/150‑1 (MQ=255)
                                                                                                       ccAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTc                                                <  2:3494659/150‑1 (MQ=255)
                                                                                                        cAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTcc                                               <  1:2461087/150‑1 (MQ=255)
                                                                                                         aGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTccc                                              <  1:2731417/150‑1 (MQ=255)
                                                                                                                 ttCCGGTTCTATCTCTTGCGGTGAACTCAACCCCGCCGTCAATAATCCGTCCAGCTGCGCGACATAATGGATGGCGTGGTCAGTGTTGTATTTGGTCAGTATATCGTCCAGCCTTTGCAACTCACTTTCGTTTAACGTTCCCGTTTTCAt                                      <  1:3110781/150‑1 (MQ=255)
                                                                                                                     ggTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCACCTCACTGTCGTTTAACGGTCCCGTTTTCATACGc                                  >  1:595767/1‑150 (MQ=255)
                                                                                                                      gTTCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCt                                 <  1:1804153/150‑1 (MQ=255)
                                                                                                                       ttCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCtt                                <  1:3667765/150‑1 (MQ=255)
                                                                                                                       ttCAATCTCTTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCtt                                <  2:4774932/150‑1 (MQ=255)
                                                                                                                        tCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCttt                               >  2:1979245/1‑150 (MQ=255)
                                                                                                                         cAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCtttt                              <  1:1025263/150‑1 (MQ=255)
                                                                                                                         cAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCtttt                              <  1:789705/150‑1 (MQ=255)
                                                                                                                          aaTCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTc                             >  1:476457/1‑150 (MQ=255)
                                                                                                                             ctctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCtt                          >  1:997386/1‑150 (MQ=255)
                                                                                                                              tctTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTtg                         <  2:3382208/150‑1 (MQ=255)
                                                                                                                               ctTGCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTtgt                        <  2:2914997/150‑1 (MQ=255)
                                                                                                                                  gCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGa                     >  2:4786612/1‑150 (MQ=255)
                                                                                                                                  gCGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGa                     <  2:578796/150‑1 (MQ=255)
                                                                                                                                   cGGAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGAATAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAt                    <  2:4419459/150‑1 (MQ=255)
                                                                                                                                    ggAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtc                   >  2:4104380/1‑150 (MQ=255)
                                                                                                                                    ggAGAACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtc                   >  1:1107781/1‑150 (MQ=255)
                                                                                                                                     gagaACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtct                  >  2:493883/1‑150 (MQ=255)
                                                                                                                                      agaACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtctc                 >  1:4276978/1‑150 (MQ=255)
                                                                                                                                      agaACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAtctc                 <  2:759778/150‑1 (MQ=255)
                                                                                                                                       gaACTCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCa                <  1:2072834/150‑1 (MQ=255)
                                                                                                                                        aaCTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCaa               <  2:3964189/150‑1 (MQ=255)
                                                                                                                                           tCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTc            <  2:3193521/150‑1 (MQ=255)
                                                                                                                                           tCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTc            <  2:212571/150‑1 (MQ=255)
                                                                                                                                           tCAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTc            >  2:2704735/1‑150 (MQ=255)
                                                                                                                                            cAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCg           <  1:1399783/150‑1 (MQ=255)
                                                                                                                                            cAACACCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCg           <  2:121188/150‑1 (MQ=255)
                                                                                                                                              acacCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGcc         >  2:4230834/1‑150 (MQ=255)
                                                                                                                                              acacCGCAGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGcc         >  2:4879678/1‑150 (MQ=255)
                                                                                                                                                acCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAg       >  2:2232693/1‑150 (MQ=255)
                                                                                                                                                   gccgTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCAc    <  1:3736348/150‑1 (MQ=255)
                                                                                                                                                   gccgTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCAc    <  1:3387489/150‑1 (MQ=255)
                                                                                                                                                     cgTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCt  <  1:4999300/150‑1 (MQ=255)
                                                                                                                                                     cgTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCt  >  1:4178297/1‑150 (MQ=255)
                                                                                                                                                     cgTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCt  <  1:2580109/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTGCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCT  >  CP009273/1984838‑1985136

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: ‑