Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,404,059 C→A 100% 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,0590CA98.1% 684.1 / 3.2 213R141S (CGC→AGC) dusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base A (100/109);  total (103/110)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 213 total aligned reads displayed.

GATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAG  >  CP009273/3403914‑3404203
                                                                                                                                                 |                                                                                                                                                
 attattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCa                                                                                                                                               <  1:51934/148‑1 (MQ=255)
 attattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCa                                                                                                                                               <  2:2377518/148‑1 (MQ=255)
  ttattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCAc                                                                                                                                              >  1:338687/1‑148 (MQ=255)
  ttattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCAc                                                                                                                                              >  1:338283/1‑148 (MQ=255)
   tattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCAcc                                                                                                                                             >  1:1722712/1‑148 (MQ=255)
     ttGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCgg                                                                                                                                           >  1:1017751/1‑148 (MQ=255)
     ttGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGAAGTACCCGGATGTCGTTAAATCTATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCgg                                                                                                                                           >  1:2001525/1‑148 (MQ=255)
      tGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGc                                                                                                                                          >  1:2855149/1‑148 (MQ=255)
      tGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGc                                                                                                                                          >  2:318631/1‑148 (MQ=255)
       gATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCt                                                                                                                                         <  1:894981/148‑1 (MQ=255)
       gATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCt                                                                                                                                         <  1:348021/148‑1 (MQ=255)
          atCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTggg                                                                                                                                      >  2:2804899/1‑148 (MQ=255)
            cAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCa                                                                                                                                    <  2:21660/148‑1 (MQ=255)
            cAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCa                                                                                                                                    >  1:3154649/1‑148 (MQ=255)
             aaTATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCAc                                                                                                                                   <  1:2112873/148‑1 (MQ=255)
             aaTATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCAc                                                                                                                                   <  1:2157250/148‑1 (MQ=255)
               tatGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCg                                                                                                                                 <  1:247072/148‑1 (MQ=255)
                  gggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGaa                                                                                                                              >  1:2340655/1‑148 (MQ=255)
                  gggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGaa                                                                                                                              <  1:1774485/148‑1 (MQ=255)
                    gTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAaca                                                                                                                            <  2:1488254/148‑1 (MQ=255)
                    gTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAaca                                                                                                                            <  1:300377/148‑1 (MQ=255)
                         ccGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTa                                                                                                                       <  2:2929618/148‑1 (MQ=255)
                          cGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTaa                                                                                                                      <  1:267193/148‑1 (MQ=255)
                          cGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTaa                                                                                                                      <  2:2362075/148‑1 (MQ=255)
                              tAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGc                                                                                                                  >  2:2634783/1‑148 (MQ=255)
                                 aaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGaa                                                                                                               <  1:3049217/148‑1 (MQ=255)
                                 aaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAag                                                                                                              <  2:195553/148‑1 (MQ=255)
                                    aGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAagag                                                                                                            <  2:1176911/148‑1 (MQ=255)
                                      tGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGAt                                                                                                          >  1:2958531/1‑148 (MQ=255)
                                      tGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGAt                                                                                                          >  2:662327/1‑148 (MQ=255)
                                             cAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCaa                                                                                                   >  1:1967475/1‑148 (MQ=255)
                                              aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAAc                                                                                                  <  2:1128839/148‑1 (MQ=255)
                                               agttCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACt                                                                                                 >  1:1888859/4‑148 (MQ=255)
                                                gCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTg                                                                                                >  1:2481107/1‑148 (MQ=255)
                                                gCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTg                                                                                                >  1:2082523/1‑148 (MQ=255)
                                                     cAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGa                                                                                           <  1:2282660/148‑1 (MQ=255)
                                                        gctCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGa                                                                                        <  1:1531649/146‑1 (MQ=255)
                                                        gTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGa                                                                                        >  2:1372676/1‑148 (MQ=255)
                                                        gTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGa                                                                                        >  2:3380034/1‑148 (MQ=255)
                                                            aGCCCTATTGCAGTACACTGATGTCGTTAAATGGATCCTTTCCGTGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAAAACCGTAACTGCGAAGAGTTTGACCAACTTGATGAATACtgt                                                                                    <  2:1717971/148‑1 (MQ=255)
                                                                ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGATCACCGTAACTGCGCAGAGATTGCCCAACTGGCTGAAGACTGTGGCa                                                                                <  2:765552/148‑1 (MQ=255)
                                                                ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCa                                                                                >  2:414833/1‑148 (MQ=255)
                                                                ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCa                                                                                >  1:262932/1‑148 (MQ=255)
                                                                 tctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAt                                                                               <  1:51240/148‑1 (MQ=255)
                                                                   ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTc                                                                             <  2:683115/148‑1 (MQ=255)
                                                                      cAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAgg                                                                          <  2:2156151/148‑1 (MQ=255)
                                                                        gTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGct                                                                        <  2:2146018/148‑1 (MQ=255)
                                                                           cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTg                                                                     <  1:441477/148‑1 (MQ=255)
                                                                              ggATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGAcc                                                                  <  1:1692571/148‑1 (MQ=255)
                                                                               gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCa                                                                 <  2:3122571/148‑1 (MQ=255)
                                                                                aTGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCAt                                                                >  2:1626192/1‑148 (MQ=255)
                                                                                 tGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCAtt                                                               >  1:2616340/1‑148 (MQ=255)
                                                                                 tGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCAtt                                                               >  2:2304689/1‑148 (MQ=255)
                                                                                   tCGTTACATCGATCCGGTCCGAGGTCGTCAATGCAGTGGACGTTCCGGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGACGAACGTACCTGCGAATAGATTGCACAACTGGCGGAAGAATGTGGCATTCAGGCTCTGACAATTCa                                                             <  2:1820333/148‑1 (MQ=255)
                                                                                   tCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCa                                                             >  2:2201859/1‑148 (MQ=255)
                                                                                   tCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCa                                                             <  1:2377558/148‑1 (MQ=255)
                                                                                    cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCAt                                                            <  1:2071235/148‑1 (MQ=255)
                                                                                      ttAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATgg                                                          <  1:2047706/148‑1 (MQ=255)
                                                                                       tAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGc                                                         >  1:1456927/1‑148 (MQ=255)
                                                                                          aTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGt                                                      >  2:2923131/1‑148 (MQ=255)
                                                                                            cGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTac                                                    >  2:2594541/1‑148 (MQ=255)
                                                                                              aTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTacac                                                  >  1:1211243/1‑148 (MQ=255)
                                                                                              aTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTacac                                                  >  1:935779/1‑148 (MQ=255)
                                                                                                ccTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACAcgc                                                >  2:2221298/1‑148 (MQ=255)
                                                                                                  ttACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTACCTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACgcgc                                              <  1:2569826/148‑1 (MQ=255)
                                                                                                  ttACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACgcgc                                              >  2:3278715/1‑148 (MQ=255)
                                                                                                  ttACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACgcgc                                              <  1:3289298/148‑1 (MQ=255)
                                                                                                       gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtt                                         <  1:341272/148‑1 (MQ=255)
                                                                                                         ggTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttg                                       <  1:57970/148‑1 (MQ=255)
                                                                                                           tcgtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttgtt                                     <  1:77118/148‑1 (MQ=255)
                                                                                                           tcgtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttgtt                                     >  1:297167/1‑148 (MQ=255)
                                                                                                               cAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAAt                                 >  1:464357/1‑148 (MQ=255)
                                                                                                                  tGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGga                              <  2:144017/148‑1 (MQ=255)
                                                                                                                    cAGTGGACGTTCCTGTTACCCTGAAGATTAGTACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGgaga                            <  1:43411/148‑1 (MQ=255)
                                                                                                                    cAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGgaga                            <  2:1640702/148‑1 (MQ=255)
                                                                                                                     aGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGaa                           <  1:1372676/148‑1 (MQ=255)
                                                                                                                      gTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAg                          >  2:2711083/1‑148 (MQ=255)
                                                                                                                      gTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAg                          <  2:2933168/148‑1 (MQ=255)
                                                                                                                       tGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGc                         <  1:630281/148‑1 (MQ=255)
                                                                                                                       tGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGc                         >  1:1711862/1‑148 (MQ=255)
                                                                                                                          aCGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGa                      >  1:1922122/1‑148 (MQ=255)
                                                                                                                            gTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGt                    >  1:1748044/1‑148 (MQ=255)
                                                                                                                             ttCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTa                   >  1:2358342/1‑148 (MQ=255)
                                                                                                                             ttCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTa                   >  1:595587/1‑148 (MQ=255)
                                                                                                                             ttCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTa                   <  1:1773845/148‑1 (MQ=255)
                                                                                                                              tCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTAc                  >  2:1090980/1‑148 (MQ=255)
                                                                                                                              tCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTAc                  >  1:1007405/1‑148 (MQ=255)
                                                                                                                                 tGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACGGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGAc               <  2:3310849/148‑1 (MQ=255)
                                                                                                                                  gTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACa              <  2:1967475/148‑1 (MQ=255)
                                                                                                                                    tACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGt            >  1:1386745/1‑148 (MQ=255)
                                                                                                                                    tACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGt            <  1:1255336/148‑1 (MQ=255)
                                                                                                                                    tACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGt            <  2:486545/148‑1 (MQ=255)
                                                                                                                                     aCCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTa           <  2:727973/148‑1 (MQ=255)
                                                                                                                                      cccTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAt          <  2:2409182/148‑1 (MQ=255)
                                                                                                                                        cTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTc        <  1:2067295/148‑1 (MQ=255)
                                                                                                                                           aaGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCggg     >  2:1273097/1‑148 (MQ=255)
                                                                                                                                            aGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGc    <  2:1456927/148‑1 (MQ=255)
                                                                                                                                              aTTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAg  <  2:297167/148‑1 (MQ=255)
                                                                                                                                              aTTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAg  <  1:2634783/148‑1 (MQ=255)
                                                                                                                                              aTTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAg  >  1:2797609/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAG  >  CP009273/3403914‑3404203

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: