Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,535,848 T→A 6.1% S289T (TCG→ACG)  pncB → nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,535,8480TA 6.1% 192.1 / 2.7 66S289T (TCG→ACG) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (28/34);  new base A (2/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

ACAGCCAGATCGCTGCGCTGGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATC  >  NC_002947/5535702‑5535994
                                                                                                                                                  |                                                                                                                                                  
aCAGCCAGATCGCTGCGCTGGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCg                                                                                                                                                  <  7:294633/149‑1 (MQ=255)
aCAGCCAGATCGCTGCGCTGGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATACg                                                                                                                                                  >  3:129964/1‑149 (MQ=255)
aCAGCCAGATCGCTGCGCTGGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATACg                                                                                                                                                  >  3:124951/1‑149 (MQ=255)
aCAGCCAGATCGCTGCGCTGGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCCACTGCATCACCATGGATCCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGATGTTCGACGGCCTGCGCAATTATTCg                                                                                                                                                  >  8:299944/1‑149 (MQ=255)
 cAGCCAGATCGCTGCGCTGGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCgg                                                                                                                                                 <  3:241917/149‑1 (MQ=255)
   gCCAGATCGCTGCGCTGGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGc                                                                                                                                               <  8:42696/149‑1 (MQ=255)
      aGATCGCTGCGCTGGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAg                                                                                                                                            >  7:303272/1‑149 (MQ=255)
         tCGCTGCGCTGGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCg                                                                                                                                         <  2:96615/149‑1 (MQ=255)
                  tgGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCAATGATTCGGGCGAGCCGGTAGCCTgg                                                                                                                                >  6:140156/1‑149 (MQ=255)
                    gACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGc                                                                                                                              <  4:59116/149‑1 (MQ=255)
                    gACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGc                                                                                                                              >  2:169475/1‑149 (MQ=255)
                    gACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCATCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCATGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGc                                                                                                                              >  2:217629/1‑149 (MQ=255)
                              tACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGc                                                                                                                     >  3:289925/1‑148 (MQ=255)
                                        ccGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATACGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCAtt                                                                                                          <  4:124951/149‑1 (MQ=255)
                                        ccGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATACGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCAtt                                                                                                          <  4:129964/149‑1 (MQ=255)
                                         cgcgGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCa                                                                                                            <  4:172289/146‑1 (MQ=255)
                                         cgcgGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCAtag                                                                                                         >  3:172289/1‑147 (MQ=255)
                                           cgGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGcatgaagc                                                                                                        >  1:346047/1‑141 (MQ=255)
                                                   tCGGCATCGTCCTGACCGACTGCATCACATTGGATGCCTTACTCGGCGATTTAGAGCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATTATTCGGGCTAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACt                                                                                                <  7:299944/148‑1 (MQ=255)
                                                     ggCATCGCCGTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTGGACGTGTTCTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCGGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTgg                                                                                              <  3:163333/148‑1 (MQ=255)
                                                     ggCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTgg                                                                                              <  4:137270/148‑1 (MQ=255)
                                                     ggCATCGCACTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGACTGCGCCATGATTCGGGCGATCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTgg                                                                                              <  1:217629/148‑1 (MQ=255)
                                                      gCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTggg                                                                                             >  6:269349/1‑148 (MQ=255)
                                                             ccTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCa                                                                                                      <  2:299365/132‑1 (MQ=255)
                                                             ccTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCa                                                                                                      >  1:299365/1‑132 (MQ=255)
                                                              cTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGAccc                                                                                    >  2:8754/1‑149 (MQ=255)
                                                                   cGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGa                                                                               <  1:120027/149‑1 (MQ=255)
                                                                   cGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGa                                                                               >  2:120027/1‑149 (MQ=255)
                                                                    gACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGAc                                                                              >  5:61427/1‑149 (MQ=255)
                                                                     aCTGCATCACGATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGAcc                                                                             <  4:123150/149‑1 (MQ=255)
                                                                     aCTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCa                                                                                                      <  1:157520/124‑1 (MQ=255)
                                                                     aCTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCa                                                                                                      >  2:157520/1‑124 (MQ=255)
                                                                     aCTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGAcc                                                                             <  4:289925/149‑1 (MQ=255)
                                                                         cATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGa                                                                         >  7:215712/1‑149 (MQ=255)
                                                                             aCCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGc                                                                      <  1:8754/148‑1 (MQ=255)
                                                                               cATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTgg                                                                   <  6:49650/149‑1 (MQ=255)
                                                                                 tGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGgtg                                                                 <  6:45988/149‑1 (MQ=255)
                                                                                   gATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAg                                                                                                                       <  8:239338/93‑1 (MQ=255)
                                                                                   gATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAg                                                                                                                       >  7:239338/1‑93 (MQ=255)
                                                                                       ccTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTc                                                            <  2:165741/148‑1 (MQ=255)
                                                                                        cTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTcc                                                           >  2:180682/1‑148 (MQ=255)
                                                                                             tCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGc                                                     >  1:101248/1‑149 (MQ=255)
                                                                                               ggCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGc                                                     <  5:269349/147‑1 (MQ=255)
                                                                                               ggCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCt                                                   >  8:237871/1‑149 (MQ=255)
                                                                                                   aTTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAAc                                               <  1:169475/149‑1 (MQ=255)
                                                                                                    tttCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTACCGACGGCCTCAAcc                                              >  7:91464/1‑149 (MQ=255)
                                                                                                        gACCTGTACTTCGCCAAGCTTTTCGCCGGCCTTCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGAc                                          <  8:303272/149‑1 (MQ=255)
                                                                                                          ccgtcaCTTCTCAACGCTGCTCGCCCGCCGCCGCCAGGACTCGTGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGAccc                                        <  4:247002/144‑1 (MQ=255)
                                                                                                            tGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCcgc                                      >  4:292140/1‑149 (MQ=255)
                                                                                                            tGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCcgc                                      >  4:129897/1‑149 (MQ=255)
                                                                                                                      cAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGAc                                          >  3:70996/1‑135 (MQ=255)
                                                                                                                      cAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGAc                                          <  4:70996/135‑1 (MQ=255)
                                                                                                                       aaGCTGTTCGTCGGTCTTCGGGATGATTCGGGCGATCCGGTAGGCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGAGCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAAt                           <  8:91464/149‑1 (MQ=255)
                                                                                                                         gCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAAttt                         >  2:336517/1‑149 (MQ=255)
                                                                                                                             ttCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTcc                      <  1:336517/148‑1 (MQ=255)
                                                                                                                             ttCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTc                       <  8:215712/147‑1 (MQ=255)
                                                                                                                               cGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCcgc                    <  5:140156/148‑1 (MQ=255)
                                                                                                                                   ggCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCt               >  3:329522/1‑149 (MQ=255)
                                                                                                                                     ccTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgc             >  5:177568/1‑149 (MQ=255)
                                                                                                                                      cTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgcg            >  3:105585/1‑149 (MQ=255)
                                                                                                                                         cgcCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGcgcg          <  7:237871/148‑1 (MQ=255)
                                                                                                                                            cATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACg      >  4:88137/1‑149 (MQ=255)
                                                                                                                                             aTGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGc     <  5:144351/149‑1 (MQ=255)
                                                                                                                                                aTTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCAt   <  3:292140/148‑1 (MQ=255)
                                                                                                                                                aTTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCAt   <  6:61427/148‑1 (MQ=255)
                                                                                                                                                aTTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATc  <  4:329522/149‑1 (MQ=255)
                                                                                                                                                 ttCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgcg            <  8:64189/138‑1 (MQ=255)
                                                                                                                                                 ttCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCAGatcgga   >  7:64189/1‑142 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACAGCCAGATCGCTGCGCTGGACTGCTGGGTACGCGAGTACCGCGGCCTGCTCGGCATCGCCCTGACCGACTGCATCACCATGGATGCCTTCCTCGGCGATTTCGACCTGTACTTCGCCAAGCTGTTCGACGGCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATC  >  NC_002947/5535702‑5535994

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