Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 972,118 T→A 5.4% L136Q (CTG→CAG)  secD → protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947972,1180TA 5.4% 237.9 / 3.6 74L136Q (CTG→CAG) secDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base A (2/2);  total (32/42)
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.76e-01

GCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACG  >  NC_002947/971972‑972261
                                                                                                                                                  |                                                                                                                                               
gccagccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgct                                                                                                                                              >  1:902445/1‑150 (MQ=255)
gccagccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgct                                                                                                                                              >  2:739475/1‑150 (MQ=255)
    gccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGaa                                                                                                                                          <  1:268833/150‑1 (MQ=255)
      ggaGGATGTAGTGCGCAAGGCACTGGGAGATGCTTACGTCGTGGCCCTGAATCTGGCCCAGATTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGt                                                                                                                                        <  2:268393/148‑1 (MQ=255)
       aaGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtg                                                                                                                                       <  2:472876/150‑1 (MQ=255)
        aGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtgg                                                                                                                                      >  2:1374566/1‑150 (MQ=255)
             gTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGGACTTCCTGCTGGCAGcggacatg                                                                                                                                 <  2:62967/150‑1 (MQ=255)
                  gcgcAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACaa                                                                                                                            >  1:1060170/1‑150 (MQ=255)
                      aaGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGcc                                                                                                                        >  1:1497829/1‑150 (MQ=255)
                          cACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAAGGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCTGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCAATGa                                                                                                                    <  1:706616/150‑1 (MQ=255)
                           aCTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACaa                                                                                                                            <  2:226804/141‑1 (MQ=255)
                           aCTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACaa                                                                                                                            >  1:226804/1‑141 (MQ=255)
                           aCTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGGGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGAc                                                                                                                   >  2:1418807/1‑150 (MQ=255)
                            cTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGATTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGGTGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACt                                                                                                                  <  1:699144/150‑1 (MQ=255)
                               ggCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGcc                                                                                                               >  1:72668/1‑150 (MQ=255)
                                    tgatTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCAt                                                                                                          <  2:902445/150‑1 (MQ=255)
                                          cgtcgtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGt                                                                                                    <  2:150237/150‑1 (MQ=255)
                                           gtcgtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCACAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTc                                                                                                   >  2:1190422/1‑150 (MQ=255)
                                              gtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTAc                                                                                                <  1:318304/150‑1 (MQ=255)
                                              gtGGCCCTGAACCTGACCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgct                                                                                                                                              >  1:1506249/1‑104 (MQ=255)
                                              gtGGCCCTGAACCTGACCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgct                                                                                                                                              <  2:1506249/104‑1 (MQ=255)
                                              gtGGCCCTGAAACTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGCGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTAc                                                                                                <  1:1570573/150‑1 (MQ=255)
                                                 gCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggaca                                                                                                                                   <  1:939830/112‑1 (MQ=255)
                                                 gCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggaca                                                                                                                                   >  2:939830/1‑112 (MQ=255)
                                                  cccTGAACCTGGCCCATACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACGTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAg                                                                                            <  2:1456802/150‑1 (MQ=255)
                                                  cccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAg                                                                                                                           <  2:1024283/119‑1 (MQ=255)
                                                  cccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAg                                                                                                                           >  1:1024283/1‑119 (MQ=255)
                                                    cTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGACACTGGGCCTGGACCTGTCCGGTGGTGCGCACTTCCTGCTGGTAGTGGACATGGACAAGGCCAAGAGTGCCCGCACGCAAGTGCCCGAAGGc                                                                                          >  2:616680/1‑150 (MQ=255)
                                                      gAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGa                                                                                        >  1:17853/1‑150 (MQ=255)
                                                      gAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATCGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGa                                                                                        >  2:1121731/1‑150 (MQ=255)
                                                       aaCCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAg                                                                                       <  2:1497829/150‑1 (MQ=255)
                                                         ccTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGt                                                                                     <  1:963055/150‑1 (MQ=255)
                                                              cccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacat                                                                                                                                  >  1:5591/1‑100 (MQ=255)
                                                              cccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacat                                                                                                                                  <  2:5591/100‑1 (MQ=255)
                                                               ccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGAc                                                                               >  1:1475026/1‑150 (MQ=255)
                                                                   aCTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtg                                                                           <  2:1475026/150‑1 (MQ=255)
                                                                         ccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcgc                                                                     <  1:1418807/150‑1 (MQ=255)
                                                                          cGCAATGGCTGCGCAACCTGGGCGCAAGCAAGATGAAACTGGGCCTGGACCAGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATAGACAAGGCCATGACTGTCCGCATGAAAGTCTAGGAAGGCGAGGTAAAGACCTTGCTGCGCa                                                                    >  2:491674/1‑150 (MQ=255)
                                                                              aTGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCGAGACCTTGCTGCGCAAaga                                                                >  2:533560/1‑150 (MQ=255)
                                                                                ggCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgc                                                              <  1:739475/150‑1 (MQ=255)
                                                                                ggCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgc                                                              <  2:1614808/150‑1 (MQ=255)
                                                                                ggCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgc                                                              >  2:905353/1‑150 (MQ=255)
                                                                                      gccaCCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTcc                                                        <  1:491674/147‑1 (MQ=255)
                                                                                           cTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGc                                                                                                                         >  1:780901/1‑80 (MQ=255)
                                                                                           cTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGc                                                                                                                         <  2:780901/80‑1 (MQ=255)
                                                                                           cTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGa                                                                                                                    >  1:779217/1‑85 (MQ=255)
                                                                                           cTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGa                                                                                                                    <  2:779217/85‑1 (MQ=255)
                                                                                            tGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtgctgc                                                                       <  1:810741/129‑1 (MQ=255)
                                                                                            tGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtgctgc                                                                       >  2:810741/1‑129 (MQ=255)
                                                                                                 gcAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                      >  2:260311/1‑93 (MQ=255)
                                                                                                 gcAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                      <  1:260311/93‑1 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtgctgc                                                                       >  1:991705/1‑116 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtgctgc                                                                       <  2:991705/116‑1 (MQ=255)
                                                                                                          aTGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCAt                                                                                                                      <  2:1134424/68‑1 (MQ=255)
                                                                                                          aTGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCAt                                                                                                                      >  1:1134424/1‑68 (MQ=255)
                                                                                                                cTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCaaa                                                                  >  2:1043243/1‑114 (MQ=255)
                                                                                                                cTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCaaa                                                                  <  1:1043243/114‑1 (MQ=255)
                                                                                                                 tGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATg                             >  1:936215/1‑150 (MQ=255)
                                                                                                                  gggCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCGAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATgg                            <  1:533560/150‑1 (MQ=255)
                                                                                                                   ggCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATggc                           <  2:1271305/150‑1 (MQ=255)
                                                                                                                     ccTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTa                                                                                                 >  1:49472/1‑78 (MQ=255)
                                                                                                                     ccTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTa                                                                                                 <  2:49472/78‑1 (MQ=255)
                                                                                                                       tGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAg                                                                                       >  1:613046/1‑86 (MQ=255)
                                                                                                                       tGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAg                                                                                       <  2:613046/86‑1 (MQ=255)
                                                                                                                         gACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGcatc                     <  1:1374566/150‑1 (MQ=255)
                                                                                                                               tCCGGTGGTGTGCACTGCCTGCGGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGTCGCCGTCTACACCTCGCTGCGCAAAGGGCGCATTGGCGGCGTCAGCCTGCCTCAGCAGTATGGCGGCATCATGCTg               <  2:185872/150‑1 (MQ=255)
                                                                                                                                 cGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTggg             >  1:530721/1‑150 (MQ=255)
                                                                                                                                  ggtggtGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGc            <  2:936215/150‑1 (MQ=255)
                                                                                                                                      gtgtGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCg        <  1:1579506/150‑1 (MQ=255)
                                                                                                                                      gtgtGCACTTCCTGCTGGAAGTGGACATGGACAAGGCAATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCg        <  1:905353/150‑1 (MQ=255)
                                                                                                                                         tccACTTCCAGCAGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATg                             >  2:540351/3‑126 (MQ=255)
                                                                                                                                         tccACTTCCAGCAGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATg                             <  1:540351/124‑1 (MQ=255)
                                                                                                                                          gCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGTAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGa    <  1:616680/150‑1 (MQ=255)
                                                                                                                                            aCTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtg                                                                           >  1:1601780/1‑77 (MQ=255)
                                                                                                                                            aCTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtg                                                                           <  2:1601780/77‑1 (MQ=255)
                                                                                                                                            aCTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACg  >  1:614851/1‑150 (MQ=255)
                                                                                                                                             cTTCCAGCAGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtgct                                                                         <  2:1210098/78‑1 (MQ=255)
                                                                                                                                             cTTCCAGCAGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtgct                                                                         >  1:1210098/1‑78 (MQ=255)
                                                                                                                                                  tgctgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgcgc                                                            >  2:573238/1‑86 (MQ=255)
                                                                                                                                                  tgctgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgcgc                                                            <  1:573238/86‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACG  >  NC_002947/971972‑972261

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