Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,866,592 +C 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.C98.3% 419.4 / 0.2 117coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (60/55);  total (61/56)
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.99e-01

Only 100 of 133 total aligned reads displayed.

CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682
                                                                                                    |                                                                                               
cATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                  >  2:207524/1‑100 (MQ=255)
 aTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c                                                                                                <  1:958824/100‑1 (MQ=255)
  tCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c                                                                                              <  1:174151/100‑1 (MQ=255)
    ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑ccc                                                                                            >  1:57213/1‑100 (MQ=255)
    ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑ccc                                                                                            >  1:2557066/1‑100 (MQ=255)
      aCCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑tccct                                                                                          <  1:2202264/100‑3 (MQ=255)
        cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccctgc                                                                                        >  1:804019/1‑96 (MQ=255)
        cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccctgc                                                                                        >  1:2461564/1‑96 (MQ=255)
         ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccctgcg                                                                                       >  2:235907/1‑95 (MQ=255)
         ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccctgcg                                                                                       >  2:617850/1‑95 (MQ=255)
         ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccctgcg                                                                                       >  1:1121927/1‑95 (MQ=255)
          tttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGt                                                                                       >  2:1119786/1‑100 (MQ=255)
          tttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGt                                                                                       <  2:2850344/100‑1 (MQ=255)
               tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCa                                                                                  >  2:156631/1‑100 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAg                                                                                 <  2:622741/100‑1 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAgg                                                                                >  1:2681919/1‑100 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGa                                                                               >  1:1340782/1‑100 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCGTACTTACTCGCGGAACAGCGCTTCGAGATTCAGCC‑CCCTGCGTTTGCAGGa                                                                               >  1:581384/1‑100 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAt                                                                              >  1:236262/1‑100 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAt                                                                              <  2:918686/100‑1 (MQ=255)
                   cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAt                                                                              >  2:2167288/1‑100 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAttt                                                                            >  1:189258/1‑100 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAttt                                                                            >  1:253404/1‑100 (MQ=255)
                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAttt                                                                            >  1:2379016/1‑100 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTcgcg                                                                        <  2:2807711/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTcgcg                                                                        <  2:874862/100‑1 (MQ=255)
                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTcgcg                                                                        <  2:876535/100‑1 (MQ=255)
                           aGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCa                                                                      >  1:1460948/1‑100 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  1:2283396/1‑100 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  1:1318620/1‑100 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  1:500670/1‑100 (MQ=255)
                             gCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCaaa                                                                    >  1:684633/1‑100 (MQ=255)
                               cAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACg                                                                  >  1:2781752/1‑100 (MQ=255)
                                 gCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGgc                                                                >  2:743888/1‑100 (MQ=255)
                                    tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             >  1:1121381/1‑100 (MQ=255)
                                    tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CTCCCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              >  1:1535873/1‑99 (MQ=255)
                                     cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            >  2:438950/1‑100 (MQ=255)
                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           >  2:2364761/1‑100 (MQ=255)
                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  1:94949/100‑1 (MQ=255)
                                       cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  1:1188520/1‑100 (MQ=255)
                                       cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  2:1032995/1‑100 (MQ=255)
                                        ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  2:459008/100‑1 (MQ=255)
                                        ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGCTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         >  1:184691/1‑100 (MQ=255)
                                          gagaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       <  2:2483557/100‑1 (MQ=255)
                                           agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      >  2:663535/1‑100 (MQ=255)
                                            ggcTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  2:2692928/3‑100 (MQ=255)
                                            gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     <  2:2251135/100‑1 (MQ=255)
                                            gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  2:1256965/1‑100 (MQ=255)
                                              cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  1:1119060/100‑1 (MQ=255)
                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  2:1866359/100‑1 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:2818032/1‑100 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACGTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:1686144/1‑100 (MQ=255)
                                                 gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcccg                                                <  1:438945/100‑3 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:2858846/100‑1 (MQ=255)
                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:551953/100‑1 (MQ=255)
                                                          aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  1:1032984/100‑1 (MQ=255)
                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      <  2:168831/100‑1 (MQ=255)
                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    >  1:673175/1‑100 (MQ=255)
                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  2:2509104/1‑100 (MQ=255)
                                                                  ttacttacTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGtt                               <  1:1166244/100‑1 (MQ=255)
                                                                   tacttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              <  1:1666343/100‑1 (MQ=255)
                                                                    acttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa                             >  2:452029/1‑100 (MQ=255)
                                                                     cttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  2:308512/1‑100 (MQ=255)
                                                                     cttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  1:2149461/1‑100 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          >  2:541816/1‑100 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  2:76191/100‑1 (MQ=255)
                                                                         ctCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  2:534368/1‑100 (MQ=255)
                                                                         ctCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        <  2:804025/100‑1 (MQ=255)
                                                                         ctCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:1456967/1‑100 (MQ=255)
                                                                          tCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       <  1:1415483/100‑1 (MQ=255)
                                                                          tCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       <  2:2379046/100‑1 (MQ=255)
                                                                           cgcgGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:2330824/1‑100 (MQ=255)
                                                                            gcgGAACAGCTCTTCGATATTCAGCC‑CCCTGCGTTTGTATTATTTTTTTTATACGGCGCAGTCCTTCAACCTTAATCTGGCGAACACGTTCACGGGTGAg                     <  2:1694334/100‑1 (MQ=255)
                                                                            gcgGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                     <  2:2485829/100‑1 (MQ=255)
                                                                             cgGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:2142193/100‑1 (MQ=255)
                                                                             cgGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:2437735/100‑1 (MQ=255)
                                                                             cgGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:1108713/100‑1 (MQ=255)
                                                                             cgGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:2550647/100‑1 (MQ=255)
                                                                             cgGAACAGCGCTTCGATATTCAG‑CGCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:2142430/100‑1 (MQ=255)
                                                                              ggAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   >  2:1399952/1‑100 (MQ=255)
                                                                              ggAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   <  2:1215391/100‑1 (MQ=255)
                                                                              ggAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   >  2:1510248/1‑100 (MQ=255)
                                                                               gAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGcc                  >  1:566824/1‑100 (MQ=255)
                                                                                aaCAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCa                 >  1:1268058/1‑100 (MQ=255)
                                                                                aaCAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCa                 <  1:2895722/100‑1 (MQ=255)
                                                                                 aCAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa                >  1:436811/1‑100 (MQ=255)
                                                                                 aCAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa                >  2:60247/1‑100 (MQ=255)
                                                                                  cAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt               >  2:1538379/1‑100 (MQ=255)
                                                                                    gcgcTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  1:437870/1‑100 (MQ=255)
                                                                                     cgcTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTc            <  1:541811/100‑1 (MQ=255)
                                                                                      gcTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:1686170/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:2364731/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:2121998/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  2:1981434/1‑100 (MQ=255)
                                                                                          cGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGCATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       >  1:28651/1‑100 (MQ=255)
                                                                                           gATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  2:2741606/1‑100 (MQ=255)
                                                                                             tatTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  1:2321749/100‑1 (MQ=255)
                                                                                             tatTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  1:2263256/100‑1 (MQ=255)
                                                                                             tatTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  1:525609/100‑1 (MQ=255)
                                                                                              cttCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc   <  1:1359264/99‑1 (MQ=255)
                                                                                                    |                                                                                               
CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682

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