Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 963,086 C→T 6.7% G364G (GGC→GGT rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913963,0860CT 6.7% 281.2 / 5.6 89G364G (GGC→GGTrpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base C (41/41);  new base T (3/3);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGA  >  NC_000913/962943‑963232
                                                                                                                                               |                                                                                                                                                  
cACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGAcc                                                                                                                                                >  2:29675/1‑148 (MQ=255)
 aCCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCg                                                                                                                                               >  8:278980/1‑148 (MQ=255)
  cccGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCgt                                                                                                                                              <  2:40107/148‑1 (MQ=255)
    cGTCCAAAGTTGTTAACGGTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTAGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGAcattgt                                                                                                                                            >  2:225235/1‑143 (MQ=255)
       ccAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTg                                                                                                                                          >  4:363897/1‑147 (MQ=255)
       ccAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTg                                                                                                                                          >  7:97092/1‑147 (MQ=255)
        cAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGaa                                                                                                                                        <  5:321496/148‑1 (MQ=255)
          aaGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTTTTGAAgg                                                                                                                                      >  7:366912/1‑148 (MQ=255)
              tgttAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCATAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTaa                                                                                                                                   >  4:309792/1‑147 (MQ=255)
                ttAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACCACAAGGGCGCCCGTGTTGAAGGTAAAAt                                                                                                                                >  6:56504/1‑148 (MQ=255)
                   aCGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGCCGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAATCTAACCCGTGGCAGCAGTTACCGGAAACCCACAACAAGCGTGACCGTGTTGAAGGTAAAATaaa                                                                                                                             >  5:122828/1‑148 (MQ=255)
                   aCGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCaa                                                                                                                             <  5:241617/148‑1 (MQ=255)
                   aCGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCa                                                                                                                              <  6:170709/147‑1 (MQ=255)
                     gTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGt                                                                                                                           <  7:203523/148‑1 (MQ=255)
                          cGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATc                                                                                                                      <  1:26080/148‑1 (MQ=255)
                                  tGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGGTGAAGGTAAAATCAAGTCTATCACTGACtt                                                                                                              >  6:67673/1‑148 (MQ=255)
                                   ggAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTc                                                                                                             >  6:299934/1‑148 (MQ=255)
                                        ttATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTAt                                                                                                        <  2:183572/148‑1 (MQ=255)
                                             gctCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTc                                                                                                    <  2:290192/145‑1 (MQ=255)
                                              ttCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCAt                                                                                                  <  8:366912/148‑1 (MQ=255)
                                                  ggATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGc                                                                                              >  6:230215/1‑148 (MQ=255)
                                                   gATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCt                                                                                             >  3:166855/1‑148 (MQ=255)
                                                     tatCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTg                                                                                           >  3:99064/1‑148 (MQ=255)
                                                       tCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTgg                                                                                          <  3:363897/147‑1 (MQ=255)
                                                            gaagaaCGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGAcggcg                                                                                    <  4:229771/148‑1 (MQ=255)
                                                              agaaCGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGAcggcgg                                                                                   >  1:363814/1‑147 (MQ=255)
                                                                 aCGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCAt                                                                                                  <  4:32916/129‑1 (MQ=255)
                                                                 aCGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCAt                                                                                                  >  3:32916/1‑129 (MQ=255)
                                                                   gtcgtcgtATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGa                                                                             <  4:166855/148‑1 (MQ=255)
                                                                     cgtcgtATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACg                                                                           >  7:87518/1‑148 (MQ=255)
                                                                        cgtATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGTGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGcc                                                                        >  3:251122/1‑148 (MQ=255)
                                                                         gtATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGcc                                                                        <  6:171000/147‑1 (MQ=255)
                                                                           aTCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACATGGGCGACCGTGTTGATGGTTAAATCACGTCTATCACTGAGTTCGGTATCTTCATCCGCTTGGAGGGCGGCATCGACGCCCTgg                                                                     >  1:338535/1‑148 (MQ=255)
                                                                           aTCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAg                                                                                                                            >  1:185918/1‑93 (MQ=255)
                                                                           aTCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAg                                                                                                                            <  2:185918/93‑1 (MQ=255)
                                                                             ctcCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAAt                                                                                                                                <  3:188793/87‑1 (MQ=255)
                                                                             ctcCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAAt                                                                                                                                >  4:188793/1‑87 (MQ=255)
                                                                             ctcCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTg                                                                                                                  <  3:244933/101‑1 (MQ=255)
                                                                             ctcCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTg                                                                                                                  >  4:244933/1‑101 (MQ=255)
                                                                             ctcCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGtt                                                                   >  7:251888/1‑148 (MQ=255)
                                                                               cccTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCa                                                                 >  2:330452/1‑148 (MQ=255)
                                                                               cccTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCa                                                                 >  2:269666/1‑148 (MQ=255)
                                                                                ccTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTg                                                                                                                  >  8:206580/1‑98 (MQ=255)
                                                                                ccTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTg                                                                                                                  <  7:206580/98‑1 (MQ=255)
                                                                                ccTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGtt                                                                   <  7:43877/145‑1 (MQ=255)
                                                                                 cTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCAcc                                                               <  8:97092/148‑1 (MQ=255)
                                                                                    ggTCTGAAACAGTGCAAAGCTAACCCGTTGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCtt                                                             <  1:225235/147‑2 (MQ=255)
                                                                                     gTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTc                                                           <  7:278980/148‑1 (MQ=255)
                                                                                      tCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACgg                                                                          >  6:183076/1‑132 (MQ=255)
                                                                                      tCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACgg                                                                          <  5:183076/132‑1 (MQ=255)
                                                                                         gAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGATCACCTGTCTGa                                                        >  1:14798/1‑147 (MQ=255)
                                                                                         gAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTGGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGa                                                        >  4:269154/1‑147 (MQ=255)
                                                                                          aaaCAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGa                                                        <  2:363814/146‑1 (MQ=255)
                                                                                          aaaCAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGAc                                                       <  1:29675/147‑1 (MQ=255)
                                                                                          aaaCAGTGCAAAGCTAACCCGTGGCAGCAGTTCGAGGAAACCCACAACATTGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGTTATCTTCCTCGGCTTGTACGGCGGCTTCGACGGCCTGGTTCACCTGTCTGAc                                                       <  4:167651/147‑1 (MQ=255)
                                                                                           aaCAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACACCAAGGGCGACCGTGTTGAAGCTAAAATCAAGTCTATCACTGACTTCGGAATCTTCATCGGCTAGGACGGCGGCATCGCCGGCCTGGTTCACCTGTCTGACAt                                                     >  5:14903/1‑148 (MQ=255)
                                                                                               gTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTcc                                                 <  4:99064/148‑1 (MQ=255)
                                                                                                 gCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTg                                               >  4:371035/1‑148 (MQ=255)
                                                                                                 gCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTg                                               >  8:40458/1‑148 (MQ=255)
                                                                                                     aGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGTGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCt                                                >  5:22745/1‑143 (MQ=255)
                                                                                                     aGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGTGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCt                                                <  6:22745/143‑1 (MQ=255)
                                                                                                     aGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAAc                                           >  2:284547/1‑148 (MQ=255)
                                                                                                       cTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGt                                         <  3:309792/148‑1 (MQ=255)
                                                                                                          aCCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGtt                                        <  1:30400/146‑1 (MQ=255)
                                                                                                          aCCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTg                                       <  5:299934/147‑1 (MQ=255)
                                                                                                              gTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGc                                  <  1:284547/148‑1 (MQ=255)
                                                                                                                 gcagcaGTTCGCGGAAACCCACAACAAGGGTGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCga                                <  6:122828/147‑1 (MQ=255)
                                                                                                                  cagcagTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCACGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACGTGTCTGACATCTCCTGGAACGTTGCAGGCgaag                              <  2:338535/148‑1 (MQ=255)
                                                                                                                      agTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGGCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGc                          <  1:269666/148‑1 (MQ=255)
                                                                                                                       gTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGaagaag                           >  8:287833/1‑146 (MQ=255)
                                                                                                                         tCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACt                                                                                                                   >  8:417846/1‑56 (MQ=255)
                                                                                                                         tCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACt                                                                                                                   <  7:417846/56‑1 (MQ=255)
                                                                                                                         tCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGt                       >  3:69736/1‑148 (MQ=255)
                                                                                                                         tCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGt                       >  2:193329/1‑148 (MQ=255)
                                                                                                                         tCGCGGAAACCCACAACAAGGGAGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAg                        >  6:269090/1‑147 (MQ=255)
                                                                                                                             ggAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCAGCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGt                   >  2:275581/1‑148 (MQ=255)
                                                                                                                              gAAACCCACAACAAGGGTGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTg                  <  4:251122/148‑1 (MQ=255)
                                                                                                                                aaCCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGaa                <  7:287833/148‑1 (MQ=255)
                                                                                                                                    cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATAc             >  2:40056/1‑147 (MQ=255)
                                                                                                                                    cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATAc             >  3:334078/1‑147 (MQ=255)
                                                                                                                                    cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCAGGAAAGTTTCAGGCGAAGAAGCAGTTCGTGAATAc             >  8:173601/1‑147 (MQ=255)
                                                                                                                                        acaaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaa        <  5:269090/148‑1 (MQ=255)
                                                                                                                                         caaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGa                 >  3:14671/1‑138 (MQ=255)
                                                                                                                                         caaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGa                 <  4:14671/138‑1 (MQ=255)
                                                                                                                                         caaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGGCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaa         >  3:276615/1‑146 (MQ=255)
                                                                                                                                          aaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGATCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAg      <  2:14798/148‑1 (MQ=255)
                                                                                                                                             ggCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGcg   <  7:40458/148‑1 (MQ=255)
                                                                                                                                             ggCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGcg   <  3:301226/148‑1 (MQ=255)
                                                                                                                                              gCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGcga  <  5:56504/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CACCCGTCCAAAGTTGTTAACGTTGGCGATGTAGTGGAAGTTATGGTTCTGGATATCGACGAAGAACGTCGTCGTATCTCCCTGGGTCTGAAACAGTGCAAAGCTAACCCGTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACTTCGGTATCTTCATCGGCTTGGACGGCGGCATCGACGGCCTGGTTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGA  >  NC_000913/962943‑963232

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