Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,005,091 T→A 5.4% C598S (TGC→AGC)  rlmL → fused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,005,0910TA 5.4% 637.8 / 20.6 205C598S (TGC→AGC) rlmLfused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (86/107);  new base A (5/6);  total (91/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 208 total aligned reads displayed.

GGCAGCGCCACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGAAACGTCTGCTGC  >  CP009273/1004947‑1005238
                                                                                                                                                |                                                                                                                                                   
ggCAGCGCCACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGca                                                                                                                                                  <  1:527809/148‑2 (MQ=255)
  cAGCGCCACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTg                                                                                                                                                >  1:1967629/1‑148 (MQ=255)
     cgcCACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCg                                                                                                                                             <  2:198385/148‑1 (MQ=255)
     cgcCACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCg                                                                                                                                             >  2:1420425/1‑148 (MQ=255)
       ccACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTg                                                                                                                                           <  1:1295722/148‑1 (MQ=255)
       ccACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATAGCCTGGCGTg                                                                                                                                           >  1:1178637/1‑148 (MQ=255)
        cACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTgg                                                                                                                                          <  1:2510296/148‑1 (MQ=255)
        cACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATAGCCTGGCGTgg                                                                                                                                          <  1:440608/148‑1 (MQ=255)
        cACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATAGCCTGGCGTgg                                                                                                                                          <  1:440630/148‑1 (MQ=255)
         aCCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGc                                                                                                                                         <  2:260712/148‑1 (MQ=255)
          ccGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCt                                                                                                                                        >  2:491020/1‑148 (MQ=255)
            gTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGc                                                                                                                                      <  1:3334733/148‑1 (MQ=255)
              gCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGt                                                                                                                                    >  1:345644/1‑148 (MQ=255)
               cATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGACCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTg                                                                                                                                   >  1:871815/1‑148 (MQ=255)
                 tGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAg                                                                                                                                 >  1:3537901/1‑148 (MQ=255)
                  gCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAgg                                                                                                                                >  1:1621558/1‑148 (MQ=255)
                     ggATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCaa                                                                                                                             >  1:319458/1‑148 (MQ=255)
                         tAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGa                                                                                                                         <  2:3452046/148‑1 (MQ=255)
                          aGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGaa                                                                                                                        <  1:367196/148‑1 (MQ=255)
                          aGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATAGCCTGGCGTGGCTGCGTGAGGCAAATGaa                                                                                                                        >  1:1861997/1‑148 (MQ=255)
                           ggCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAAc                                                                                                                       <  2:2985555/148‑1 (MQ=255)
                           ggCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAAc                                                                                                                       <  1:1354948/148‑1 (MQ=255)
                           ggCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAAc                                                                                                                       >  2:2149961/1‑148 (MQ=255)
                           ggCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAAc                                                                                                                       >  2:3430611/1‑148 (MQ=255)
                            gCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACa                                                                                                                      <  2:505779/148‑1 (MQ=255)
                            gCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACa                                                                                                                      >  1:3290380/1‑148 (MQ=255)
                            gCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACa                                                                                                                      <  2:505096/148‑1 (MQ=255)
                             cGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAg                                                                                                                     <  2:3078070/148‑1 (MQ=255)
                              ggTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGt                                                                                                                    <  1:2585347/148‑1 (MQ=255)
                                 gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCg                                                                                                                 <  2:3234530/148‑1 (MQ=255)
                                  cACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCga                                                                                                                <  1:1715279/148‑1 (MQ=255)
                                   aCGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGATGCGTGAGGCAAATGAACAGTTCgat                                                                                                               <  1:3012551/148‑1 (MQ=255)
                                    cGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCgatc                                                                                                              <  2:3270278/148‑1 (MQ=255)
                                    cGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCgatc                                                                                                              <  2:1508445/148‑1 (MQ=255)
                                    cGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCgatc                                                                                                              <  1:1695155/148‑1 (MQ=255)
                                    cGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCgatc                                                                                                              >  2:2889899/1‑148 (MQ=255)
                                           ccaccaccGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCtt                                                                                                       <  1:2202240/148‑1 (MQ=255)
                                            caccaccGTGGATATGTCGCGTACTTATCTGGAATGGGCTGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCttt                                                                                                      <  2:746076/148‑1 (MQ=255)
                                            caccaccGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCttt                                                                                                      <  2:2527639/148‑1 (MQ=255)
                                             accaccGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTa                                                                                                     <  2:3354781/148‑1 (MQ=255)
                                               caccGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATAGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATc                                                                                                   <  1:126113/148‑1 (MQ=255)
                                                accGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCg                                                                                                  <  1:2890732/148‑1 (MQ=255)
                                                accGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCg                                                                                                  >  2:140808/1‑148 (MQ=255)
                                                 ccGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGa                                                                                                 >  2:182877/1‑148 (MQ=255)
                                                 ccGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGa                                                                                                 >  2:3425406/1‑148 (MQ=255)
                                                  cGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGAt                                                                                                >  2:3160503/1‑148 (MQ=255)
                                                  cGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGAt                                                                                                <  2:2302300/148‑1 (MQ=255)
                                                     ggATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATcag                                                                                             >  1:3312820/1‑146 (MQ=255)
                                                     ggATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCg                                                                                             >  1:1906441/1‑148 (MQ=255)
                                                               cgTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTtct                                                                                   <  2:1706231/148‑1 (MQ=255)
                                                                gTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTtctc                                                                                  <  1:968978/148‑1 (MQ=255)
                                                                 tACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTctct                                                                                 >  2:200630/1‑148 (MQ=255)
                                                                  aCTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTa                                                                                >  1:2806917/1‑148 (MQ=255)
                                                                   cTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATAGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTaa                                                                               <  2:2392876/148‑1 (MQ=255)
                                                                          ggAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCaaaa                                                                        <  2:2203750/148‑1 (MQ=255)
                                                                          ggAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCaaaa                                                                        <  1:2335144/148‑1 (MQ=255)
                                                                          ggAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCaaaa                                                                        <  1:2747409/148‑1 (MQ=255)
                                                                           gAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAAc                                                                       >  1:3482995/1‑148 (MQ=255)
                                                                           gAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAAc                                                                       >  2:3448275/1‑148 (MQ=255)
                                                                             aTGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGa                                                                     <  2:1863761/148‑1 (MQ=255)
                                                                              tGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGaa                                                                    <  2:2177848/148‑1 (MQ=255)
                                                                               gggCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAAt                                                                   <  1:1130778/148‑1 (MQ=255)
                                                                                   aGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGaa                                                               >  1:713518/1‑148 (MQ=255)
                                                                                    gAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAg                                                              <  1:1420425/148‑1 (MQ=255)
                                                                                    gAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAg                                                              >  1:2833897/1‑148 (MQ=255)
                                                                                     aaCGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGa                                                             >  1:3328352/1‑148 (MQ=255)
                                                                                     aaCGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATAGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGa                                                             >  2:1642221/1‑148 (MQ=255)
                                                                                        gCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGc                                                          >  1:1879915/1‑148 (MQ=255)
                                                                                        gCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGc                                                          >  1:598484/1‑148 (MQ=255)
                                                                                         cAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCg                                                         >  1:2463558/1‑148 (MQ=255)
                                                                                            ccTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGttt                                                      >  2:79595/1‑148 (MQ=255)
                                                                                                  tCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGtt                                                >  1:1795847/1‑148 (MQ=255)
                                                                                                  tCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGtt                                                <  2:545491/148‑1 (MQ=255)
                                                                                                   cTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTc                                               <  2:1304759/148‑1 (MQ=255)
                                                                                                           ccTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGat                                       <  2:3328352/148‑1 (MQ=255)
                                                                                                           ccTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGat                                       >  2:4796/1‑148 (MQ=255)
                                                                                                           ccTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGat                                       <  2:1266246/148‑1 (MQ=255)
                                                                                                            cTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCCGCGCGatc                                      >  2:3156838/1‑148 (MQ=255)
                                                                                                            cTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGatc                                      <  2:108365/148‑1 (MQ=255)
                                                                                                            cTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGatc                                      >  2:1865584/1‑148 (MQ=255)
                                                                                                              gACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGatcat                                    <  1:1362938/148‑1 (MQ=255)
                                                                                                               aCCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGatcatc                                   >  2:3487805/1‑148 (MQ=255)
                                                                                                                ccGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCt                                  <  2:520605/148‑1 (MQ=255)
                                                                                                                 cGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTg                                 >  1:1003425/1‑148 (MQ=255)
                                                                                                                       tGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCtg                           >  1:227825/1‑148 (MQ=255)
                                                                                                                        gcgcATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCtga                          <  2:2139622/148‑1 (MQ=255)
                                                                                                                         cgcATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCtgat                         >  2:2689237/1‑148 (MQ=255)
                                                                                                                            aTCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGaa                      <  2:70443/148‑1 (MQ=255)
                                                                                                                             tCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGaaa                     >  2:2790891/1‑148 (MQ=255)
                                                                                                                             tCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGaaa                     >  2:2461411/1‑148 (MQ=255)
                                                                                                                               gCCTGATTCAGGCCGATGGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGa                   <  2:383384/148‑1 (MQ=255)
                                                                                                                                     ttCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAATAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGaa             <  1:3514170/148‑1 (MQ=255)
                                                                                                                                     ttCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGaa             <  2:3537835/148‑1 (MQ=255)
                                                                                                                                      tCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGaaa            >  2:1588966/1‑148 (MQ=255)
                                                                                                                                         ggCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGAAACGt         <  2:2124175/148‑1 (MQ=255)
                                                                                                                                         ggCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGAAACGt         >  1:2053672/1‑148 (MQ=255)
                                                                                                                                         ggCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGAAACGt         <  2:1141413/148‑1 (MQ=255)
                                                                                                                                           ccGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGAAACGTct       <  2:871815/148‑1 (MQ=255)
                                                                                                                                              aTTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGAAACGTctgct    >  2:2675620/1‑148 (MQ=255)
                                                                                                                                               ttGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGAAACGTctgctg   >  1:2017454/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGCAGCGCCACCGTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTGCCTGGCGTGGCTGCGTGAGGCAAATGAACAGTTCGATCTGATCTTTATCGATCCGCCAACCTTCTCTAACTCAAAACGAATGGAAGATGCGTTTGATGTTCAGCGCGATCATCTGGCGCTGATGAAAGATTTGAAACGTCTGCTGC  >  CP009273/1004947‑1005238

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: