Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,592:1 +C 100% coding (960/993 nt) 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,866,5881.C97.7% 307.8 / ‑0.0 88coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (48/38);  total (49/39)
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.92e-01

CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682
                                                                                                    |                                                                                                
cATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                   <  2:948406/100‑1 (MQ=255)
        cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccctgc                                                                                         <  1:1029408/100‑5 (MQ=255)
         ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccctgcg                                                                                        >  1:216031/1‑95 (MQ=255)
              ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGc                                                                                   >  1:1483666/1‑100 (MQ=255)
               tttgCACATAAAAGTCCAGCCTCGCTTGAGATTGGCCTTGCTGACAGATGCTTACTTACTCGCGGAACTGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCa                                                                                  <  2:1900471/98‑1 (MQ=255)
               tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCa                                                                                  <  2:1094876/100‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAg                                                                                 <  1:1116269/100‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAg                                                                                 <  1:2190941/100‑1 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAgg                                                                                >  1:1735230/1‑100 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAgg                                                                                <  2:779759/100‑1 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGa                                                                               >  2:722618/1‑100 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGa                                                                               >  2:1289250/1‑100 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGATTCGATATTCAGCC‑CCC‑TGCGTTTGCAgg                                                                                >  2:652091/1‑99 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGAt                                                                              <  1:2028204/100‑1 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGAttt                                                                            >  1:749922/1‑100 (MQ=255)
                       gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTcg                                                                          >  1:1125598/1‑100 (MQ=255)
                       gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTcg                                                                          >  1:976322/1‑100 (MQ=255)
                       gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTcg                                                                          >  2:1793748/1‑100 (MQ=255)
                       gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTcg                                                                          >  1:1586872/1‑100 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTcgcg                                                                        <  2:1013835/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTcgcg                                                                        <  1:2136218/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑CCCATGCGTTTGCAGGATTTcgcg                                                                        <  2:259104/100‑1 (MQ=255)
                          aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCgcgc                                                                       <  1:283884/100‑1 (MQ=255)
                          aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCgcgc                                                                       <  1:1179904/100‑1 (MQ=255)
                           aGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCa                                                                      <  1:668405/100‑1 (MQ=255)
                             gCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCaaa                                                                    >  1:2223119/1‑100 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCaa                                                                     >  2:1449548/1‑98 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAAc                                                                   <  2:2174537/100‑1 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  2:1267874/1‑100 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:576920/1‑100 (MQ=255)
                                aGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCATGATTTCGCGCAAACgg                                                                 >  2:2081849/1‑100 (MQ=255)
                                aGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACgg                                                                 <  1:1457271/100‑1 (MQ=255)
                                   cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              <  1:1386252/100‑1 (MQ=255)
                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAgn                                                           <  1:1053382/100‑2 (MQ=255)
                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  1:1589405/100‑1 (MQ=255)
                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  1:926867/100‑1 (MQ=255)
                                       cTTGAGACTGGCTTTTCTGACAGATGCTTAATTACTCGCGGAACAGCGCTTCGATATTAAG‑CGCCC‑TGAGTTTGCAGAATTTCGCGAAAACGGCGCAg                                                            >  2:1515975/1‑98 (MQ=255)
                                        ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  1:483721/100‑1 (MQ=255)
                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATCTTCAGCC‑CCC‑TGCGTTTACAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  2:213572/1‑100 (MQ=255)
                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  2:1209370/1‑100 (MQ=255)
                                          gagaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       <  2:1103313/100‑1 (MQ=255)
                                           agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       <  2:1667733/99‑1 (MQ=255)
                                            gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  1:179591/1‑100 (MQ=255)
                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  2:508600/100‑1 (MQ=255)
                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:1010269/100‑1 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 <  1:1247915/100‑1 (MQ=255)
                                                 gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  1:339165/1‑100 (MQ=255)
                                                  ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  2:2226229/1‑100 (MQ=255)
                                                   cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              >  2:1420508/1‑100 (MQ=255)
                                                    tttCTGATAGATGCTTACTGACTCGCGGAACACCGCTTCGATATTCAGCC‑CCC‑TGCTTTTGCAGGATTTGGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              <  1:1159331/99‑1 (MQ=255)
                                                     ttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATc                                            >  1:1467490/1‑100 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  1:1944058/1‑100 (MQ=255)
                                                         gACAGATGCTTACGTACTCGCGGAACAGCGCTTCGATAATCAGCC‑CCC‑TGCGTTCGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  2:1380925/1‑100 (MQ=255)
                                                          aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCATGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  1:2081863/100‑1 (MQ=255)
                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  2:1327235/1‑100 (MQ=255)
                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    >  2:2119152/1‑100 (MQ=255)
                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  2:1738284/100‑1 (MQ=255)
                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    >  1:327855/1‑100 (MQ=255)
                                                                 cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  1:1798063/1‑100 (MQ=255)
                                                                   tacttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              >  1:1439627/1‑100 (MQ=255)
                                                                   tacttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              <  1:1833219/100‑1 (MQ=255)
                                                                     cttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  2:58871/1‑100 (MQ=255)
                                                                      ttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                           <  1:722618/100‑1 (MQ=255)
                                                                      ttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                           >  1:1233173/1‑100 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  2:1735219/100‑1 (MQ=255)
                                                                        actCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  1:2055179/100‑1 (MQ=255)
                                                                        actCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  2:666645/100‑1 (MQ=255)
                                                                        actCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  2:794632/100‑1 (MQ=255)
                                                                         ctCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:1566856/1‑100 (MQ=255)
                                                                          tCGCGGAACAGCGCTTCGATATTCAGCC‑CCC‑TGAGTCTGCAGGATTTCGCGCAACCGGCGCAGGCCCTCCACCTAAATCTGGCGACCACGTCCACGGGt                        >  2:655129/1‑99 (MQ=255)
                                                                           cgcgGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  1:1367347/1‑100 (MQ=255)
                                                                           cgcgGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  1:1128680/1‑100 (MQ=255)
                                                                           cgcgGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  1:2182215/1‑100 (MQ=255)
                                                                           cgcgGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:1110759/1‑100 (MQ=255)
                                                                             cgGAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  1:104351/100‑1 (MQ=255)
                                                                              ggAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   >  2:209040/1‑100 (MQ=255)
                                                                               gAACAGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGcc                  >  2:1737181/1‑100 (MQ=255)
                                                                                   aGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              >  1:2203024/1‑100 (MQ=255)
                                                                                   aGCGCTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              <  1:2226244/100‑1 (MQ=255)
                                                                                    gcgcTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  1:635662/1‑100 (MQ=255)
                                                                                     cgcTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTc            >  2:636066/1‑100 (MQ=255)
                                                                                      gcTTCGATATTCAGCC‑CCC‑TGAGTTTGCAGGATTTCGCGCAAACGGCGCAGGCATTCAACCTGAATCTGGCGCACACGTTCACGGGTGAGGACact               >  2:2023210/1‑94 (MQ=255)
                                                                                       cTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1663715/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:2023224/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  2:1012858/1‑100 (MQ=255)
                                                                                       cTTCGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGGTCACGGGTGAGGCCAATTTCAc          <  2:1300520/100‑1 (MQ=255)
                                                                                          cGATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  1:1938154/100‑1 (MQ=255)
                                                                                           gATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  2:1539937/1‑100 (MQ=255)
                                                                                           gATATTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGACTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAGTTCACGAcc      >  2:1113222/1‑100 (MQ=255)
                                                                                              atTCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc   >  2:697150/1‑100 (MQ=255)
                                                                                               ttCAGCC‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  1:1089873/1‑100 (MQ=255)
                                                                                                    |                                                                                                
CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑CCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682

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