Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,641,274 G→C 5.5% T139T (ACC→ACG) ‡ cmoA ← tRNA 5‑methoxyuridine(34)/uridine‑5‑oxyacetic acid(34) metyltransferase CmoA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,641,2740GC 5.5% 225.8 / 4.6 73T139T (ACC→ACG) ‡cmoAtRNA 5‑methoxyuridine(34)/uridine‑5‑oxyacetic acid(34) metyltransferase CmoA
Reads supporting (aligned to +/- strand):  ref base G (29/40);  new base C (2/2);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

AATCCAGGTGCAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTGGCGGCAGCGCTCGACCATGGCGGCGGAGTT  >  NC_002947/1641126‑1641420
                                                                                                                                                    |                                                                                                                                                  
aaTCCCGGTGCAGTTCATTGAGCAGGTCCTGCTCTTGGTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAgg                                                                                                                                                    >  2:253962/1‑149 (MQ=255)
  tCCAGGTGCAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGt                                                                                                                                                  <  2:188988/149‑1 (MQ=255)
    cAGGTGCAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGa                                                                                                                                                <  5:50388/149‑1 (MQ=255)
     aGGTGCAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTg                                                                                                                                                 >  4:1236/1‑147 (MQ=255)
      ggTGCAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAg                                                                                                                                              <  3:1236/149‑1 (MQ=255)
       gTGCAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGt                                                                                                                                             <  4:75202/149‑1 (MQ=255)
       gTGCAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGt                                                                                                                                             <  8:309492/149‑1 (MQ=255)
        tGCAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGtt                                                                                                                                            >  2:128105/1‑149 (MQ=255)
         gCAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGtt                                                                                                                                            >  2:343314/1‑148 (MQ=255)
          cAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCa                                                                                                                                          >  1:93817/1‑149 (MQ=255)
             ttCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCCGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGGGAAGTTCATCg                                                                                                                                       <  1:253962/149‑1 (MQ=255)
             ttCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCg                                                                                                                                       >  3:311428/1‑149 (MQ=255)
                aTTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCa                                                                                                                                    >  2:102715/1‑149 (MQ=255)
                aTTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCa                                                                                                                                    >  3:19169/1‑149 (MQ=255)
                 ttGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCAc                                                                                                                                   >  1:160882/1‑149 (MQ=255)
                    aGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTa                                                                                                                                 <  5:164682/148‑1 (MQ=255)
                         gTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCACCATCGCCTGGCGGATTCGCCCGACCAGTTCTAGGCGCTGGTCAGGGGCCATGAACTGCAGGGTGAAGTTCATCGCCACTACCGacg                                                                                                                           >  2:312915/1‑147 (MQ=255)
                               ctctTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCt                                                                                                                     >  5:94643/1‑149 (MQ=255)
                                tctTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAg                                                                                                                            <  5:191178/141‑1 (MQ=255)
                                tctTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAg                                                                                                                            >  6:191178/1‑141 (MQ=255)
                                tctTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTc                                                                                                                    >  2:239866/1‑149 (MQ=255)
                                     ttCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAagg                                                                                                               <  1:102715/149‑1 (MQ=255)
                                         tCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTa                                                                                                                                 >  6:134027/1‑127 (MQ=255)
                                         tCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTa                                                                                                                                 <  5:134027/127‑1 (MQ=255)
                                         tCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAagggcag                                                                                                           <  5:164317/149‑1 (MQ=255)
                                           ggCGATGCGCAGCTTCTCCGACAGAATCCGCGCGCGGGGGGGCAGCAGCGCCTGGGGTATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAAGTGATGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAagggcagg                                                                                                          <  2:160882/148‑1 (MQ=255)
                                           ggCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAagggcaggg                                                                                                         >  6:15969/1‑149 (MQ=255)
                                            gCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAagggcagggc                                                                                                        <  1:239866/149‑1 (MQ=255)
                                              gAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCaa                                                                                                      <  8:10665/149‑1 (MQ=255)
                                                 gcgcAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGa                                                                                                   <  1:128105/149‑1 (MQ=255)
                                                     aGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATc                                                                                               >  1:352834/1‑149 (MQ=255)
                                                     aGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATc                                                                                               <  1:82470/149‑1 (MQ=255)
                                                     aGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATc                                                                                               <  4:36204/149‑1 (MQ=255)
                                                       cTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCgg                                                                                             >  8:220118/1‑149 (MQ=255)
                                                          ctcCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCtt                                                                                          >  1:45301/1‑149 (MQ=255)
                                                             ggcttcaaTCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCa                                                                                       <  1:190162/143‑1 (MQ=255)
                                                             cGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCa                                                                                       <  2:93817/149‑1 (MQ=255)
                                                             cGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCa                                                                                       <  1:251724/149‑1 (MQ=255)
                                                              gACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTTCACCCTGCAGTTCATCGc                                                                                                                                      <  8:92164/101‑1 (MQ=38)
                                                              gACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTTCACCCTGCAGTTCATCGc                                                                                                                                      >  7:92164/1‑101 (MQ=38)
                                                              gACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAg                                                                                      >  3:153887/1‑149 (MQ=255)
                                                              gACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGACTCGAAGGGCAGGGCAAGGATATCGGCTTcccg                                                                                      >  3:5823/1‑147 (MQ=255)
                                                                cAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGc                                                                                                                      <  1:324909/115‑1 (MQ=255)
                                                                cAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGc                                                                                                                      >  2:324909/1‑115 (MQ=255)
                                                                   aaTCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACtt                                                                                 <  4:311428/149‑1 (MQ=255)
                                                                     tCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTg                                                                                >  1:54674/1‑148 (MQ=255)
                                                                          gcgcCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTgg                                                                          <  7:220118/149‑1 (MQ=255)
                                                                              cgccgGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGgcag                                                                       >  7:77974/1‑148 (MQ=255)
                                                                                     cagcagCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCt                                                               <  2:45301/149‑1 (MQ=255)
                                                                                      agcagcGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTg                                                              <  6:94643/149‑1 (MQ=255)
                                                                                       gcagcGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTgg                                                             <  2:54674/149‑1 (MQ=255)
                                                                                             ccTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATc                                                       <  4:19169/149‑1 (MQ=255)
                                                                                                  cGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATc                                                  <  5:15969/149‑1 (MQ=255)
                                                                                                    gATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCt                                                <  1:343314/149‑1 (MQ=255)
                                                                                                         gCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCggg                                                                                                                       <  4:236482/73‑1 (MQ=255)
                                                                                                         gCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCggg                                                                                                                       >  3:236482/1‑73 (MQ=255)
                                                                                                         gCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTgcg                                             >  7:69983/1‑147 (MQ=255)
                                                                                                         gCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTgcg                                             >  2:248069/1‑147 (MQ=255)
                                                                                                              aGCTGTTCTCGGCGCTGGTAAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAg                                      <  6:210620/149‑1 (MQ=255)
                                                                                                               gCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAgg                                     >  1:353609/1‑149 (MQ=255)
                                                                                                                          cgcTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAagggcagg                                                                                                          >  4:264902/1‑69 (MQ=255)
                                                                                                                          cgcTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAagggcagg                                                                                                          <  3:264902/69‑1 (MQ=255)
                                                                                                                                tCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTGGCGGCAGCGCTc                    >  3:261674/1‑149 (MQ=255)
                                                                                                                                  aGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTGGCGGCAGCGCTCg                   <  2:24382/148‑1 (MQ=255)
                                                                                                                                  aGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTGGCGGCAGCGCTCGa                  <  4:153887/149‑1 (MQ=255)
                                                                                                                                  aGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTGGCGGCAGCGCTCGa                  <  2:128109/149‑1 (MQ=255)
                                                                                                                                     ggCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTggcg                             <  4:264571/135‑1 (MQ=255)
                                                                                                                                     ggCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTggcg                             >  3:264571/1‑135 (MQ=255)
                                                                                                                                        gATGAACTTCACCCTGCAGTTCATCGCCCCTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTgc                                              >  4:291607/1‑115 (MQ=38)
                                                                                                                                        gATGAACTTCACCCTGCAGTTCATCGCCCCTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTgc                                              <  3:291607/115‑1 (MQ=38)
                                                                                                                                            aaCTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTGGCGGCAGCGCTCGACCATggcgg         <  8:77974/148‑1 (MQ=255)
                                                                                                                                                 cAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTGGCGGCAGCGCTCGACCATGGCGGCGGa     <  8:69983/147‑1 (MQ=255)
                                                                                                                                                 cAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTGGCGGCAGCGCTCGACCATGGCGGCGGAGt   <  4:261674/149‑1 (MQ=255)
                                                                                                                                                  aGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTGGCGGCAGCGCTCGACCATGGCGGCGGAGtt  <  1:326516/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATCCAGGTGCAGTTCATTGAGCAGGTCCTGCTCTTGTTCATCGGCGAAGCGCAGCTTCTCCGACAGAATCAGCGCGCCGCCGGGCAGCAGCGCCTGGCGGATTCGCCCGAGCAGTTCTAGGCGCTGGTCAGGGGCGATGAACTGCAGGGTGAAGTTCATCGCCACTACCGAGGCGGGCTCGAAGGGCAGGGCAAGGATATCGGCTTCCAGTACTTGCACTGGCAGCAACTCCTGGAACATCGAATCCTGCGCAGTGAGGTACTGGCGGCAGCGCTCGACCATGGCGGCGGAGTT  >  NC_002947/1641126‑1641420

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: