Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,535,236 C→T 6.8% P85S (CCC→TCC)  pncB → nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,535,2360CT 6.8% 219.6 / 4.0 73P85S (CCC→TCC) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base T (3/2);  total (41/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.47e-01

CCGACGTCGAATGGGAATTCCGCTGCCGCAACGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAG  >  NC_002947/5535090‑5535381
                                                                                                                                                  |                                                                                                                                                  
ccGACGTCGAATGGGAATTCCGCTGCCGCAACGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGccc                                                                                                                                                  <  4:39842/149‑1 (MQ=255)
        gAATGGGAATTCCGCTGCCGCAACGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACtttt                                                                                                                                           <  3:92062/148‑1 (MQ=255)
            gggAATTCCGCTGCCGCAACGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGc                                                                                                                                      >  6:174343/1‑149 (MQ=255)
              gAATTCCGCTGCCGCAACGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCtt                                                                                                                                    >  7:13139/1‑149 (MQ=255)
                     gCTGCCGCAACGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTAc                                                                                                                                        <  4:17789/138‑1 (MQ=255)
                     gCTGCCGCAACGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTAc                                                                                                                                        >  3:17789/1‑138 (MQ=255)
                          cGCAACGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGAGGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGtt                                                                                                                       >  1:99186/1‑149 (MQ=255)
                             aaCGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCg                                                                                                                    >  2:8656/1‑149 (MQ=255)
                             aaCGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCg                                                                                                                    >  2:3587/1‑149 (MQ=255)
                              aCGGTGAGGACCTGCGCCCCTACCTCGGCGAAGTTCGCAATCAGCTCGAACGGCTTGCCGACATCACCCTGGACGATGGCCGGCTGGCCGTCCTCGAACGCAGCAGCGTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCGGtgacgg                                                                                                                   >  4:49925/1‑144 (MQ=255)
                                  tGAGGACCTGCGCCACTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTAGCCGACCTCACCCAGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCTCCTGTTCCGCTTc                                                                                                                >  4:215040/1‑148 (MQ=255)
                                         cTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTCTCTGCGCCTGTTCCGCTTCAATCTGc                                                                                                         <  1:8656/149‑1 (MQ=255)
                                         cTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTCTCTGCGCCTGTTCCGCTTCAATCTGc                                                                                                         <  1:3587/149‑1 (MQ=255)
                                         cTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGc                                                                                                         <  7:185880/148‑1 (MQ=255)
                                         cTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGc                                                                                                         >  8:185880/1‑148 (MQ=255)
                                         cTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCg                                                                                                        >  2:270418/1‑149 (MQ=255)
                                            cgcCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGtt                                                                                                      <  8:13139/148‑1 (MQ=255)
                                                        ggCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGAGGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCgg                                                                                         <  2:99186/149‑1 (MQ=255)
                                                        ggCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCgg                                                                                         >  3:305891/1‑149 (MQ=255)
                                                        ggCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCg                                                                                          <  5:174343/148‑1 (MQ=255)
                                                         gCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTgc                                                                                                                              >  5:148849/1‑111 (MQ=255)
                                                         gCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTgc                                                                                                                              <  6:148849/111‑1 (MQ=255)
                                                           gAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACtttt                                                                                                                                           <  8:217487/97‑1 (MQ=255)
                                                           gAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACtttt                                                                                                                                           >  7:217487/1‑97 (MQ=255)
                                                              aTTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCg                                                                                    >  1:68874/1‑148 (MQ=255)
                                                               ttCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCgg                                                                                         <  5:18819/142‑1 (MQ=255)
                                                               ttCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCgg                                                                                         >  6:18819/1‑142 (MQ=255)
                                                               ttCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGa                                                                                   >  7:3024/1‑148 (MQ=255)
                                                                   cAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATg                                                                              <  1:5600/149‑1 (MQ=255)
                                                                      tCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGACTGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGAt                                                                            >  1:249669/1‑148 (MQ=255)
                                                                         gCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGc                                                                        >  2:283845/1‑149 (MQ=255)
                                                                             gAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGt                                                                     >  4:216385/1‑148 (MQ=255)
                                                                             gAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGAGCACGTCGGCATCGAAAATGATCAGCTGtt                                                                    >  7:101826/1‑149 (MQ=255)
                                                                                 ggCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGc                                                                                                         >  1:104579/1‑108 (MQ=255)
                                                                                 ggCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGc                                                                                                         <  2:104579/108‑1 (MQ=255)
                                                                                  gCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCTAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTgc                                                               <  3:216385/149‑1 (MQ=255)
                                                                                  gCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACCCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTgc                                                               >  1:100221/1‑149 (MQ=255)
                                                                                               aCCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTc                                                   <  8:207182/148‑1 (MQ=255)
                                                                                               aCCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGAGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTcc                                                  <  8:101826/149‑1 (MQ=255)
                                                                                               aCCCAGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTc                                                   <  1:270418/148‑1 (MQ=255)
                                                                                                 ccTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCGGCTGTTCCTGCGCCTGAAAGATCCGt                                                <  3:215040/149‑1 (MQ=255)
                                                                                                 ccTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGt                                                >  2:171745/1‑149 (MQ=255)
                                                                                                 ccTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGt                                                <  2:68874/149‑1 (MQ=255)
                                                                                                 ccTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGt                                                >  7:228785/1‑149 (MQ=255)
                                                                                                  cTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGaaaa                                                                                >  1:32314/1‑116 (MQ=255)
                                                                                                  cTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGaaaa                                                                                <  2:32314/116‑1 (MQ=255)
                                                                                                  cTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTg                                               >  1:139538/1‑149 (MQ=255)
                                                                                                     gACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTgg                                              <  2:139538/147‑1 (MQ=255)
                                                                                                      aCGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTg                                           >  8:274952/1‑149 (MQ=255)
                                                                                                       cGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTg                                           >  3:20308/1‑148 (MQ=255)
                                                                                                       cGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGc                                          >  7:177628/1‑149 (MQ=255)
                                                                                                        gATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCa                                         <  8:3024/149‑1 (MQ=255)
                                                                                                                         ccTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGaa                         <  4:305891/148‑1 (MQ=255)
                                                                                                                         ccTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAg                        >  6:84089/1‑149 (MQ=255)
                                                                                                                           tCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATc                                  <  6:240168/137‑1 (MQ=255)
                                                                                                                           tCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATc                                  >  5:240168/1‑137 (MQ=255)
                                                                                                                           tCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTg                      >  8:48781/1‑149 (MQ=255)
                                                                                                                             gAACGCATCAGCTTCCTCAAGTCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGaaa                                                       >  4:25877/1‑114 (MQ=255)
                                                                                                                             gAACGCATCAGCTTCCTCAAGTCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGaaa                                                       <  3:25877/114‑1 (MQ=255)
                                                                                                                             gAACGCATCAGCTTCCTCAAGTCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGaaa                                                       <  5:100472/114‑1 (MQ=255)
                                                                                                                             gAACGCATCAGCTTCCTCAAGTCCGACTTTTTACGCTTTCTGCGCCTGTTCCGATTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGaaa                                                       >  6:100472/1‑114 (MQ=255)
                                                                                                                              aaCGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCg                   <  7:274952/149‑1 (MQ=255)
                                                                                                                                   aTCAGCTTCCTCAAGTCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCgctgct              >  3:51740/1‑149 (MQ=255)
                                                                                                                                   aTCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCgctgct              >  3:137441/1‑149 (MQ=255)
                                                                                                                                       gCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGcc          <  2:249669/149‑1 (MQ=255)
                                                                                                                                           ccTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATc                                                                           >  6:162207/1‑80 (MQ=255)
                                                                                                                                           ccTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATc                                                                           <  5:162207/80‑1 (MQ=255)
                                                                                                                                            cTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATc                                  <  8:181303/120‑1 (MQ=255)
                                                                                                                                            cTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATc                                  >  7:181303/1‑120 (MQ=255)
                                                                                                                                             tCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCcatcat     >  8:31579/1‑148 (MQ=255)
                                                                                                                                             tCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGAGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCcatcat     >  8:117747/1‑148 (MQ=255)
                                                                                                                                               aaGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCAtcatca   <  7:31579/148‑1 (MQ=255)
                                                                                                                                               aaGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCAtcatca   <  1:283845/148‑1 (MQ=255)
                                                                                                                                               aaGCCCGACTATTTAAGCTTACTGCGCCTCTTCCGCTTCAATCTGCGTTACGTCCACCTCGGAATCGAAAATGATCAGCTGTTCCTCCGCAAGAAAGGACCGAGGCTGCATGTGATCCTGTTCGAACTGCCGCTCATGGCCATCATCAg  >  6:289599/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCGACGTCGAATGGGAATTCCGCTGCCGCAACGGTGAGGACCTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAG  >  NC_002947/5535090‑5535381

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: