Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,903 C→T 100% E217K (GAG→AAG)  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,9030CT98.8% 552.3 / ‑2.9 164E217K (GAG→AAG) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (74/87);  total (77/87)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 164 total aligned reads displayed.

TATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAA  >  NC_000913/2866807‑2866999
                                                                                                |                                                                                                
tATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATCTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGtt                                                                                               >  2:2101791/1‑100 (MQ=255)
 aTCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTa                                                                                              >  1:1830313/1‑100 (MQ=255)
  tCTTCCGGACCGTTCTCTTTTTTTTTTTTTTTTATGTCCAGCAACGCGTTTTCGGATTCACTACCCAGCGGGGTGTCTACCGAGGTTATTCGCTTGtt                                                                                               >  2:1959714/1‑98 (MQ=255)
   cTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTTTTTTTTTCTACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAg                                                                                            >  2:3504482/1‑100 (MQ=255)
    ttCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTACGa                                                                                           >  1:3130157/1‑100 (MQ=255)
     tCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGAc                                                                                          <  2:2310157/100‑1 (MQ=255)
        ggACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGaa                                                                                       <  1:1991726/100‑1 (MQ=255)
          aCCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGc                                                                                     <  1:3928142/100‑1 (MQ=255)
           ccGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCa                                                                                    >  1:3033091/1‑100 (MQ=255)
               tctcTTTTTCATCGGCCAGGAGGTCCAGCAACGCGTTTTCAGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATAc                                                                                 <  2:1938694/99‑1 (MQ=255)
                 tctTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGTCGAAGCATACGGc                                                                              <  2:2436910/100‑1 (MQ=255)
                 tctTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGc                                                                              >  1:3740904/1‑100 (MQ=255)
                   tttttCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTg                                                                            >  1:3669941/1‑100 (MQ=255)
                     tttCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGAc                                                                          <  2:2508209/100‑1 (MQ=255)
                       tCATCGGCCAGGATGTCCAGCAACGCTTTTTTTTTTTTTTTTCCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGGCGAAGCATACGGCTGActt                                                                        >  2:924567/1‑98 (MQ=255)
                        cATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTAGTTAAGACGAAGCATACGACTGACGtc                                                                       >  2:1684083/1‑100 (MQ=255)
                         aTCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGtca                                                                      <  1:1604077/100‑1 (MQ=255)
                           cGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGtcatc                                                                    <  1:2809546/100‑1 (MQ=255)
                            ggCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGtcatca                                                                   <  2:2817266/100‑1 (MQ=255)
                             gCCAGGGTGTGCAGCAACGCTTTTTCGGAATCAGCACCTAGCGGGGTGTCTACCGGGGTAATGCGATTGTTAAGACGAAGCATACGGCTGACGtcatc                                                                    <  2:325536/98‑1 (MQ=255)
                                aGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTg                                                               <  1:18164/100‑1 (MQ=255)
                                 ggATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTgg                                                              <  2:1762196/100‑1 (MQ=255)
                                  gATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGc                                                             <  1:2952003/100‑1 (MQ=255)
                                  gATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGc                                                             <  2:431787/100‑1 (MQ=255)
                                  gATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGc                                                             <  1:1318150/100‑1 (MQ=255)
                                   aTGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCt                                                            >  1:527848/1‑100 (MQ=255)
                                        cAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTGATCAACTGGCTTATc                                                        >  2:2302836/1‑99 (MQ=255)
                                        cAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTggataagcc                                                       <  1:2540849/100‑8 (MQ=255)
                                        cAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATcc                                                       <  1:3151622/100‑1 (MQ=255)
                                        cAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATcc                                                       <  2:2352029/100‑1 (MQ=255)
                                        cAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATcc                                                       >  1:1159588/1‑100 (MQ=255)
                                        cAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATcc                                                       <  1:1961204/100‑1 (MQ=255)
                                        cAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATcc                                                       <  1:2302848/100‑1 (MQ=255)
                                          gCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGATGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAg                                                     <  2:2111642/100‑1 (MQ=255)
                                          gCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCa                                                      <  2:665536/99‑1 (MQ=255)
                                          gCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAg                                                     >  1:2607963/1‑100 (MQ=255)
                                          gCAACGCTTTTTCGGAATCACCACCCAGCGGGATGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAg                                                     >  1:3743766/1‑100 (MQ=255)
                                           cAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGt                                                    >  2:3584526/1‑100 (MQ=255)
                                           cAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGt                                                    <  2:1509866/100‑1 (MQ=255)
                                            aaCGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGtt                                                   <  1:2477789/100‑1 (MQ=255)
                                            aaCGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGtt                                                   <  2:798067/100‑1 (MQ=255)
                                             aCGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTg                                                  >  1:1910369/1‑100 (MQ=255)
                                              cGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACAAAGCATACGGCTGACGTCATCAACTGGATAAGCCAGtt                                                   <  2:298530/98‑1 (MQ=255)
                                              cGCTTTGTCGGAATCACCACCCAGCGGGGTGTCTATCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGc                                                 <  1:2190022/100‑1 (MQ=255)
                                               gCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGATGTCATCAACTGGCTTATCCAGTTGct                                                >  1:2760044/1‑100 (MQ=255)
                                               gCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGATAAGCCAGTTgat                                                <  1:1690484/100‑3 (MQ=255)
                                               gCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGATAAGCCAGTTgat                                                <  2:2931680/100‑3 (MQ=255)
                                                cTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGATAAGCCAGTTgatg                                               <  2:1990395/100‑4 (MQ=255)
                                                 tttttCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGATAAGCCAGTTgatga                                              <  1:2927012/100‑5 (MQ=255)
                                                  ttttCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGATAAGCCAGTTgatgac                                             <  1:3102645/100‑6 (MQ=255)
                                                   tttCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGc                                            >  1:241464/1‑100 (MQ=255)
                                                      cGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGAt                                         <  2:63621/100‑1 (MQ=255)
                                                      cGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGAt                                         >  1:445744/1‑100 (MQ=255)
                                                      cGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTACTAAACTGGCTTATCCCGTTGCTCTGc                                            >  2:3244017/1‑97 (MQ=255)
                                                       tcggtagacaccCCGCTGGGTGTCTACCGAGGTAATGCACTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGAtc                                        >  2:2421847/9‑100 (MQ=255)
                                                       ggAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGAtc                                        <  1:1514033/100‑1 (MQ=255)
                                                       ggAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGAtc                                        <  1:445900/100‑1 (MQ=255)
                                                        gAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGAtct                                       <  1:1997498/100‑1 (MQ=255)
                                                        gAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGAtct                                       <  1:3871916/100‑1 (MQ=255)
                                                          aTCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATctct                                     <  1:1229719/100‑1 (MQ=255)
                                                             accaccCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTcc                                  >  1:687489/1‑100 (MQ=255)
                                                              ccaccCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCg                                 <  2:137665/100‑1 (MQ=255)
                                                               caccCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGc                                <  2:1927896/100‑1 (MQ=255)
                                                                accCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCa                               >  2:3749148/1‑100 (MQ=255)
                                                                accCAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCAGCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTcc                                  <  1:1439508/97‑1 (MQ=255)
                                                                 cccAGCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCAc                              <  1:358643/100‑1 (MQ=255)
                                                                     gCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTgg                          >  2:2686150/1‑100 (MQ=255)
                                                                     gCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTgg                          >  1:1409657/1‑100 (MQ=255)
                                                                     gCGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTgg                          >  1:3357312/1‑100 (MQ=255)
                                                                      cGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGt                         >  2:3446178/1‑100 (MQ=255)
                                                                      cGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGt                         <  2:933924/100‑1 (MQ=255)
                                                                      cGGGGTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGt                         <  2:1551178/100‑1 (MQ=255)
                                                                       ggggTGTCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCa                               <  2:506321/93‑1 (MQ=255)
                                                                          gtgtCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCAt                     <  1:385220/100‑1 (MQ=255)
                                                                           tgtCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATg                    <  1:2946446/100‑1 (MQ=255)
                                                                           tgtCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATg                    <  1:175457/100‑1 (MQ=255)
                                                                            gtCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCACCTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATgg                   >  2:2237328/1‑100 (MQ=255)
                                                                            gtCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATgg                   <  2:1419800/100‑1 (MQ=255)
                                                                             tCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGt                  >  2:3499383/1‑100 (MQ=255)
                                                                             tCTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGt                  <  1:3383370/100‑1 (MQ=255)
                                                                              cTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTc                 >  1:559754/1‑100 (MQ=255)
                                                                              cTACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCAGCACTTGGTTCATGGt                  >  1:883995/1‑99 (MQ=255)
                                                                               tACCGAGGTAATGCGCTTGTTGAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTcc                >  2:2889126/1‑100 (MQ=255)
                                                                               tACCGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTcc                <  2:925928/100‑1 (MQ=255)
                                                                               tACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTcc                >  2:103930/1‑100 (MQ=255)
                                                                                 ncGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAg              >  1:2077976/2‑100 (MQ=255)
                                                                                 ccGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAg              <  1:364887/100‑1 (MQ=255)
                                                                                  cGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGc             >  1:823212/1‑100 (MQ=255)
                                                                                  cGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGc             >  1:1225190/1‑100 (MQ=255)
                                                                                  cGAGGTAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGc             >  2:1591292/1‑100 (MQ=255)
                                                                                       tAATGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTGATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATg        <  2:1539281/100‑1 (MQ=255)
                                                                                         aTGCGCTTGTTAAGACGAAGCATACGGCTGCCGTCATCCCCTGGCTTATCCAGTTGCTTAGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATggg      <  2:2149917/100‑1 (MQ=255)
                                                                                         aTGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATggg      <  2:3669913/100‑1 (MQ=255)
                                                                                          tGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGa     >  1:471376/1‑100 (MQ=255)
                                                                                          tGCGCTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGa     >  1:1008894/1‑100 (MQ=255)
                                                                                           gcgcTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGAc    <  1:3036375/100‑1 (MQ=255)
                                                                                           gcgcTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGAc    >  1:1578774/1‑100 (MQ=255)
                                                                                            cgcTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACa   <  2:1983443/100‑1 (MQ=255)
                                                                                             gcTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACaa  >  2:452505/1‑100 (MQ=255)
                                                                                             gcTTGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACaa  >  1:3347701/1‑100 (MQ=255)
                                                                                                |                                                                                                
TATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCTACCGAGGTAATGCGCTCGTTAAGACGAAGCATACGGCTGACGTCATCAACTGGCTTATCCAGTTGCTCTGCGATCTCTTCCGCACTTGGTTCATGGTCCAGCTTATGGGACAA  >  NC_000913/2866807‑2866999

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