Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,214,605 G→A G520E (GGG→GAG)  rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,214,6050GA98.3% 387.2 / ‑2.7 120G520E (GGG→GAG) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (60/58);  total (60/60)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 123 total aligned reads displayed.

CTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGAC  >  NC_000913/3214508‑3214701
                                                                                                 |                                                                                                
cTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTggagg                                                                                                <  2:68453/100‑1 (MQ=255)
 tGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTggagga                                                                                               >  1:937928/1‑100 (MQ=255)
   aTGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGAtt                                                                                             <  1:1119986/100‑1 (MQ=255)
    tGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGGTGATGAAGATTCGCATCTGGAGGAttt                                                                                            <  2:1759363/100‑1 (MQ=255)
    tGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGAttt                                                                                            <  2:702577/100‑1 (MQ=255)
    tGCCGGAAGACAAGATCCGCAAAGGGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGAtt                                                                                             <  1:1119549/99‑1 (MQ=255)
      ccGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCa                                                                                          >  1:1262259/1‑100 (MQ=255)
       cGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCAt                                                                                         <  2:2928752/100‑1 (MQ=255)
        ggAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATc                                                                                        <  2:1088882/100‑1 (MQ=255)
        ggAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATc                                                                                        >  2:1536774/1‑100 (MQ=255)
        ggAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATc                                                                                        <  1:1775811/100‑1 (MQ=255)
        ggAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGAGGAAGATTCGCATCTGGAGGATTTCATc                                                                                        >  2:1899937/1‑100 (MQ=255)
         gAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCg                                                                                       >  2:2207436/1‑100 (MQ=255)
         gAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCg                                                                                       <  2:1490159/100‑1 (MQ=255)
           aGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                                                                     <  2:7754/100‑1 (MQ=255)
           aGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAg                                                                                     >  1:2031227/1‑100 (MQ=255)
            gACAAGATCCGCAAAGTGCTTAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGTATTTCATCGAgg                                                                                    <  1:952414/100‑1 (MQ=255)
             aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGa                                                                                   >  1:1938160/1‑100 (MQ=255)
             aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGa                                                                                   >  2:2872463/1‑100 (MQ=255)
             aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGa                                                                                   >  1:455385/1‑100 (MQ=255)
              cAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGAt                                                                                  >  1:1965365/1‑100 (MQ=255)
              cAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGAt                                                                                  <  2:1971405/100‑1 (MQ=255)
               aaGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATa                                                                                 <  1:519155/100‑1 (MQ=255)
               aaGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATa                                                                                 <  1:3036815/100‑1 (MQ=255)
                aGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATac                                                                                <  1:1042967/100‑1 (MQ=255)
                aGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATCATGAAGATTCGCATCTGGAGGATTTCATCGAGGATac                                                                                <  2:1755105/100‑1 (MQ=255)
                 gATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATacc                                                                               <  2:2712328/100‑1 (MQ=255)
                 gATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATacc                                                                               >  2:2830197/1‑100 (MQ=255)
                  aTCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATacca                                                                              <  2:30405/100‑1 (MQ=255)
                   tCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATaccac                                                                             >  2:1965517/1‑100 (MQ=255)
                      gCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCt                                                                          >  1:3251253/1‑100 (MQ=255)
                        aaaGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCg                                                                        >  1:1643876/1‑100 (MQ=255)
                         aaGTGCTGAAGATCGCCAAAGATCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTc                                                                         <  2:2590036/98‑1 (MQ=255)
                          aGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAg                                                                      <  2:2129404/100‑1 (MQ=255)
                           gTGCTGAAGATCGCCAAAGGGCCAATCTCCATGGAAACGCCGATCGGTGATGATTAAGCTTCGCATCTTGAGGATTTCATCGAGGATACCACCCTCGAGc                                                                     <  2:1490886/100‑1 (MQ=255)
                           gTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGc                                                                     <  2:1203008/100‑1 (MQ=255)
                           gTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGc                                                                     <  1:599189/100‑1 (MQ=255)
                            tGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCt                                                                    >  2:2342021/1‑100 (MQ=255)
                              cTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGc                                                                  >  2:243657/1‑100 (MQ=255)
                               tGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGcc                                                                 <  2:1590766/100‑1 (MQ=255)
                               tGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGcc                                                                 <  2:3125083/100‑1 (MQ=255)
                                gAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCg                                                                >  2:320618/1‑100 (MQ=255)
                                gAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCg                                                                >  1:3148374/1‑100 (MQ=255)
                                 aaGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGc                                                               <  2:3214700/100‑1 (MQ=255)
                                  aGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCt                                                              >  1:514736/1‑100 (MQ=255)
                                  aGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCt                                                              <  2:831286/100‑1 (MQ=255)
                                   gATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTg                                                             >  1:1605612/1‑100 (MQ=255)
                                   gATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTg                                                             <  2:1305235/100‑1 (MQ=255)
                                    aTCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTgg                                                            <  1:2165252/100‑1 (MQ=255)
                                    aTCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTgg                                                            <  1:1456671/100‑1 (MQ=255)
                                     tCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGa                                                           >  1:2300599/1‑100 (MQ=255)
                                       gCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGAtt                                                         >  1:2630858/1‑100 (MQ=255)
                                        ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTTGAGCTGCCGCTGGATTc                                                        >  2:364943/1‑100 (MQ=255)
                                        ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGa                                                           >  2:1020219/1‑97 (MQ=255)
                                        ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTc                                                        <  1:333494/100‑1 (MQ=255)
                                           nagagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGc                                                     >  1:3018250/2‑100 (MQ=255)
                                           aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGc                                                     <  1:2583690/100‑1 (MQ=255)
                                            agagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCg                                                    >  2:1284054/1‑100 (MQ=255)
                                            agagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCg                                                    <  2:3013582/100‑1 (MQ=255)
                                               gCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGacc                                                 <  2:1570483/100‑1 (MQ=255)
                                                  aaTCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGaccacc                                              >  1:3074263/1‑100 (MQ=255)
                                                  aaTCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGaccacc                                              >  1:1970891/1‑100 (MQ=255)
                                                    tctcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACTCTCGAGCTGCCGCTCGATTCTGCGACCACCGa                                            >  2:2049933/1‑100 (MQ=255)
                                                       ccATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGTGTTGTCGAATGATACGGCGACCACCGAAAg                                         <  2:300932/100‑1 (MQ=255)
                                                       ccATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGATGGATTCGGCGACCACCGagag                                         <  2:97127/100‑1 (MQ=255)
                                                        cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGc                                        >  1:124201/1‑100 (MQ=255)
                                                        cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGc                                        <  1:50167/100‑1 (MQ=255)
                                                        cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCCGCCGCTGGATTCTGCGACCACCGAAAGc                                        >  2:2598917/1‑100 (MQ=255)
                                                           ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtg                                     <  2:1605595/100‑1 (MQ=255)
                                                            gAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgc                                    >  1:1299642/1‑100 (MQ=255)
                                                             aaaCTCCGAGCGGTGATGATGAAGATTCGCATGTGGAGGATTGCATCGAGGATACCACCCTCGTTTTGTTTAATGATACTGCGACCACCGAGAGCCtgcg                                   <  2:2522822/100‑1 (MQ=255)
                                                                cGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTTCCGCAGGATTCTGCGACCACCGAGAGCCtgcg                                   <  2:1073432/97‑1 (MQ=255)
                                                                cGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgc                                >  1:1821692/1‑100 (MQ=255)
                                                                cGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgc                                <  1:243663/100‑1 (MQ=255)
                                                                cGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgc                                >  1:873895/1‑100 (MQ=255)
                                                                 gCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgtg                               <  1:2207462/100‑3 (MQ=255)
                                                                   cGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATGCCACCCTCTTTTTGTTGAATGATACGGCGACCACCGAGAGCCTGCGTGCgg                              <  2:2978408/99‑1 (MQ=255)
                                                                     aTCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCaa                           <  2:937916/100‑1 (MQ=255)
                                                                      tCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAAc                          <  2:2031202/100‑1 (MQ=255)
                                                                        ggTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGc                        >  2:1831083/1‑100 (MQ=255)
                                                                         gTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTTGAGCTGCCGCAGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCa                       <  2:635836/100‑1 (MQ=255)
                                                                         gTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGTGCTGCCGCTGGATTCGGCGACCACCGAGAGCCtgcgtgcgttgccgca                       <  2:455377/100‑9 (MQ=255)
                                                                              gatgaAGATTCGCATCTGGAGGATTTCATCGAGGATACCAACCTCGAGCTGCCGCTGGATTCTTCGACCACCGAAAGCCTGCGCGCGGCCACGCccga                    >  2:1126071/1‑98 (MQ=255)
                                                                              gatgaAGATTCGCATCTGGAGGATTTCATCAAGAATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCCACGCacgacg                  >  2:2141938/1‑100 (MQ=255)
                                                                               atgaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGt                 >  1:2525852/1‑100 (MQ=255)
                                                                                tgaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTg                >  1:1114174/1‑100 (MQ=255)
                                                                                tgaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTg                >  1:848739/1‑100 (MQ=255)
                                                                                  aaGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgct              >  2:1924876/1‑100 (MQ=255)
                                                                                   aGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGTTGCCGCTGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTgctg             <  2:575297/100‑1 (MQ=255)
                                                                                   aGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctg             >  2:1924093/1‑100 (MQ=255)
                                                                                    gATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCTTGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctgg            <  2:2518719/100‑1 (MQ=255)
                                                                                    gATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctgg            <  1:2109784/100‑1 (MQ=255)
                                                                                      ttCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggct          >  1:328297/1‑100 (MQ=255)
                                                                                        cGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGctggctgg        <  1:2872498/100‑1 (MQ=255)
                                                                                        cGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGctggctgg        >  1:32920/1‑100 (MQ=255)
                                                                                          cATCTGGAGGATTTCATCGAGGATACCACCCTCGTGCTGCTTCATGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGcc      <  2:1940095/100‑1 (MQ=255)
                                                                                          cATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGcc      <  2:3265931/100‑1 (MQ=255)
                                                                                            tCTGGAGGATTTCATCGAGGATACCACCCTCGATCTGCCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTg    <  2:2365407/100‑1 (MQ=255)
                                                                                             cTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCATGCGTGCGGCAACGCACGACGTGCTGGATGGCCTGa   >  1:136944/1‑100 (MQ=255)
                                                                                              tGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGAc  >  1:1081235/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGAC  >  NC_000913/3214508‑3214701

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: