Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,214,415 A→C I457L (ATT→CTT)  rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,214,4150AC98.4% 195.3 / NA 64I457L (ATT→CTT) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (26/37);  minor base G (0/1);  total (26/38)
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.98e-01

GGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAG  >  NC_000913/3214277‑3214553
                                                                                                                                          |                                                                                                                                          
ggaggggtgaAAGCGGTTGATAAATTCGAATTCCGCCGGGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTga                                                                                                                                        <  2:852687/136‑1 (MQ=255)
 gTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATGACGGTGCATATGCTTgag                                                                                                                                       >  1:367115/1‑143 (MQ=255)
           atccgTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATcaaca                                                                                                                             <  1:1113798/139‑1 (MQ=255)
           aGCGGTTGATAAATTCGAATACCGCCGTGGTGACAAGTTCTCCACCGACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATcaaca                                                                                                                             <  2:1231372/143‑1 (MQ=255)
            gggcttGATAAATTCGTATACCGCCGTGGTTACAAGTTATCCACCTACGCAACGTGGTGGATCCGTCAGGCGATGAACCGCTCTATCGCGGATCAGGAGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATcaacaa                                                                                                                            <  1:43245/139‑1 (MQ=255)
            gCGGTTGATAAACTCGAATCCCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATcaacaa                                                                                                                            <  2:766991/143‑1 (MQ=255)
              ggTTCATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGc                                                                                                                          <  2:790501/143‑1 (MQ=255)
                    tAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAATCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCAACATCCGTCTTCCGGTGCCTATGCTTGAGACCAACAGCACGCTccccc                                                                                                                    >  2:1111435/1‑139 (MQ=255)
                     aaaTTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCg                                                                                                                   >  1:633190/1‑143 (MQ=255)
                     aaaTTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCg                                                                                                                   >  1:843941/1‑143 (MQ=255)
                             ctACCGCCGTGGTGACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTtctc                                                                                                           <  1:1396669/142‑1 (MQ=255)
                              tACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCg                                                                                                          <  1:233277/143‑1 (MQ=255)
                               actgccgTGATTACACGCGCTCCACCTCCCCAACCTGGTGGATCCGTCCGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGGGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGc                                                                                                         <  2:554312/140‑1 (MQ=255)
                                 cgccgTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCa                                                                                                       >  1:8405/1‑143 (MQ=255)
                                     gTGGTTACAAGTTCTCCACCTACGCAACCTGGTCGATCCGTCAGGCGATCACCCGCTCTATCGCGGAGCAGGCCCGCACACCCCGTCTTCCGCTGCCTATGCTTGAGAACATCAACAAGATGAACGGTATTTACCGCCAGacg                                                                                                   >  2:553019/1‑141 (MQ=255)
                                       ggTTACAACTCCTCCACCTACGCAACCTGGTGGCTCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGAtgct                                                                                                 <  1:487246/143‑1 (MQ=255)
                                        gTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCCGGTGCTTTTGATGACTTGCAGTGATTAATATCGTTGCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGAtgctg                                                                                                <  2:642974/143‑1 (MQ=255)
                                        gTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGCCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGAtgctg                                                                                                <  1:1210152/143‑1 (MQ=255)
                                         ttACAAGTTCTCCACCTACGCGACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGAtgctgc                                                                                               <  2:914554/143‑1 (MQ=255)
                                                ttCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGAt                                                                                        >  1:1305269/1‑143 (MQ=255)
                                                ttCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGAt                                                                                        >  1:794528/1‑143 (MQ=255)
                                                  ctcCACCTACGCAACCTGGTGGATTCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATgg                                                                                      <  2:9486/143‑1 (MQ=255)
                                                  ctcCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATgg                                                                                      <  2:164037/143‑1 (MQ=255)
                                                        cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTg                                                                                >  1:775345/1‑143 (MQ=255)
                                                          aCGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCAGTATTGCGGTGCATATGCTTGAGACCATCAACCAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAAATGGACCGTGaa                                                                              >  2:1158478/1‑143 (MQ=255)
                                                           cGCAACCTGGTGGCTCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAAc                                                                             <  1:274803/143‑1 (MQ=255)
                                                            tgactcTGGTGGATCCGTCAGGCGATCACCCGCTCGAGCGCGGATCGGGCGCGCACCTTCCGTATTCCGGTGCATATGCTTGAGACCATTAACAAGCTCAACCGCATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAAcc                                                                            <  2:377141/138‑1 (MQ=255)
                                                            gCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAGGCTCAACCGTATTTCTCGCCCGATGCTGCACGCGATGGGCCGTGCAcc                                                                            >  1:930618/1‑143 (MQ=255)
                                                              aaCCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGa                                                                          <  1:486333/143‑1 (MQ=255)
                                                                ccTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATGCGTATTCCGGTGCGTATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAAACgcag                                                                        >  2:129575/1‑140 (MQ=255)
                                                                 cTGTTGTATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGc                                                                       <  2:721455/143‑1 (MQ=255)
                                                                  tggtggATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGACCCGACggc                                                                      >  1:1116749/1‑141 (MQ=255)
                                                                  tggtggATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGcc                                                                      >  1:650751/1‑143 (MQ=255)
                                                                  tggtggATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGcc                                                                      <  1:1301453/143‑1 (MQ=255)
                                                                    gtggATCCGTCAGGCGGTCACCCGGTCTATCCCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCgg                                                                    <  2:353960/143‑1 (MQ=255)
                                                                     tggATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGAGAACCGACGCCGga                                                                   >  1:154948/1‑143 (MQ=255)
                                                                          ccGGCAGGCGATCACCCGCGCTAACGCGGATCAGGCGCGCACCATCCGTATTCCGGGGCATAAGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAAc                                                              <  2:1251400/143‑1 (MQ=255)
                                                                          ccGGCAGCCGATCTCCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAAc                                                              <  2:989389/143‑1 (MQ=255)
                                                                          ccGCCAGGCCACCACCCGCTCTCTCGCTGATCAGTCCCCTCCCATCCGTATTCCGGTCCATATGCTTGAGCCCCTCCACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAAc                                                              <  2:469555/143‑1 (MQ=255)
                                                                            gTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGCTGCTGCCCGCGCTGGGCCGTGCACCCACGCCGCAAGAACTg                                                            >  1:1227205/1‑143 (MQ=255)
                                                                             tCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTgg                                                           <  1:966839/143‑1 (MQ=255)
                                                                             ccAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTgg                                                           <  2:459790/142‑1 (MQ=255)
                                                                              cAGGCGATCACCAGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTCCATAGGCTTGAGACAATCAACGAGCGCAAGAGGAGTTATCACCAGATGCGGGAAGCGATGGGCAGAGACACGGCGCGGGAAGGACAggg                                                          >  2:1147384/1‑142 (MQ=255)
                                                                                       cctcGCTCTATCGCTGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGAACATCAACAAGCTCAACCGAATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTAt                                                 <  1:1147352/140‑1 (MQ=255)
                                                                                       aCCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCGACAAGCTCAACCGTATTTCTCGCCCGATGCTGCACGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGGTGAACGTAt                                                 >  1:785071/1‑143 (MQ=255)
                                                                                            ctctGTCGGGGATCAGGCGCGCACGATCCGTATTCCGGGGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGa                                            <  1:319968/143‑1 (MQ=255)
                                                                                             gctATCGCGGTTCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGAt                                           <  2:509942/142‑1 (MQ=255)
                                                                                               tctCGCGGATCAGGCGCGTACCGGCCGGCTTCCCGTGCATATGCTTGAGCCCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATCGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGc                                         <  2:349301/141‑1 (MQ=255)
                                                                                                 tCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCg                                       >  2:123181/1‑143 (MQ=255)
                                                                                                    cgGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCGCAACCGTATTCCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGaa                                    <  1:80027/143‑1 (MQ=255)
                                                                                                    ccgATCAGGCGCGCACCATCCGTATCCCGGTGCATATGCTTGAGACCAGCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAACAACTGGCTGAACGTATGCTGATGCCGGaa                                    <  1:80035/141‑1 (MQ=255)
                                                                                                      gATCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACACGCGCAACCGTATTTCTCGCCAGAGGCTGCAAGAGATGGGCCGCGAACCGACGCGGGAACAAACGACTGAACGCACGCAGCTGCCgggcga                                  >  2:467712/1‑139 (MQ=255)
                                                                                                        tCAGGCGCGCACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACa                                >  1:1176553/1‑143 (MQ=255)
                                                                                                            gcgcgcACCATCCGTATTCCGGTGCATATGCTTGAGGCCATCAACAAGCTCAACCGTATTTCTCGCCAGCTGCTGCAAGAGATGGGGCGTGAACCGACGCCGGAAGCACTGGCTGAACGCCTGCTGCGGCCGGACGACAAGac                            >  2:91123/1‑142 (MQ=255)
                                                                                                            gcgcgcACCATCCGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAATAACGAGCTGAACGGATGCCGATGCCGGAAGACAGGaa                            >  2:84824/1‑142 (MQ=255)
                                                                                                            gcgcgcACCATCCCTATTCCGGTGCATATGCCTGAGACCATCAACGAGCTCACCCGTATGTCTCGCCAGATGCTGCAAGATATGGGCCGTGAACCGACGCCGGAAGAAATGGCTGAACGTATGCTGATGCCGGAAGACAAGAt                            <  1:402148/143‑1 (MQ=255)
                                                                                                                       ccGTATTCCGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGc                 <  2:843941/143‑1 (MQ=255)
                                                                                                                           atgccGGTGCATATGGTTGAGACCATCAGCGAGCTGGACCGGAGTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGaa             <  2:393249/140‑1 (MQ=255)
                                                                                                                              ccGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGAt          >  1:990980/1‑143 (MQ=255)
                                                                                                                               cGGTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATc         <  2:1305269/143‑1 (MQ=255)
                                                                                                                               cGGTGCAGAGGCTGGAGACCATCAACAAGCTCAGCCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGGGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATc         <  2:1191465/143‑1 (MQ=255)
                                                                                                                                ggTGCATATGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCg        <  2:8405/143‑1 (MQ=255)
                                                                                                                                      tatGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCCAGAGATGGGCCGTGGACCGACGCCGGAAGAACCGGCTGACCGTATGCTCATGCGGGGAGACAAGAACCGCGAAGTGCCTATCATCGCCaaac  >  2:1399839/1‑142 (MQ=255)
                                                                                                                                      tatGCTTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACCAGATCCGCAGAGTGCCGAAGAGCGCCCAag  >  1:130164/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                          
GGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCATCCGTATTCCGGTGCATATGATTGAGACCATCAACAAGCTCAACCGTATTTCTCGCCAGATGCTGCAAGAGATGGGCCGTGAACCGACGCCGGAAGAACTGGCTGAACGTATGCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAG  >  NC_000913/3214277‑3214553

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