Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,867,269 T→A 13.4% I95F (ATC→TTC)  rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,2690TA13.4% 262.1 / 33.0 82I95F (ATC→TTC) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base T (31/40);  new base A (4/7);  total (35/47)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

GCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCA  >  NC_000913/2867173‑2867366
                                                                                                |                                                                                                 
gCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatca                                                                                                 >  1:66012/1‑99 (MQ=255)
gCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatca                                                                                                 <  2:1319233/99‑1 (MQ=255)
 ccccAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatcat                                                                                                >  2:40699/1‑99 (MQ=255)
  cccAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatcatc                                                                                               <  2:844458/99‑1 (MQ=255)
     aGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcgg                                                                                            >  1:1546992/1‑99 (MQ=255)
      ggTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCcggc                                                                                           >  2:503641/1‑99 (MQ=255)
        ttGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcggcgg                                                                                         >  2:1242520/1‑99 (MQ=255)
           cccTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCga                                                                                      <  2:2034547/99‑1 (MQ=255)
           cccTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCga                                                                                      <  2:2172331/99‑1 (MQ=255)
               ctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAgg                                                                                  <  2:2287465/99‑1 (MQ=255)
                ttCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGc                                                                                 <  2:1606342/99‑1 (MQ=255)
                  cGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGa                                                                               <  2:1544185/99‑1 (MQ=255)
                   gATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGAc                                                                              >  1:128119/1‑99 (MQ=255)
                       aGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACAtct                                                                          <  1:1289502/99‑1 (MQ=255)
                       aGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACAtct                                                                          <  2:1634350/99‑1 (MQ=255)
                          tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCa                                                                       <  1:2475822/99‑1 (MQ=255)
                          tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCa                                                                       <  1:2040233/99‑1 (MQ=255)
                          tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCa                                                                       >  1:2376139/1‑99 (MQ=255)
                          tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCa                                                                       <  1:1448713/99‑1 (MQ=255)
                          tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCa                                                                       <  1:1733474/99‑1 (MQ=255)
                           ccAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCAc                                                                      <  2:2579767/99‑1 (MQ=255)
                            cAGCAACGACAGACCACGATTGCCATACCGGCGGGCAATTTTTACCACCAGACGCAAGCTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCAcc                                                                     >  2:349881/1‑98 (MQ=255)
                             aGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGc                                                                    >  1:387429/1‑99 (MQ=255)
                               cAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAg                                                                  <  2:973613/99‑1 (MQ=255)
                               cAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAg                                                                  <  2:318840/99‑1 (MQ=255)
                               cAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAg                                                                  <  1:1945163/99‑1 (MQ=255)
                                 aCGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTg                                                                <  1:1681834/99‑1 (MQ=255)
                                 aCGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTg                                                                >  2:849272/1‑99 (MQ=255)
                                 aCGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTg                                                                >  1:2400002/1‑99 (MQ=255)
                                  cGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgc                                                               <  1:668986/99‑1 (MQ=255)
                                    ccAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcgc                                                             <  2:230204/99‑1 (MQ=255)
                                    ccAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcgc                                                             <  1:1912511/99‑1 (MQ=255)
                                      aGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGa                                                           <  1:363259/99‑1 (MQ=255)
                                       gACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc                                                          >  1:2269115/1‑99 (MQ=255)
                                         ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcgc                                                        >  2:1559592/1‑99 (MQ=255)
                                          cACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcgcg                                                       <  2:1223483/99‑1 (MQ=255)
                                             gATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaa                                                    >  1:1129117/1‑99 (MQ=255)
                                              aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaaa                                                   >  1:1416626/1‑99 (MQ=255)
                                              aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaaa                                                   >  1:1297318/1‑99 (MQ=255)
                                              aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaaa                                                   >  2:1896564/1‑99 (MQ=255)
                                                   cATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATaaa                                              >  2:1610612/1‑99 (MQ=255)
                                                     tAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAAct                                            <  1:1004937/99‑1 (MQ=255)
                                                     tAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAAct                                            <  1:1277153/99‑1 (MQ=255)
                                                      aaCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActt                                           <  1:1201421/99‑1 (MQ=255)
                                                      aaCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActt                                           <  1:1994932/99‑1 (MQ=255)
                                                       aCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtc                                          >  1:1965463/1‑99 (MQ=255)
                                                       aCGGCGAGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtc                                          <  1:724514/99‑1 (MQ=255)
                                                         ggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtctt                                        <  2:799106/99‑1 (MQ=255)
                                                          gcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttc                                       <  2:738531/99‑1 (MQ=255)
                                                          gcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttc                                       <  1:1557140/99‑1 (MQ=255)
                                                          gcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttc                                       <  1:335076/99‑1 (MQ=255)
                                                          gcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttc                                       <  2:1077535/99‑1 (MQ=255)
                                                           cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttct                                      <  1:663630/99‑1 (MQ=255)
                                                            gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGTGCGACGCGCAAAATAAActtcttctt                                     <  1:840863/99‑1 (MQ=255)
                                                            gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttctt                                     <  2:1661018/99‑1 (MQ=255)
                                                            gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttctt                                     <  1:2622652/99‑1 (MQ=255)
                                                            gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttctt                                     >  1:691018/1‑99 (MQ=255)
                                                            gggCAATTATTACTACCAGACGCTTGTTACTCTCGATGATCAGGCCGCGAGAGGCGACACCGCTACTTAGTCCGCGACGCGCTAAATAAActtcttctt                                     >  2:2288045/1‑99 (MQ=255)
                                                              gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCg                                   <  1:2364635/99‑1 (MQ=255)
                                                              gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCg                                   <  1:2314785/99‑1 (MQ=255)
                                                              gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCg                                   >  2:735150/1‑99 (MQ=255)
                                                                aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGc                                 <  1:1923766/99‑1 (MQ=255)
                                                                 aTTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGcc                                >  1:2186898/1‑99 (MQ=255)
                                                                  tttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCg                               >  1:1502606/1‑99 (MQ=255)
                                                                      tACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTaa                           <  2:1644430/99‑1 (MQ=255)
                                                                         caccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAg                        >  1:1538674/1‑99 (MQ=255)
                                                                          accaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgt                       >  2:2603298/1‑99 (MQ=255)
                                                                            caGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgtgg                     >  1:1494907/1‑99 (MQ=255)
                                                                            caGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgtgg                     <  2:2248974/99‑1 (MQ=255)
                                                                                cGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGaa                 <  1:2030914/99‑1 (MQ=255)
                                                                                cGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGaa                 <  2:1457310/99‑1 (MQ=255)
                                                                                cGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGaa                 <  1:920138/99‑1 (MQ=255)
                                                                                  cAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATa               >  2:2200884/1‑99 (MQ=255)
                                                                                  cAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATa               >  1:1082221/1‑99 (MQ=255)
                                                                                  cAAGTTACTCTAGATCATCCGGCGGCGCGCGGCGACATGTCCCCGCAGTGCGCGACGCGCAAAATAAACTCCTTCTTCCGCCGTTAACAGCGGTGAATa               >  2:661747/1‑99 (MQ=255)
                                                                                   aaGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATaa              <  1:1374037/99‑1 (MQ=255)
                                                                                    aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAAc             >  1:1160886/1‑99 (MQ=255)
                                                                                    aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAAc             >  2:1653331/1‑99 (MQ=255)
                                                                                    aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAAc             >  2:1943106/1‑99 (MQ=255)
                                                                                      ttACTCTTGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCa           >  2:439234/1‑99 (MQ=255)
                                                                                      ttACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCa           <  2:1780730/99‑1 (MQ=255)
                                                                                         ctctCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAtc        >  1:1141266/1‑99 (MQ=255)
                                                                                          tctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAtct       >  1:635088/1‑99 (MQ=255)
                                                                                          tctcGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAtct       <  1:503641/99‑1 (MQ=255)
                                                                                             cGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCAc    <  2:973427/99‑1 (MQ=255)
                                                                                              gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCAcc   >  2:623553/1‑99 (MQ=255)
                                                                                               atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCa  >  1:1628820/1‑99 (MQ=255)
                                                                                                |                                                                                                 
GCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCA  >  NC_000913/2867173‑2867366

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