Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 604,470 C→G 5.7% P101R (CCA→CGA) ‡ nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947604,4700CG 5.7% 225.4 / 3.9 70P101R (CCA→CGA) ‡nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base C (29/37);  new base G (2/2);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTG  >  NC_002947/604325‑604615
                                                                                                                                                 |                                                                                                                                                 
gAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGt                                                                                                                                                >  6:303650/1‑149 (MQ=255)
  aGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTc                                                                                                                                               >  2:193186/1‑148 (MQ=255)
   gcgcAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGa                                                                                                                                             <  4:165783/149‑1 (MQ=255)
      cAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAgct                                                                                                                                          <  1:13866/149‑1 (MQ=255)
            cGGGTGGAAAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAATCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctcttcca                                                                                                                                    <  5:264562/149‑7 (MQ=255)
                  gATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGtgctgc                                                                                                                               <  1:193186/148‑1 (MQ=255)
                  gATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGtgctgc                                                                                                                               <  7:61175/148‑1 (MQ=255)
                   aTAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTgcgc                                                                                                                             >  2:340920/1‑149 (MQ=255)
                   aTAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTgcgc                                                                                                                             >  5:145002/1‑149 (MQ=255)
                    tAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgct                                                                                                                                          >  6:2800/1‑135 (MQ=255)
                    tAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgct                                                                                                                                          <  5:2800/135‑1 (MQ=255)
                    tAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCAACCCAGTCGAgct                                                                                                                                          <  8:204169/135‑1 (MQ=255)
                    tAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCAACCCAGTCGAgct                                                                                                                                          >  7:204169/1‑135 (MQ=255)
                      aCGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTg                                                                                                                          <  6:291160/149‑1 (MQ=255)
                        gagttCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTc                                                                                                                        <  3:5552/146‑1 (MQ=255)
                        gATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTc                                                                                                                        <  1:340920/149‑1 (MQ=255)
                              aCCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTg                                                                                                                  >  1:16425/1‑149 (MQ=255)
                               ccGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCAGGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTgg                                                                                                                 <  6:315503/149‑1 (MQ=255)
                                   tggcctcGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATGGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGa                                                                                                               <  4:253794/141‑1 (MQ=255)
                                     tCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCt                                                                                                                           <  8:216026/133‑1 (MQ=255)
                                     tCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCt                                                                                                                           >  7:216026/1‑133 (MQ=255)
                                     tCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTg                                                                                                           <  6:145002/149‑1 (MQ=255)
                                       gACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGACTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGa                                                                                                          >  5:36341/1‑148 (MQ=255)
                                        aCGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGAt                                                                                                         >  1:292243/1‑148 (MQ=255)
                                        aCGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATc                                                                                                        <  5:303650/149‑1 (MQ=255)
                                           gTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCaaa                                                                                                     <  5:195565/149‑1 (MQ=255)
                                            tGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAc                                                                                                    >  7:347330/1‑149 (MQ=255)
                                             gCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACTAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCTCCTGTCAACCTGGGAGCTGATCAAAcg                                                                                                   >  3:37184/1‑149 (MQ=255)
                                               cTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAcgcg                                                                                                 >  1:297220/1‑149 (MQ=255)
                                                  tttCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCg                                                                                              <  3:203583/149‑1 (MQ=255)
                                                        cgAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCGGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGc                                                                                        <  1:211874/149‑1 (MQ=255)
                                                                     ggggTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGcgtcgt                                                                           >  8:180958/1‑149 (MQ=255)
                                                                         tCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGAt                                                                        >  6:231065/1‑148 (MQ=255)
                                                                               cAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGaa                                                                 <  2:297220/149‑1 (MQ=255)
                                                                                aaTCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAAc                                                                    <  6:122283/145‑1 (MQ=255)
                                                                                aaTCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAAc                                                                    >  5:122283/1‑145 (MQ=255)
                                                                                  tCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGgt                                                              <  8:347330/149‑1 (MQ=255)
                                                                                     aGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGATCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGtgtg                                                           <  5:26328/149‑1 (MQ=255)
                                                                                          gAAATCGACAAAGCGCTGGTACCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACt                                                      >  4:62938/1‑149 (MQ=255)
                                                                                           aaaTCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTGCCGCGTCGTGATCAACGAAGGTGTGGAACTg                                                     <  4:176717/149‑1 (MQ=255)
                                                                                            aaTCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTgg                                                    <  2:292243/149‑1 (MQ=255)
                                                                                             aTCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGc                                                   >  2:128727/1‑149 (MQ=255)
                                                                                               cGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGcc                                                  >  6:333801/1‑148 (MQ=255)
                                                                                                gACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCaa                                                <  1:270491/149‑1 (MQ=255)
                                                                                                   aaaGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGAt                                                                                                         >  4:73428/1‑89 (MQ=255)
                                                                                                   aaaGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGAt                                                                                                         <  3:73428/89‑1 (MQ=255)
                                                                                                    aaGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTc                                            >  5:222803/1‑149 (MQ=255)
                                                                                                     aGCTCTGTTGCATTTCATGACCATCGCACTGTAAGAGCTCGACACAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCt                                           <  4:37184/149‑1 (MQ=255)
                                                                                                       cgcTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTc                                         <  2:16425/149‑1 (MQ=255)
                                                                                                              gCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCa                                  >  7:295311/1‑149 (MQ=255)
                                                                                                               ccTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCa                                  <  7:180958/148‑1 (MQ=255)
                                                                                                               ccTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCAc                                 >  3:228847/1‑149 (MQ=255)
                                                                                                                   gCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCa                                                                      <  6:265790/108‑1 (MQ=255)
                                                                                                                   gCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCa                                                                      >  5:265790/1‑108 (MQ=255)
                                                                                                                    cATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGAc                             <  8:295311/148‑1 (MQ=255)
                                                                                                                     aTGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGTCCAAGGTCTTCGGCGCCTCTGACgg                           <  6:36341/149‑1 (MQ=255)
                                                                                                                     aTGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACgg                           >  4:268356/1‑149 (MQ=255)
                                                                                                                        aCCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGcc                         <  5:333801/148‑1 (MQ=255)
                                                                                                                         ccATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCcac                       <  3:62938/149‑1 (MQ=255)
                                                                                                                          cATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCaa                                                                     <  3:300563/102‑1 (MQ=255)
                                                                                                                          cATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCaa                                                                     >  4:300563/1‑102 (MQ=255)
                                                                                                                            tCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTgg                                                                                                                 <  4:11757/56‑1 (MQ=255)
                                                                                                                            tCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTgg                                                                                                                 >  3:11757/1‑56 (MQ=255)
                                                                                                                               cACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTc                 <  5:231065/149‑1 (MQ=255)
                                                                                                                                    gAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCaa            >  8:107835/1‑149 (MQ=255)
                                                                                                                                         gCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCg       >  7:262671/1‑149 (MQ=255)
                                                                                                                                             gACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAgg                                                               <  2:80729/89‑1 (MQ=255)
                                                                                                                                             gACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAgg                                                               >  1:80729/1‑89 (MQ=255)
                                                                                                                                             gACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGTCCACAAGTTCGTCAACGGCGTACt   >  3:40633/1‑149 (MQ=255)
                                                                                                                                              aCCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTg  <  7:107835/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTG  >  NC_002947/604325‑604615

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: