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.C96.7% 313.6 / 3.7 90coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (44/43);  total (45/45)
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

ATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866490‑2866682
                                                                                                   |                                                                                               
aTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c                                                                                                <  1:2335064/100‑1 (MQ=255)
aTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c                                                                                                <  2:900820/100‑1 (MQ=255)
aTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c                                                                                                <  1:776788/100‑1 (MQ=255)
 tCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c                                                                                              <  1:975684/100‑1 (MQ=255)
  cGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cc                                                                                             <  2:915545/100‑1 (MQ=255)
   ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑ccc                                                                                            >  2:933875/1‑100 (MQ=255)
    gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccc                                                                                           <  1:173091/100‑2 (MQ=255)
       cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccctgc                                                                                        >  1:272771/1‑96 (MQ=255)
       cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccctgc                                                                                        >  1:1140331/1‑96 (MQ=255)
        ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑cccctgcg                                                                                       >  1:1667669/1‑95 (MQ=255)
            aTTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTg                                                                                    <  2:2365330/100‑1 (MQ=255)
             ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGc                                                                                   <  2:978995/100‑1 (MQ=255)
             ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGc                                                                                   <  1:1810176/100‑1 (MQ=255)
              tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCa                                                                                  <  2:1631620/100‑1 (MQ=255)
              tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCa                                                                                  <  2:1326668/100‑1 (MQ=255)
               gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACATCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAg                                                                                 <  1:505262/100‑1 (MQ=255)
               gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTCTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAg                                                                                 >  2:1269253/1‑100 (MQ=255)
                 gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGa                                                                               >  2:1046866/1‑100 (MQ=255)
                 gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGa                                                                               >  2:2132334/1‑100 (MQ=255)
                  cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAt                                                                              >  2:866560/1‑100 (MQ=255)
                  cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAt                                                                              >  2:2106886/1‑100 (MQ=255)
                  cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAt                                                                              <  2:1425746/100‑1 (MQ=255)
                    caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAttt                                                                            >  1:1772915/1‑100 (MQ=255)
                    caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAttt                                                                            >  1:1897655/1‑100 (MQ=255)
                    caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGAttt                                                                            >  2:2305352/1‑100 (MQ=255)
                      gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTcg                                                                          >  1:1054783/1‑100 (MQ=255)
                      gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTcg                                                                          <  1:1877646/100‑1 (MQ=255)
                      gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACGCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGGTTGCAGGATTTcg                                                                          >  1:595412/1‑100 (MQ=255)
                       aaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTcgc                                                                         <  2:1111033/100‑1 (MQ=255)
                        aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTcgcg                                                                        <  1:1834787/100‑1 (MQ=255)
                             ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   <  2:15206/100‑1 (MQ=255)
                             ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  2:1648041/1‑100 (MQ=255)
                             ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  2:1380056/1‑100 (MQ=255)
                             ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCGTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  2:141019/1‑100 (MQ=255)
                                 ccTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGgcg                                                               <  1:1046870/100‑1 (MQ=255)
                                 ccTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGgcg                                                               >  1:2077007/1‑100 (MQ=255)
                                  cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑CCCTGCGTTTGCAGGATTTTGCGCAAACGGCGCa                                                             <  1:1597011/100‑1 (MQ=255)
                                   tCGCTTGAGACTGGCCTTTCTGACAGATTCTTACTTACTCGCCGAACAGCGCTTCGATAGTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             <  2:2357168/100‑1 (MQ=255)
                                   tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                             <  1:1764586/100‑1 (MQ=255)
                                     gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  1:241098/100‑1 (MQ=255)
                                      cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  2:136314/1‑100 (MQ=255)
                                        tGAGACTGGCCTTTCTGACAGATGCTTACTGACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        <  2:1927906/100‑1 (MQ=255)
                                         gagaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       >  1:1174897/1‑100 (MQ=255)
                                           gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  1:668854/1‑100 (MQ=255)
                                             cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   >  2:159931/1‑100 (MQ=255)
                                             cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  1:1535670/100‑1 (MQ=255)
                                              tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:866562/100‑1 (MQ=255)
                                               ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  2:1486620/1‑100 (MQ=255)
                                               ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 <  1:1164477/100‑1 (MQ=255)
                                               ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:2279731/1‑100 (MQ=255)
                                                gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGCTATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  2:2289505/1‑100 (MQ=255)
                                                 ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  2:2021153/1‑100 (MQ=255)
                                                  cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              <  1:1122570/100‑1 (MQ=255)
                                                   tttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAAt                                             <  1:1378433/100‑1 (MQ=255)
                                                   tttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAAt                                             >  1:10748/1‑100 (MQ=255)
                                                    ttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATc                                            <  1:707199/100‑1 (MQ=255)
                                                     tCTGACAGATTTTTACTTACTCGCTGAACAGCCCTTCGATATTTAG‑CGCCCTGCGTTTGCAGGATTTCGCGCAACCGGCCCAGGCCTTCAACCTGAATCt                                           <  2:2035092/100‑1 (MQ=255)
                                                      cTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTg                                          >  1:484796/1‑100 (MQ=255)
                                                      cTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCACGATTTCACGCAAACGGCGCGGGCCTTCAACCTGAATCTg                                          >  2:1676103/1‑100 (MQ=255)
                                                       tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGTAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  2:1678139/1‑100 (MQ=255)
                                                       tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         <  1:1261318/100‑1 (MQ=255)
                                                       tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  2:1564216/1‑100 (MQ=255)
                                                        gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  1:1878657/1‑100 (MQ=255)
                                                        gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  2:1809482/1‑100 (MQ=255)
                                                         aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  1:1776577/100‑1 (MQ=255)
                                                         aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  1:1628130/100‑1 (MQ=255)
                                                          cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  2:780401/1‑100 (MQ=255)
                                                            gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  1:136315/100‑1 (MQ=255)
                                                            gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    >  1:580901/1‑100 (MQ=255)
                                                             aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   <  1:423001/100‑1 (MQ=255)
                                                              tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAAAGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  <  2:1426148/100‑1 (MQ=255)
                                                               gCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg                                 <  1:2202490/100‑1 (MQ=255)
                                                                gttacttacTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGAAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  1:2376287/2‑100 (MQ=255)
                                                                cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  1:2080373/1‑100 (MQ=255)
                                                                 ttacttacTCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGtt                               >  2:634664/1‑100 (MQ=255)
                                                                  tacttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              >  1:563554/1‑100 (MQ=255)
                                                                  tacttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              <  1:1264823/100‑1 (MQ=255)
                                                                     ttactCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                           <  1:1465052/100‑1 (MQ=255)
                                                                      tactCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          >  1:1335325/1‑100 (MQ=255)
                                                                        ctCGCGGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  2:801868/99‑1 (MQ=255)
                                                                          cgcgGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCTCAGGCCTTCAACCTGAATCTTGCGAACACGTTCACGGGTGa                      >  1:1901709/1‑100 (MQ=255)
                                                                          cgcgGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       <  2:2234579/99‑1 (MQ=255)
                                                                           gcgGAACAGCGCTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                     <  2:1975450/100‑1 (MQ=255)
                                                                                   gcgcTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             <  1:1983832/100‑1 (MQ=255)
                                                                                   gcgcTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  1:401870/1‑100 (MQ=255)
                                                                                      cTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1963242/100‑1 (MQ=255)
                                                                                      cTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:2305373/100‑1 (MQ=255)
                                                                                      cTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  2:1963226/1‑100 (MQ=255)
                                                                                      cTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  1:616983/1‑100 (MQ=255)
                                                                                      cTTCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  1:617162/1‑100 (MQ=255)
                                                                                        tCGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        <  1:1380064/100‑1 (MQ=255)
                                                                                         cGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGTCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  1:681845/100‑1 (MQ=255)
                                                                                         cGATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  1:89775/100‑1 (MQ=255)
                                                                                          gATATTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      <  1:1190380/100‑1 (MQ=255)
                                                                                          gATATTCAG‑C‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTGACGGGGGAGGCCAATTTCACGAcc      >  2:365185/1‑99 (MQ=255)
                                                                                            tatTCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    >  2:458317/1‑100 (MQ=255)
                                                                                              ttCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc   >  2:1416763/1‑99 (MQ=255)
                                                                                              ttCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  <  1:1020999/100‑1 (MQ=255)
                                                                                              ttCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  2:1538487/1‑100 (MQ=255)
                                                                                              ttCAGCC‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  <  1:2280870/100‑1 (MQ=255)
                                                                                                   |                                                                                               
ATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑CCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866490‑2866682

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