Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,404,059 C→A 6.6% 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,0590CA 6.6% 1035.5 / 32.2 270R141S (CGC→AGC) dusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base C (117/135);  new base A (12/6);  total (129/141)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 272 total aligned reads displayed.

CCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAA  >  CP009273/3403910‑3404207
                                                                                                                                                     |                                                                                                                                                    
cccAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTc                                                                                                                                                      >  2:1582523/1‑150 (MQ=255)
  cAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGc                                                                                                                                                    >  1:399576/1‑150 (MQ=255)
   aGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCa                                                                                                                                                   >  2:2153520/1‑150 (MQ=255)
    gATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAc                                                                                                                                                  >  1:3470377/1‑150 (MQ=255)
     attattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAcc                                                                                                                                                 >  2:3014475/1‑150 (MQ=255)
     attattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAcc                                                                                                                                                 >  1:1741820/1‑150 (MQ=255)
      ttattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCg                                                                                                                                                <  2:3830398/150‑1 (MQ=255)
           gATATCAATATGGGTTGCCCGGCTCAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTgg                                                                                                                                           >  2:1247638/1‑150 (MQ=255)
           gATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTgg                                                                                                                                           >  1:2812177/1‑150 (MQ=255)
            atatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg                                                                                                                                          <  2:3606780/150‑1 (MQ=255)
            atatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGAAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTggg                                                                                                                                          >  1:2448481/1‑150 (MQ=255)
              atCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCa                                                                                                                                        >  2:1275505/1‑150 (MQ=255)
               tCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACACTGAAGATTCGCACCGGCTGGGCAc                                                                                                                                       >  1:2745597/1‑150 (MQ=255)
                cAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCAcc                                                                                                                                      <  1:2769893/150‑1 (MQ=255)
                 aaTATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg                                                                                                                                     >  1:2515198/1‑150 (MQ=255)
                 aaTATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg                                                                                                                                     >  2:4159063/1‑150 (MQ=255)
                  atatGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg                                                                                                                                    >  2:2332242/1‑150 (MQ=255)
                   tatGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGa                                                                                                                                   >  2:582007/1‑150 (MQ=255)
                    atGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGaa                                                                                                                                  >  2:2748393/1‑150 (MQ=255)
                     tGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAac                                                                                                                                 <  2:4134116/150‑1 (MQ=255)
                     tGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAac                                                                                                                                 >  2:3755897/1‑150 (MQ=255)
                      gggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAaca                                                                                                                                >  1:134643/1‑150 (MQ=255)
                       ggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAacac                                                                                                                               <  2:2285941/150‑1 (MQ=255)
                        ggtGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACAcc                                                                                                                              <  2:1783543/148‑1 (MQ=255)
                             ccGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAAc                                                                                                                         >  1:679749/1‑150 (MQ=255)
                              cGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACt                                                                                                                        <  1:1915148/150‑1 (MQ=255)
                                gCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGc                                                                                                                      >  1:3798851/1‑150 (MQ=255)
                                 cTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCg                                                                                                                     <  2:711026/150‑1 (MQ=255)
                                  tAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGACCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGa                                                                                                                    <  1:3879905/150‑1 (MQ=255)
                                   aaaaaaGTGAATCGCAAGCTCGAAGGTTCAGCCCTCTTGCAGTACACTGATGTCTTTAAATCGATACATACCGAGGTCTGCAATGCATTGGACGTTACTTTTAACCTGAATATTCGCACCTGCTTTGCACCGGAACACCGTAACTGCGaa                                                                                                                   >  2:1951087/1‑150 (MQ=255)
                                    aaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAag                                                                                                                  >  1:2697036/1‑150 (MQ=255)
                                      aaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAagag                                                                                                                >  1:3067785/1‑150 (MQ=255)
                                       aaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAgaga                                                                                                               >  1:1437604/1‑150 (MQ=255)
                                         gTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGAtt                                                                                                             <  1:3341554/150‑1 (MQ=255)
                                            aaTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGcc                                                                                                          >  2:3271041/1‑150 (MQ=255)
                                             atagCAATCTCTCAGGTTCAGCCCTATTGCAGTACCCGGATGTCGTTAAATCTATCATTACAGATGTCGTCAATGCAGTGGACGTTCCAGATACCCTGAAGATTCGCACCGGCTGGGCACCGGAAAACCGTAACTGAGAAGAGATTGacc                                                                                                         >  2:2654377/4‑150 (MQ=255)
                                              tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCa                                                                                                        <  2:1969886/150‑1 (MQ=255)
                                               cGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCaa                                                                                                       >  2:4630855/1‑150 (MQ=255)
                                                gCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAAc                                                                                                      <  1:1645720/150‑1 (MQ=255)
                                                 cAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTTGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACt                                                                                                     >  2:1431523/1‑150 (MQ=255)
                                                 cAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACt                                                                                                     >  1:4586090/1‑150 (MQ=255)
                                                  aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTg                                                                                                    >  2:3683295/1‑150 (MQ=255)
                                                      tCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTg                                                                                                >  1:2567631/1‑150 (MQ=255)
                                                                aGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTgg                                                                                      <  1:4266944/150‑1 (MQ=255)
                                                                aGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTgg                                                                                      >  2:3240486/1‑150 (MQ=255)
                                                                   ccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAt                                                                                   <  2:3474868/150‑1 (MQ=255)
                                                                   ccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAt                                                                                   <  2:331352/150‑1 (MQ=255)
                                                                    ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                                  <  2:2732736/150‑1 (MQ=255)
                                                                    ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                                  <  2:4045999/150‑1 (MQ=255)
                                                                         gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGc                                                                             >  1:931835/1‑150 (MQ=255)
                                                                         gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGc                                                                             >  1:2248545/1‑150 (MQ=255)
                                                                           aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGCACACCGTAACTGCGAAGAGATTGCCCAACTCGCTGAAGACTGTGGCATTCAGGctc                                                                           >  2:2478163/1‑150 (MQ=255)
                                                                           aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGctc                                                                           <  1:4545547/150‑1 (MQ=255)
                                                                             tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTg                                                                         <  1:2301113/150‑1 (MQ=255)
                                                                               cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGAc                                                                       <  1:863110/150‑1 (MQ=255)
                                                                                 cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTTCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCa                                                                     <  2:257156/150‑1 (MQ=255)
                                                                                 cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCa                                                                     <  1:3681708/150‑1 (MQ=255)
                                                                                   gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCAtt                                                                   <  1:63253/150‑1 (MQ=255)
                                                                                        cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATgg                                                              >  1:1015158/1‑150 (MQ=255)
                                                                                         gTTAAATCGATCCTTACCGAGGTCGTCAATNCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGc                                                             >  1:4282135/1‑150 (MQ=255)
                                                                                         gTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGc                                                             >  1:736530/1‑150 (MQ=255)
                                                                                         gTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGc                                                             <  2:1918181/150‑1 (MQ=255)
                                                                                           tAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCg                                                           >  2:3567404/1‑150 (MQ=255)
                                                                                            aaaTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGt                                                          >  1:2419188/1‑150 (MQ=255)
                                                                                              aTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTac                                                        <  2:1486936/150‑1 (MQ=255)
                                                                                                 gATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACAcg                                                     <  1:2979655/150‑1 (MQ=255)
                                                                                                    ccTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACgcgc                                                  >  1:2858236/1‑150 (MQ=255)
                                                                                                      ttACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCt                                                <  1:944138/150‑1 (MQ=255)
                                                                                                      ttACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCt                                                >  1:1530649/1‑150 (MQ=255)
                                                                                                           gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttg                                           >  1:3919281/1‑150 (MQ=255)
                                                                                                           gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttg                                           >  1:1318371/1‑150 (MQ=255)
                                                                                                            aGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGACGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttgt                                          >  2:4069066/1‑150 (MQ=255)
                                                                                                                cgtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCaa                                      >  2:129667/1‑150 (MQ=255)
                                                                                                                cgtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCaa                                      <  1:2609266/150‑1 (MQ=255)
                                                                                                                  tcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATg                                    >  1:3330401/1‑150 (MQ=255)
                                                                                                                  tcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATg                                    <  2:4166649/150‑1 (MQ=255)
                                                                                                                  tcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTAAACGCGCCTGTTTGTTCAATg                                    <  2:1363548/150‑1 (MQ=255)
                                                                                                                    aaTGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGga                                  >  2:874987/1‑150 (MQ=255)
                                                                                                                         aGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGc                             >  1:3011478/1‑150 (MQ=255)
                                                                                                                          gTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCt                            <  1:497301/150‑1 (MQ=255)
                                                                                                                            ggACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGa                          <  1:1717965/150‑1 (MQ=255)
                                                                                                                               cGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTa                       <  1:235782/150‑1 (MQ=255)
                                                                                                                                gTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTAc                      >  1:1362221/1‑150 (MQ=255)
                                                                                                                                gTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTAc                      >  2:4284254/1‑150 (MQ=255)
                                                                                                                                 ttCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACg                     >  1:1631721/1‑150 (MQ=255)
                                                                                                                                   ccTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGAc                   <  1:4484170/150‑1 (MQ=255)
                                                                                                                                    cTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACa                  <  1:1670939/150‑1 (MQ=255)
                                                                                                                                     tGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAg                 <  1:2147493/150‑1 (MQ=255)
                                                                                                                                      gTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGt                <  2:3123307/150‑1 (MQ=255)
                                                                                                                                        tACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAt              >  1:1113539/1‑150 (MQ=255)
                                                                                                                                         aCCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAtt             >  1:3574296/1‑150 (MQ=255)
                                                                                                                                         aCCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAtt             <  1:3986153/150‑1 (MQ=255)
                                                                                                                                         aCCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAtt             <  1:4528467/150‑1 (MQ=255)
                                                                                                                                            cTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCgg          <  2:1886330/150‑1 (MQ=255)
                                                                                                                                               aaGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCa       <  2:2885979/150‑1 (MQ=255)
                                                                                                                                                aGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAg      <  2:1077731/150‑1 (MQ=255)
                                                                                                                                                   ttCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTa   <  1:2214217/150‑1 (MQ=255)
                                                                                                                                                   ttCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTa   <  2:2852988/150‑1 (MQ=255)
                                                                                                                                                    tCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTaa  <  1:2976894/150‑1 (MQ=255)
                                                                                                                                                    tCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTaa  <  2:3437785/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAA  >  CP009273/3403910‑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: