Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,404,059 C→A 86.4% R141S (CGC→AGC)  dusB → tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,404,0590CA86.4% 713.6 / 120.0 295R141S (CGC→AGC) dusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base A (122/133);  total (144/151)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 299 total aligned reads displayed.

CCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAA  >  CP009273/3403911‑3404207
                                                                                                                                                    |                                                                                                                                                    
ccAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTc                                                                                                                                                      <  1:277810/149‑1 (MQ=255)
 cAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCg                                                                                                                                                     >  2:1876233/1‑149 (MQ=255)
   gATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCa                                                                                                                                                   >  2:1194891/1‑149 (MQ=255)
   gATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCa                                                                                                                                                   <  2:780719/149‑1 (MQ=255)
   gATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCAGAAGATTAGCa                                                                                                                                                   >  2:4389452/1‑149 (MQ=255)
    attattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCAc                                                                                                                                                  >  2:4005388/1‑149 (MQ=255)
     ttattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAcc                                                                                                                                                 >  2:3371025/1‑149 (MQ=255)
        ttGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGc                                                                                                                                              <  1:3400425/149‑1 (MQ=255)
        ttGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGc                                                                                                                                              <  2:4374663/149‑1 (MQ=255)
         tGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCt                                                                                                                                             <  2:2407909/149‑1 (MQ=255)
         tGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCt                                                                                                                                             >  2:1809863/1‑149 (MQ=255)
          gATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTg                                                                                                                                            >  2:3322053/1‑149 (MQ=255)
          gATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTg                                                                                                                                            >  1:96259/1‑149 (MQ=255)
           atatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTgg                                                                                                                                           >  2:2505038/1‑149 (MQ=255)
             atCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGc                                                                                                                                         >  1:345566/1‑149 (MQ=255)
             atCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGc                                                                                                                                         >  1:673935/1‑149 (MQ=255)
              tCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCa                                                                                                                                        >  2:1260184/1‑149 (MQ=255)
                  tatGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCgg                                                                                                                                    >  2:1011832/1‑149 (MQ=255)
                  tatGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCgg                                                                                                                                    >  1:2562/1‑149 (MQ=255)
                        ttGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACAcc                                                                                                                              <  1:2181444/149‑1 (MQ=255)
                        ttGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACAcc                                                                                                                              <  1:1846722/149‑1 (MQ=255)
                          gCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGt                                                                                                                            >  1:4087899/1‑149 (MQ=255)
                              ggCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACt                                                                                                                        <  1:2747651/149‑1 (MQ=255)
                               gCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTg                                                                                                                       <  1:2919627/149‑1 (MQ=255)
                                 tAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCg                                                                                                                     <  2:4715640/149‑1 (MQ=255)
                                  aaaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGa                                                                                                                    >  1:1340848/1‑149 (MQ=255)
                                   aaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGaa                                                                                                                   >  2:625236/1‑149 (MQ=255)
                                       aGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAgaga                                                                                                               <  1:996568/149‑1 (MQ=255)
                                             tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGccc                                                                                                         <  1:4323047/149‑1 (MQ=255)
                                             tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGccc                                                                                                         <  1:3037647/149‑1 (MQ=255)
                                             tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGccc                                                                                                         <  2:727674/149‑1 (MQ=255)
                                               gCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCaa                                                                                                       <  1:4585438/149‑1 (MQ=255)
                                               gCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCaa                                                                                                       >  2:1966445/1‑149 (MQ=255)
                                                cAAGCGCGCAGGTTCAGCCCTCGAGCAGTACCCGGATGTCGTGACATCGATCCGTACCGAGGTCGTCCATGCAGTGGACGTTCCTGTTACCCGGAAGATTAGCACCGGCTGGGCACCGGAACACCGTACCTGCGACGAGATTGCCCacc                                                                                                      <  2:2303411/149‑3 (MQ=255)
                                                 aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACt                                                                                                     >  2:4488318/1‑149 (MQ=255)
                                                 aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACt                                                                                                     >  2:3163148/1‑149 (MQ=255)
                                                  aGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTg                                                                                                    <  1:4246930/149‑1 (MQ=255)
                                                  aGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTg                                                                                                    >  1:26911/1‑149 (MQ=255)
                                                            ttCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACt                                                                                          <  2:3636679/149‑1 (MQ=255)
                                                             tCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACtg                                                                                         <  2:1454855/149‑1 (MQ=255)
                                                              cAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACtgt                                                                                        <  2:533778/149‑1 (MQ=255)
                                                               aGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACtgtg                                                                                       >  2:890310/1‑149 (MQ=255)
                                                                    tctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                                  <  1:3990821/149‑1 (MQ=255)
                                                                     cgtGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATGAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTc                                                                                 <  1:3538967/147‑1 (MQ=255)
                                                                      ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCa                                                                                <  2:3889509/149‑1 (MQ=255)
                                                                           gTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGctc                                                                           <  2:1727250/149‑1 (MQ=255)
                                                                            tACCCGGATGTCGTTACATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGctct                                                                          <  2:4126971/149‑1 (MQ=255)
                                                                                   aTGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCAtt                                                                   <  2:2942257/149‑1 (MQ=255)
                                                                                     gTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTTGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCa                                                                 <  2:3950554/149‑1 (MQ=255)
                                                                                     gTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCa                                                                 <  1:1710756/149‑1 (MQ=255)
                                                                                     gTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCa                                                                 <  1:4117054/149‑1 (MQ=255)
                                                                                     gTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCa                                                                 <  1:1074632/149‑1 (MQ=255)
                                                                                     gTCGTTAAATCGATCCTTACCGAGGTCGTCAAGGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCa                                                                 <  2:150356/149‑1 (MQ=255)
                                                                                      tCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCAt                                                                <  2:2790370/149‑1 (MQ=255)
                                                                                      tCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCAt                                                                >  1:3860543/1‑149 (MQ=255)
                                                                                       cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATg                                                               <  2:2550572/149‑1 (MQ=255)
                                                                                       cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATg                                                               <  1:1428512/149‑1 (MQ=255)
                                                                                         ttAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGc                                                             <  1:1867710/149‑1 (MQ=255)
                                                                                         ttAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGc                                                             <  1:4257557/149‑1 (MQ=255)
                                                                                         ttAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGc                                                             >  1:2852261/1‑149 (MQ=255)
                                                                                          aaaaTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCAGTTACCCAGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGcc                                                            <  2:3471093/148‑1 (MQ=255)
                                                                                            aaTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGt                                                          <  1:3221880/149‑1 (MQ=255)
                                                                                             aTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTa                                                         >  2:4149392/1‑149 (MQ=255)
                                                                                              tCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTac                                                        <  1:3593218/149‑1 (MQ=255)
                                                                                               cGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTaca                                                       >  1:3956853/1‑149 (MQ=255)
                                                                                               cGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTaca                                                       >  1:1038801/1‑149 (MQ=255)
                                                                                                   ccTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACAcgcg                                                   >  2:42363/1‑149 (MQ=255)
                                                                                                    cTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACgcgc                                                  >  1:1002085/1‑149 (MQ=255)
                                                                                                      tACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCt                                                >  2:2081322/1‑149 (MQ=255)
                                                                                                      tACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCt                                                <  2:39673/149‑1 (MQ=255)
                                                                                                         cGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtt                                             <  1:4442242/149‑1 (MQ=255)
                                                                                                          gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGTCGTACACGCGCCtgttt                                            >  2:4605913/1‑149 (MQ=255)
                                                                                                          gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgttt                                            >  1:1206114/1‑149 (MQ=255)
                                                                                                          gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgttt                                            >  1:2280482/1‑149 (MQ=255)
                                                                                                             gtcgtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttgtt                                         >  1:3222683/1‑149 (MQ=255)
                                                                                                                gtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCaa                                      <  1:2817712/149‑1 (MQ=255)
                                                                                                                 tcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAAt                                     >  1:2402914/1‑149 (MQ=255)
                                                                                                                 tcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAAt                                     >  1:1187230/1‑149 (MQ=255)
                                                                                                                      gCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGgaga                                >  2:2032208/1‑149 (MQ=255)
                                                                                                                      gCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGgaga                                <  2:4612200/149‑1 (MQ=255)
                                                                                                                      gCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGgaga                                <  2:3500627/149‑1 (MQ=255)
                                                                                                                      gCAGTGGACGGTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGgaga                                <  1:3610327/149‑1 (MQ=255)
                                                                                                                       cAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGaa                               <  2:2319468/149‑1 (MQ=255)
                                                                                                                         gTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGc                             >  2:2660115/1‑149 (MQ=255)
                                                                                                                           ggACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTg                           >  1:806026/1‑149 (MQ=255)
                                                                                                                            gACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGa                          <  1:4617136/149‑1 (MQ=255)
                                                                                                                            gACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGa                          <  1:931308/149‑1 (MQ=255)
                                                                                                                            gACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGa                          >  2:486311/1‑149 (MQ=255)
                                                                                                                             aCGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAg                         >  1:3814413/1‑149 (MQ=255)
                                                                                                                              cGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGt                        >  1:984977/1‑149 (MQ=255)
                                                                                                                              cGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGt                        >  1:596676/1‑149 (MQ=255)
                                                                                                                               ggtCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTa                       <  1:3697692/147‑1 (MQ=255)
                                                                                                                               gTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTa                       >  2:1090144/1‑149 (MQ=255)
                                                                                                                                ttCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTAc                      >  1:4668373/1‑149 (MQ=255)
                                                                                                                                ttCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTAc                      >  2:3061430/1‑149 (MQ=255)
                                                                                                                                        aCCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAt              <  2:3784501/149‑1 (MQ=255)
                                                                                                                                         cccTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAtt             <  1:627400/149‑1 (MQ=255)
                                                                                                                                         cccTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAtt             <  1:877091/149‑1 (MQ=255)
                                                                                                                                         cccTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAtt             <  1:476396/149‑1 (MQ=255)
                                                                                                                                            tGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCgg          >  1:610642/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAA  >  CP009273/3403911‑3404207

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: