Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,214,148 G→A G368S (GGT→AGT)  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,1480GA98.8% 302.4 / ‑4.5 85G368S (GGT→AGT) rpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (35/49);  total (36/49)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AATTGCGATGAACAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTG  >  NC_000913/3214003‑3214295
                                                                                                                                                 |                                                                                                                                                   
aaTTGCGATGAACAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGa                                                                                                                                                 >  2:739417/1‑150 (MQ=255)
  ttGCGATGAACAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGc                                                                                                                                              >  1:551765/1‑151 (MQ=255)
   tGCGATGAACAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGc                                                                                                                                              <  1:776434/150‑1 (MQ=255)
   tGCGATGAACAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCg                                                                                                                                             >  1:1066839/1‑151 (MQ=255)
     cGATGAACAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGa                                                                                                                                            >  2:933517/1‑150 (MQ=255)
        tGAACAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGc                                                                                                                                        >  1:876438/1‑151 (MQ=255)
          aaCAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGccc                                                                                                                                      >  2:1143034/1‑151 (MQ=255)
          aaCAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGccc                                                                                                                                      <  2:23975/151‑1 (MQ=255)
            cAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCccgc                                                                                                                                    <  2:706873/151‑1 (MQ=255)
            cAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCccgc                                                                                                                                    <  1:1146196/151‑1 (MQ=255)
             aaGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCccgcc                                                                                                                                   <  1:86220/151‑1 (MQ=255)
             aaGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCccgcc                                                                                                                                   <  1:818347/151‑1 (MQ=255)
               gCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGt                                                                                                                                 >  2:889095/1‑151 (MQ=255)
                  gTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCg                                                                                                                              <  2:555636/151‑1 (MQ=255)
                  gTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCg                                                                                                                              >  1:206735/1‑151 (MQ=255)
                  gTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCg                                                                                                                              >  1:1104570/1‑151 (MQ=255)
                  gTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCg                                                                                                                              <  2:89655/151‑1 (MQ=255)
                          aaaaaCTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCgaagaa                                                                                                                         >  1:480261/1‑148 (MQ=255)
                          aaaaaCTGCACGATGTCTCTGAAAAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGAGCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGaagaaaga                                                                                                                      <  2:752743/151‑1 (MQ=255)
                           aaaaCTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATAACAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAagag                                                                                                                     <  1:557404/151‑1 (MQ=255)
                             aaCTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGaagaaag                                                                                                                       >  2:446611/1‑147 (MQ=255)
                                  cACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGtt                                                                                                               <  1:773886/150‑1 (MQ=255)
                                   aCGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTg                                                                                                              <  1:250894/150‑1 (MQ=255)
                                   aCGATGTCTCTGAAGAAGTGCATCGCGCCCGGCAAAAACTGCAGCAGATTGAAGAAGAACCCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGtt                                                                                                               <  1:860829/149‑1 (MQ=255)
                                    cGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGaa                                                                                                            <  2:818269/151‑1 (MQ=255)
                                         tctcTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGa                                                                                                        <  2:301006/150‑1 (MQ=255)
                                           tctGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACt                                                                                                     >  1:282512/1‑151 (MQ=255)
                                                agaagTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGt                                                                                                >  1:623698/1‑151 (MQ=255)
                                                agaagTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGt                                                                                                >  2:885401/1‑151 (MQ=255)
                                                   agTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTg                                                                                             <  1:695620/151‑1 (MQ=255)
                                                       cATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGtt                                                                                          >  2:913624/1‑150 (MQ=255)
                                                       cATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTa                                                                                         <  1:149988/151‑1 (MQ=255)
                                                          cgcgCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTAttt                                                                                      <  1:212092/151‑1 (MQ=255)
                                                              cccTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTAt                                                                                  <  2:926850/151‑1 (MQ=255)
                                                              cccTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTAt                                                                                  >  2:920492/1‑151 (MQ=255)
                                                              cccTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCACCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTAt                                                                                  >  1:1067638/1‑151 (MQ=255)
                                                               ccTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATc                                                                                 >  2:652518/1‑151 (MQ=255)
                                                                 tGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGc                                                                               >  1:1188612/1‑151 (MQ=255)
                                                                  gCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCt                                                                              <  1:69130/151‑1 (MQ=255)
                                                                     aaaaCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAg                                                                           <  2:60016/151‑1 (MQ=255)
                                                                      aaaCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTa                                                                             >  2:505603/1‑148 (MQ=255)
                                                                      aaaCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGa                                                                          <  1:980340/151‑1 (MQ=255)
                                                                       aaCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGa                                                                          >  2:750497/1‑150 (MQ=255)
                                                                        aCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGaa                                                                         >  1:819052/1‑150 (MQ=255)
                                                                                    tGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCg                                                             <  1:240963/150‑1 (MQ=255)
                                                                                     gaagaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGt                                                            <  1:333797/150‑1 (MQ=255)
                                                                                     gaagaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGt                                                            <  1:728622/150‑1 (MQ=255)
                                                                                     gaagaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTg                                                           >  1:109983/1‑151 (MQ=255)
                                                                                      aagaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTgg                                                          <  2:879395/151‑1 (MQ=255)
                                                                                      aagaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTgg                                                          >  1:164260/1‑151 (MQ=255)
                                                                                      aagaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTgg                                                          <  2:873555/151‑1 (MQ=255)
                                                                                       agaataaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGc                                                         >  2:839060/1‑151 (MQ=255)
                                                                                           gaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCAACAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGc                                                     <  1:786120/151‑1 (MQ=255)
                                                                                              aCCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGt                                                  <  2:650497/151‑1 (MQ=255)
                                                                                               ccGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGtt                                                 <  1:750468/151‑1 (MQ=255)
                                                                                                cGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTc                                                <  1:1154335/151‑1 (MQ=255)
                                                                                                  gCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCt                                              >  1:428980/1‑151 (MQ=255)
                                                                                                      gACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTg                                            <  2:615280/149‑1 (MQ=255)
                                                                                                       aCCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGAcc                                         >  2:1022235/1‑151 (MQ=255)
                                                                                                       aCCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGAcc                                         >  2:601658/1‑151 (MQ=255)
                                                                                                         cATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTg                                       <  2:85773/151‑1 (MQ=255)
                                                                                                           tCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGAt                                     >  2:854421/1‑151 (MQ=255)
                                                                                                                cAGGTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTc                                   <  2:746125/148‑1 (MQ=255)
                                                                                                                  ggTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGaa                              <  1:562532/151‑1 (MQ=255)
                                                                                                                  ggTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGa                               >  1:1156558/1‑150 (MQ=255)
                                                                                                                   gTTAAAGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGa                               >  2:880980/1‑149 (MQ=255)
                                                                                                                       aaGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCaa                         <  2:382033/151‑1 (MQ=255)
                                                                                                                        aGATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAAc                        >  1:53465/1‑151 (MQ=255)
                                                                                                                         gATATCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACa                       <  2:465610/151‑1 (MQ=255)
                                                                                                                            atCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCg                    <  1:161383/151‑1 (MQ=255)
                                                                                                                            atCAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCg                    <  1:311738/151‑1 (MQ=255)
                                                                                                                              cAACCGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGATTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCgg                   <  2:596199/150‑1 (MQ=255)
                                                                                                                                 ccGTCGTATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCt                <  1:57296/150‑1 (MQ=255)
                                                                                                                                      gtATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaa          <  1:162396/151‑1 (MQ=255)
                                                                                                                                       tATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaa          >  1:658345/1‑150 (MQ=255)
                                                                                                                                       tATGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaa          >  1:28102/1‑150 (MQ=255)
                                                                                                                                        aTGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa         <  2:114704/150‑1 (MQ=255)
                                                                                                                                        aTGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaa          >  2:770191/1‑149 (MQ=255)
                                                                                                                                        aTGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAg        <  1:920436/151‑1 (MQ=255)
                                                                                                                                         tGTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAg        <  1:612193/150‑1 (MQ=255)
                                                                                                                                          gTCCATCAGTTAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCTAACTTACGTCTGGTTATTTCTAGCGCTCAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCGACATCGGTCTGATGaaa         <  1:71976/148‑1 (MQ=255)
                                                                                                                                          gTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGc       <  2:32972/150‑1 (MQ=255)
                                                                                                                                          gTCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGc       <  1:430471/150‑1 (MQ=255)
                                                                                                                                           tCCATCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCgg     >  2:1004412/1‑151 (MQ=255)
                                                                                                                                              aTCAGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTg  <  1:108703/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AATTGCGATGAACAAGCCGTGGTCGGAAAAACTGCACGATGTCTCTGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTG  >  NC_000913/3214003‑3214295

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