Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 539,942 C→A 7.2% P774T (CCG→ACG)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947539,9420CA 7.2% 219.7 / 6.2 69P774T (CCG→ACG) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base A (3/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

CGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCT  >  NC_002947/539795‑540086
                                                                                                                                                   |                                                                                                                                                
cgcgGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTcc                                                                                                                                                 >  7:6648/1‑149 (MQ=255)
 gcgGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCg                                                                                                                                                >  6:65156/1‑149 (MQ=255)
      ggtgTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACATGACCAAGTACACCAGTTCGAACCAGAACACCTTCATCAACCAGCGTCCGCTGGt                                                                                                                                           >  5:89211/1‑149 (MQ=255)
       gtgtGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGt                                                                                                                                           <  8:73364/148‑1 (MQ=255)
         gtgATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAg                                                                                                                                        >  4:30371/1‑149 (MQ=255)
              cGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAgg                                                                                                                                   >  8:249577/1‑149 (MQ=255)
               gACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCAGGTgcggcagt                                                                                                                                   <  8:6648/148‑8 (MQ=255)
                      tCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAg                                                                                                                           >  4:68578/1‑149 (MQ=255)
                        gACGCCAGCCGTATCGTTGTTCTCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGt                                                                                                                         >  7:227824/1‑149 (MQ=255)
                        gACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGt                                                                                                                                           >  1:224613/1‑131 (MQ=255)
                        gACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGt                                                                                                                                           <  2:224613/131‑1 (MQ=255)
                        gACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAgg                                                                                                                          <  5:65156/148‑1 (MQ=255)
                           gCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGtt                                                                                                                        >  8:7862/1‑147 (MQ=255)
                             cAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGc                                                                                                                    <  6:57087/149‑1 (MQ=255)
                                  gTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                               <  5:184447/149‑1 (MQ=255)
                                   tATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgtggt                                                                                                               >  1:186900/1‑148 (MQ=255)
                                     tCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGTGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGa                                                                                                                                         <  6:35847/120‑1 (MQ=255)
                                     tCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGa                                                                                                             >  2:311185/1‑148 (MQ=255)
                                           ttCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGTGTGACAAGGTTCAGCGTGGTGACATCATg                                                                                                      >  3:228205/1‑149 (MQ=255)
                                             cgcgTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATgg                                                                                                     <  7:249577/148‑1 (MQ=255)
                                             cgcgTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGc                                                                                                    >  5:2655/1‑149 (MQ=255)
                                              gcgTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGc                                                                                                    >  8:218938/1‑148 (MQ=255)
                                              gcgTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGc                                                                                                    <  7:218938/148‑1 (MQ=255)
                                                    ccGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGc                                                                                             <  2:234418/149‑1 (MQ=255)
                                                      gacgacgaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGccc                                                                                           <  3:30371/149‑1 (MQ=255)
                                                      gacgacgaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGccc                                                                                           >  8:248578/1‑149 (MQ=255)
                                                         gacgaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgt                                                                                                                  <  8:58964/123‑1 (MQ=255)
                                                         gacgaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCgt                                                                                                                  >  7:58964/1‑123 (MQ=255)
                                                         gacgaAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTc                                                                                        >  8:191887/1‑149 (MQ=255)
                                                                  gAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCTAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACAt                                                                               <  7:7862/149‑1 (MQ=255)
                                                                  gAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACAt                                                                               >  2:38397/1‑149 (MQ=255)
                                                                    aaCCTGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGGTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATgg                                                                             <  7:318647/149‑1 (MQ=255)
                                                                    aaCCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATgg                                                                             <  2:186900/149‑1 (MQ=255)
                                                                    aaCCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATgg                                                                             <  1:311185/149‑1 (MQ=255)
                                                                     aCCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATggg                                                                            <  7:191887/149‑1 (MQ=255)
                                                                                   tgtgGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTgg                                                              <  3:68578/149‑1 (MQ=255)
                                                                                             tACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACAt                                                    >  5:39664/1‑149 (MQ=255)
                                                                                               cAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAAc                                                      <  1:38397/145‑1 (MQ=255)
                                                                                               cAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATgc                                                  >  4:232187/1‑149 (MQ=255)
                                                                                                       ccAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATcgcg                                          >  6:293401/1‑149 (MQ=255)
                                                                                                         aaGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGtt                                        >  1:159068/1‑149 (MQ=255)
                                                                                                           gTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCa                                      >  3:124713/1‑149 (MQ=255)
                                                                                                              cacCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATgg                                   >  5:42018/1‑149 (MQ=255)
                                                                                                              cacCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATgg                                   >  2:271773/1‑149 (MQ=255)
                                                                                                                 ccGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATgg                                                                             >  5:58123/1‑104 (MQ=255)
                                                                                                                 ccGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACAt                                                    <  6:327634/129‑1 (MQ=255)
                                                                                                                 ccGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACAt                                                    >  5:327634/1‑129 (MQ=255)
                                                                                                                  cGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTg                               <  2:159068/149‑1 (MQ=255)
                                                                                                                       gAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAAc                           <  1:271773/148‑1 (MQ=255)
                                                                                                                           cAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCtt                      <  7:248578/149‑1 (MQ=255)
                                                                                                                             gAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTc                     >  1:314568/1‑148 (MQ=255)
                                                                                                                             gAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTc                     >  7:128092/1‑148 (MQ=255)
                                                                                                                                 acCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCtcgaaat                 >  5:108682/1‑142 (MQ=255)
                                                                                                                                     gCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCg              >  3:40343/1‑147 (MQ=255)
                                                                                                                                      cATCAACCAGCTTCCGCTGGTGAGCAAGTGTGACAAGGTTCAGCGTGTTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGTGTCAGAACATGCGCATAGCGTTCATTGCGTGTAACGGCTTCAACTTCGaa            <  6:89211/148‑1 (MQ=255)
                                                                                                                                      cATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGaa            <  3:232187/148‑1 (MQ=255)
                                                                                                                                        tCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGaa            >  2:282491/1‑146 (MQ=255)
                                                                                                                                        tCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGaa            <  1:282491/146‑1 (MQ=255)
                                                                                                                                         tcaCCAGCGGACGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATgg                                                                             <  5:23125/78‑1 (MQ=255)
                                                                                                                                         tcaCCAGCGGACGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATgg                                                                             >  6:23125/3‑80 (MQ=255)
                                                                                                                                         tcaCCAGCGGACGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGaa            >  7:26543/3‑145 (MQ=255)
                                                                                                                                         tcaCCAGCGGACGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGaa            <  8:26543/143‑1 (MQ=255)
                                                                                                                                         tcaCCAGCGGACGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACt        >  1:39515/3‑149 (MQ=255)
                                                                                                                                           aCCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTc       >  5:196008/1‑148 (MQ=255)
                                                                                                                                             cAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTc       <  8:78706/146‑1 (MQ=255)
                                                                                                                                             cAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTc       >  7:78706/1‑146 (MQ=255)
                                                                                                                                             cAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCAt    >  7:15794/1‑149 (MQ=255)
                                                                                                                                              aGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATc   <  4:40343/149‑1 (MQ=255)
                                                                                                                                               gCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACCTGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCt  <  8:15794/149‑1 (MQ=255)
                                                                                                                                               gCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCt  <  8:128092/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCCGACGACGAAGTGGAAACCGGCGAAGCGGGTGTGGATATCTACAACCTGACCAAGTACACCCGTTCGAACCAGAACACCTGCATCAACCAGCGTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCT  >  NC_002947/539795‑540086

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: ‑