Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,214,605 G→A 80.1% 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,6050GA80.1% 240.6 / 51.9 119G520E (GGG→GAG) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base G (10/14);  new base A (48/47);  total (58/61)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

GCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCG  >  NC_000913/3214507‑3214703
                                                                                                  |                                                                                                  
gCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTgggg                                                                                                   >  2:1384623/1‑100 (MQ=255)
gCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTgggg                                                                                                   <  1:1384621/100‑1 (MQ=255)
 cTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTggagg                                                                                                  <  2:712964/100‑1 (MQ=255)
  tGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTggagga                                                                                                 <  1:500593/100‑1 (MQ=255)
    aTGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGAtt                                                                                               <  1:1572318/100‑1 (MQ=255)
      gCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTc                                                                                             >  2:468383/1‑100 (MQ=255)
       ccGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCa                                                                                            >  2:164437/1‑100 (MQ=255)
       ccGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCa                                                                                            >  2:2268160/1‑100 (MQ=255)
       ccGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCa                                                                                            >  1:1898494/1‑100 (MQ=255)
       ccGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCa                                                                                            >  1:1137947/1‑100 (MQ=255)
        cGGAAGACAAGATCCGCAAAATGCTGGAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGGGATGATGAAGATTCGCATCTGGAGGATTTCAt                                                                                           >  1:1775216/1‑100 (MQ=255)
         ggAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATc                                                                                          <  1:866652/100‑1 (MQ=255)
          gAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCg                                                                                         <  2:1296525/100‑1 (MQ=255)
            aGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                                                                       <  1:1487017/100‑1 (MQ=255)
            aGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAg                                                                                       >  2:36958/1‑100 (MQ=255)
            aGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAg                                                                                       >  2:2385091/1‑100 (MQ=255)
               cAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGAt                                                                                    <  1:2454601/100‑1 (MQ=255)
               cAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGAt                                                                                    >  2:2555988/1‑100 (MQ=255)
               cAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGAt                                                                                    >  2:2454612/1‑100 (MQ=255)
                aaGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATa                                                                                   >  2:2164773/1‑100 (MQ=255)
                 aGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATac                                                                                  >  2:1491237/1‑100 (MQ=255)
                 aGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATac                                                                                  <  1:289809/100‑1 (MQ=255)
                 aGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATac                                                                                  >  2:1203380/1‑100 (MQ=255)
                  gATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATacc                                                                                 >  1:1654860/1‑100 (MQ=255)
                  gATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAACATTCGCATCTGGAGGATTTCATCGAGGATacc                                                                                 <  1:1190436/100‑1 (MQ=255)
                   aTCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATacca                                                                                <  1:1209992/100‑1 (MQ=255)
                     ccGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATaccacc                                                                              <  1:1330238/100‑1 (MQ=255)
                     ccGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATaccacc                                                                              >  2:1631829/1‑100 (MQ=255)
                       gCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATATGAAGGATTTCATAGAGGATACCACCCt                                                                            >  2:2464895/1‑100 (MQ=255)
                        cAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAGTCGCATCTGGAGGATTTCATCGAGGATACCACCCTc                                                                           <  2:1121118/100‑1 (MQ=255)
                         aaaGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCg                                                                          <  1:2727186/100‑1 (MQ=255)
                           aGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAAAGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGat                                                                        <  1:1991144/100‑2 (MQ=255)
                            gTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGc                                                                       <  1:2581026/100‑1 (MQ=255)
                            gTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGc                                                                       <  2:1394604/100‑1 (MQ=255)
                              gCTGAAGATCGCCAAAGAGCCAATCTCCATGGAGACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTg                                                                     >  1:326794/1‑100 (MQ=255)
                              gCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTg                                                                     <  2:1317797/100‑1 (MQ=255)
                               cTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCAAGCTGc                                                                    >  2:1190293/1‑100 (MQ=255)
                                tGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGcc                                                                   <  2:1919163/100‑1 (MQ=255)
                                   aGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCt                                                                <  2:1836302/100‑1 (MQ=255)
                                      tCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGa                                                             <  1:939695/100‑1 (MQ=255)
                                       cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGAt                                                            >  1:2514999/1‑100 (MQ=255)
                                         ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTc                                                          >  2:772009/1‑100 (MQ=255)
                                           aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTg                                                        <  1:1688050/100‑1 (MQ=255)
                                           aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTg                                                        <  2:92041/100‑1 (MQ=255)
                                            aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGc                                                       >  2:1324918/1‑100 (MQ=255)
                                            aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGc                                                       <  1:556904/100‑1 (MQ=255)
                                             agagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCg                                                      <  2:1292269/100‑1 (MQ=255)
                                                gCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGacc                                                   >  1:2418251/1‑100 (MQ=255)
                                                   aaTCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGaccacc                                                <  2:2782730/100‑1 (MQ=255)
                                                   aaTCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGaccacc                                                <  2:2134903/100‑1 (MQ=255)
                                                      ctcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaa                                             >  2:2528125/1‑100 (MQ=255)
                                                      ctcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaa                                             <  2:1083059/100‑1 (MQ=255)
                                                      ctcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaa                                             >  1:77655/1‑100 (MQ=255)
                                                        ccATGGAAACGCCGATCGGTGATGATGAAGATTCGCTTCTGTGGGATTTCATCGAGGATACCACACTCTTGTTGCTGCAGGATACGGCGACCACCGagag                                           <  1:947802/100‑1 (MQ=255)
                                                        ccATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAg                                           >  1:553349/1‑100 (MQ=255)
                                                         cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGc                                          >  1:1090532/1‑100 (MQ=255)
                                                          aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGcc                                         <  2:2379288/100‑1 (MQ=255)
                                                           tGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCt                                        >  2:2644849/1‑100 (MQ=255)
                                                            ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtg                                       >  1:1596551/1‑100 (MQ=255)
                                                              aaaCGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCggcg                                     >  1:72677/1‑100 (MQ=255)
                                                                 cGCCGATCGGTGATGATGAAGATTCGCATCTGGAGGACTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgc                                  >  2:510417/1‑100 (MQ=255)
                                                                  gCCTATCGGTGATGATGAAGATTCGCATCTGGAGGATGTCATCGAGGATACCACCCTCGTTTTGTCGCTGGATTCGGCGACCACCGAAAGCCtgcgtgcg                                 <  1:2668933/100‑1 (MQ=255)
                                                                    ctaTCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGATCTGCCGCAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGc                               <  1:1536635/98‑1 (MQ=255)
                                                                    cGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGTGCTGCTGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGc                               <  1:772007/100‑1 (MQ=255)
                                                                    cGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGc                               <  2:671974/100‑1 (MQ=255)
                                                                    cGATCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGc                               <  2:361531/100‑1 (MQ=255)
                                                                    cGAACGGTGATGATGAAGATTCGCATCTGGAGGATTTCAGCGAGGATACCACCCTCGAGCTGTCGCTGGATTCGGCGACCACCGAGAGCCTGCGTGCGGc                               <  1:2131746/100‑1 (MQ=255)
                                                                      aTCGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCaa                             >  2:878703/1‑100 (MQ=255)
                                                                        cGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACg                           >  2:1919873/1‑100 (MQ=255)
                                                                        cGGTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACg                           <  2:1819080/100‑1 (MQ=255)
                                                                         ggTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgcgttgccgc                          <  2:553350/100‑8 (MQ=255)
                                                                          gTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCa                         <  2:1081072/100‑1 (MQ=255)
                                                                          gTGATGATGAAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCa                         <  2:2750865/100‑1 (MQ=255)
                                                                           tgatgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCTTTCTTTTTAATGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCa                         <  1:1485857/99‑1 (MQ=255)
                                                                              tgatgaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCccgac                     >  1:205290/1‑100 (MQ=255)
                                                                              tgatgaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCacgac                     >  2:2769257/1‑100 (MQ=255)
                                                                              tgatgaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCCGCGGGCGGCAACGCacga                      >  1:421733/1‑99 (MQ=255)
                                                                              tgatgaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCAGCGACCACCGAAAGCCTGCGTGCGGCAACGCacgac                     >  2:2765495/1‑100 (MQ=255)
                                                                               gatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCacgacg                    >  1:1052024/1‑100 (MQ=255)
                                                                               gatgaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCacgacg                    >  2:663023/1‑100 (MQ=255)
                                                                                  gaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGTGCTTTCGCATGATTCTGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTgc                 <  1:2279657/100‑1 (MQ=255)
                                                                                  gaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgc                 >  2:27664/1‑100 (MQ=255)
                                                                                  gaAGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgc                 >  2:1724640/1‑100 (MQ=255)
                                                                                   aaGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgct                >  2:1175641/1‑100 (MQ=255)
                                                                                    aGATTCGGATCTGGAGTATTTCATCGAGGATACCACCCTCGATTTTTCTAAGGATTCGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTgctg               <  1:2251681/100‑1 (MQ=255)
                                                                                    aGATTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGGCTGCGTGCGGCAACGCACGACGTgct                >  1:1929777/1‑99 (MQ=255)
                                                                                     gATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctgg              <  2:2307226/100‑1 (MQ=255)
                                                                                      tttCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggc             <  1:825575/99‑1 (MQ=255)
                                                                                      aTTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCCACGCACGACgcgcg                >  1:1098936/1‑96 (MQ=255)
                                                                                      aTTCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggc             <  2:659957/100‑1 (MQ=255)
                                                                                       ttCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGTGCTGTCTATGGATACGGCGACCACCGAGAGCCTGCGTGCGGCAACGCACGACGTgctggct            <  1:61207/100‑1 (MQ=255)
                                                                                       ttCGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTgctggca            >  2:825577/1‑99 (MQ=255)
                                                                                         cGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGctggctgg          >  2:387663/1‑100 (MQ=255)
                                                                                         cGCATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGctggctgg          >  2:2057743/1‑100 (MQ=255)
                                                                                           cATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGcc        <  2:24934/100‑1 (MQ=255)
                                                                                           cATCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGcc        <  2:2225855/100‑1 (MQ=255)
                                                                                             tCTGGAGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTg      >  1:2431631/1‑100 (MQ=255)
                                                                                              cTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGa     <  2:530444/100‑1 (MQ=255)
                                                                                                 ggggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcg  >  2:2894463/1‑100 (MQ=255)
                                                                                                 ggggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACcg  >  2:254458/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
GCTGATGCCGGAAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGCAACGCACGACGTGCTGGCTGGCCTGACCG  >  NC_000913/3214507‑3214703

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: ‑