Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,535,315 G→A 6.6% R95H (CGC→CAC)  pncB → nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,535,3150GA 6.6% 173.5 / 2.2 60R95H (CGC→CAC) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base A (2/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCA  >  NC_002947/5535179‑5535453
                                                                                                                                         |                                                                                                                                          
cTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGt                                                                                                                               >  1:1095710/1‑150 (MQ=255)
cTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGt                                                                                                                               <  2:1095710/150‑1 (MQ=255)
      cccTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCCAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTg                                                                                                                         <  1:837131/150‑1 (MQ=255)
      cccTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTg                                                                                                                         >  1:642695/1‑150 (MQ=255)
      cccTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCACTTCAATCTGCGTTACGTg                                                                                                                         >  1:1126803/1‑150 (MQ=255)
         tACCTCGGCGAAATTCGCAATCAGCTCGAACAGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCAc                                                                                                                      >  2:1193188/1‑150 (MQ=255)
          aCCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACg                                                                                                                     >  2:1633959/1‑150 (MQ=255)
          aCCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACg                                                                                                                     >  2:1632409/1‑150 (MQ=255)
            cTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTc                                                                                                                   <  2:642695/150‑1 (MQ=255)
             tCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTc                                                                                                                                        >  2:1170472/1‑128 (MQ=255)
             tCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTc                                                                                                                                        <  1:1170472/128‑1 (MQ=255)
             tCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCg                                                                                                                  >  1:1146306/1‑150 (MQ=255)
              cGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCgg                                                                                                                 <  2:1446366/150‑1 (MQ=255)
               ggCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGc                                                                                                                <  1:1049573/150‑1 (MQ=255)
               ggCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGc                                                                                                                >  2:1071896/1‑150 (MQ=255)
                gCGAAATTCGCAATCAGCTCGCAAGGCTAGAAAGACTTCACCTGCTTCTTAGACAGCTGTGTGTCCTCTGACGCATCATCTTCCTCAAGCCCGACTCTTTACCCTTTCTGCGCCTCTCACGCTTAAATCTGCGTTACGGGCGCGCCGGCa                                                                                                               <  2:1391145/150‑1 (MQ=255)
                    aaTTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTTGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGa                                                                                                           <  1:1071896/150‑1 (MQ=255)
                            aTCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTGCGCTTCAAGCTGCGTTACGTGCACGTCGGCATCGAAAATGATc                                                                                                   <  1:1306838/150‑1 (MQ=255)
                            aTCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATc                                                                                                   >  1:1780672/1‑150 (MQ=255)
                             tCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCa                                                                                                  <  1:1633959/150‑1 (MQ=255)
                             tCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCa                                                                                                  <  1:1632409/150‑1 (MQ=255)
                             tCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCa                                                                                                  <  1:317380/150‑1 (MQ=255)
                                gCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCt                                                                                               >  2:1722876/1‑150 (MQ=255)
                                        ggCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTgc                                                                                       <  2:1146306/150‑1 (MQ=255)
                                             gCCGACCTCACCCTGGACGATGGCCAACTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTg                                                                                  >  2:838264/1‑150 (MQ=255)
                                                 gACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGaaa                                                                               >  1:937711/1‑150 (MQ=255)
                                                          cTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTgg                                                                      >  1:589257/1‑150 (MQ=255)
                                                           tGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGc                                                                     <  2:1238714/150‑1 (MQ=255)
                                                            ggACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCt                                                                    <  1:1173100/150‑1 (MQ=255)
                                                             gACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTg                                                                   >  1:1237859/1‑150 (MQ=255)
                                                                gATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCACTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCAt                                                                <  2:1126803/150‑1 (MQ=255)
                                                                  tGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTCCGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCAtgt                                                              >  2:278456/1‑150 (MQ=255)
                                                                     ccAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCACTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGAt                                                           >  2:1737657/1‑150 (MQ=255)
                                                                        gCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCt                                                        >  1:484750/1‑150 (MQ=255)
                                                                          tGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGt                                                      >  2:424314/1‑150 (MQ=255)
                                                                              cTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGa                                                  >  1:1765531/1‑150 (MQ=255)
                                                                               ttCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGaa                                                 <  1:1193188/150‑1 (MQ=255)
                                                                                tCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAg                                                <  1:838264/150‑1 (MQ=255)
                                                                                  cTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTg                                              >  1:1783367/1‑150 (MQ=255)
                                                                                   tCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGTGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGGAGTGc                                             >  2:130348/1‑150 (MQ=255)
                                                                                           aTCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGctgctg                                     <  1:424314/150‑1 (MQ=255)
                                                                                           aTCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGctgctg                                     >  2:775913/1‑150 (MQ=255)
                                                                                             cAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGc                                   >  2:83900/1‑150 (MQ=255)
                                                                                                 ttCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCcatc                               <  1:1722876/150‑1 (MQ=255)
                                                                                                 ttCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCcatc                               <  2:1783367/150‑1 (MQ=255)
                                                                                                   ccTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCcatcat                             <  2:1780672/150‑1 (MQ=255)
                                                                                                     tCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGctgctg                                     >  1:235454/1‑140 (MQ=255)
                                                                                                     tCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGctgctg                                     <  2:235454/140‑1 (MQ=255)
                                                                                                     tCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCAtcatca                           >  2:1929737/1‑150 (MQ=255)
                                                                                                     tCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCAtcatca                           >  2:1279147/1‑150 (MQ=255)
                                                                                                      cAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAg                          >  1:1245098/1‑150 (MQ=255)
                                                                                                       aaGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGc                         <  1:212521/150‑1 (MQ=255)
                                                                                                           ccGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAg                     >  2:1204478/1‑150 (MQ=255)
                                                                                                            cGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGc                         <  2:1690520/145‑1 (MQ=255)
                                                                                                            cGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGc                         >  1:1690520/1‑145 (MQ=255)
                                                                                                                  tttACGCTTTCTGCGCCTTTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCaa              <  1:217337/150‑1 (MQ=255)
                                                                                                                    tACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCGTCAGCGTAGTACTCGAcc            >  2:1080831/1‑150 (MQ=255)
                                                                                                                       gCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGATGGCCATCATCAGCGAAGTACGCAACCGCa         >  1:1459417/1‑150 (MQ=255)
                                                                                                                         tttCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAAt       >  1:1101574/1‑150 (MQ=255)
                                                                                                                              gcgcCTGTTCCACTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCa  <  1:1737657/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CTGCGCCCCTACCTCGGCGAAATTCGCAATCAGCTCGAACGGCTTGCCGACCTCACCCTGGACGATGGCCAGCTGGCCTTCCTCGAACGCATCAGCTTCCTCAAGCCCGACTTTTTACGCTTTCTGCGCCTGTTCCGCTTCAATCTGCGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCA  >  NC_002947/5535179‑5535453

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