Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,205,438 T→C intergenic (‑328/‑287) purH ← / → rrsE phosphoribosylaminoimidazolecarboxamideformyltra nsferase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,205,4380TC100.0% 673.8 / NA 210intergenic (‑328/‑287)purH/rrsEphosphoribosylaminoimidazolecarboxamideformyltra nsferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (118/92);  total (118/92)

Only 100 of 222 total aligned reads displayed.

TAAAGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCA  >  NC_000913/4205348‑4205511
                                                                                          |                                                                         
  aaGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCAt                                                                       >  2:5782804/1‑93 (MQ=255)
  aaGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGACCCCAt                                                                       >  2:5869456/1‑93 (MQ=255)
   aGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATc                                                                      >  1:1214997/1‑93 (MQ=255)
   aGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATc                                                                      <  1:3578361/93‑1 (MQ=255)
   aGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATc                                                                      <  2:4465440/93‑1 (MQ=255)
   aGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATc                                                                      <  1:1963716/93‑1 (MQ=255)
   aGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATc                                                                      <  1:6344086/93‑1 (MQ=255)
    gCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCg                                                                     >  2:2964663/1‑93 (MQ=255)
    gCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCg                                                                     <  1:2507783/93‑1 (MQ=255)
      gggCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGac                                                                   >  2:5718573/1‑93 (MQ=255)
       ggCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGaca                                                                  <  1:5628638/93‑1 (MQ=255)
        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                 >  2:4989062/1‑93 (MQ=255)
        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                 >  2:4917018/1‑93 (MQ=255)
        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                 <  1:5399367/93‑1 (MQ=255)
        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                 <  1:5235982/93‑1 (MQ=255)
          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                               >  1:6257267/1‑93 (MQ=255)
          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                               >  2:719839/1‑93 (MQ=255)
          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                               >  2:1247412/1‑93 (MQ=255)
          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgc                                                               >  2:1181499/1‑92 (MQ=255)
            ttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcggcg                                                             >  2:3082726/1‑93 (MQ=255)
              gCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGa                                                           <  1:3926915/93‑1 (MQ=255)
              gCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGa                                                           <  2:6281632/93‑1 (MQ=255)
              gCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGa                                                           <  2:3825537/93‑1 (MQ=255)
               cTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAt                                                          >  1:3771261/1‑93 (MQ=255)
               cTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAt                                                          <  2:6464275/93‑1 (MQ=255)
                 ggggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTTCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgt                                                        >  2:4699142/1‑93 (MQ=255)
                 ggggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgt                                                        >  2:1016230/1‑93 (MQ=255)
                 ggggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgt                                                        >  1:3718342/1‑93 (MQ=255)
                 ggggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgg                                                        >  1:3127323/1‑92 (MQ=255)
                  gggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGGACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgtg                                                       >  1:6303784/1‑93 (MQ=255)
                  gggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGGtgtg                                                       >  1:4216823/1‑93 (MQ=255)
                  gggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgtg                                                       >  2:1050451/1‑93 (MQ=255)
                     aaTGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAAt                                                    <  1:878212/93‑1 (MQ=255)
                       tGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCa                                                  >  2:6549960/1‑93 (MQ=255)
                         cagcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACt                                                <  1:4855576/93‑1 (MQ=255)
                         cagcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACt                                                <  1:1824311/93‑1 (MQ=255)
                          agcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACtt                                               >  2:1953506/1‑93 (MQ=255)
                          agcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACtt                                               >  2:4403903/1‑93 (MQ=255)
                           gcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTc                                              >  1:4276375/1‑93 (MQ=255)
                             agtcagtcATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcac                                            >  1:1496623/1‑93 (MQ=255)
                               tcagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcacac                                          >  1:5700145/1‑93 (MQ=255)
                               tcagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcacac                                          <  1:608477/93‑1 (MQ=255)
                                 agtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaa                                        >  2:686582/1‑93 (MQ=255)
                                 agtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaa                                        >  1:1561760/1‑93 (MQ=255)
                                  gtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaaa                                       >  2:5747476/1‑93 (MQ=255)
                                   tcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAAc                                      >  2:4778244/1‑93 (MQ=255)
                                   tcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAAc                                      <  2:3302074/93‑1 (MQ=255)
                                   tcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACAAAAAc                                      >  1:4688059/1‑93 (MQ=255)
                                    caTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACa                                     <  2:5495914/93‑1 (MQ=255)
                                    caTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACa                                     >  2:5929434/1‑93 (MQ=255)
                                     aTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAg                                    >  1:3267447/1‑93 (MQ=255)
                                     aTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAg                                    >  2:1816058/1‑93 (MQ=255)
                                     aTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAg                                    >  1:4582078/1‑93 (MQ=255)
                                      tttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAGGGCGGATGTGGATCACTTCACACAAACAGc                                   >  2:2916528/1‑93 (MQ=255)
                                      tttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCGCACAAACAGc                                   >  1:298602/1‑93 (MQ=255)
                                      tttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGc                                   >  2:459580/1‑93 (MQ=255)
                                      tttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGc                                   >  1:4521559/1‑93 (MQ=255)
                                       ttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGcc                                  <  1:1690418/93‑1 (MQ=255)
                                         ttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcgg                                >  2:5224935/1‑93 (MQ=255)
                                          tCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggt                               >  2:5369532/1‑93 (MQ=255)
                                           cTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGGATGACTTCACACAAACAGCcggtt                              >  1:2827943/1‑93 (MQ=255)
                                           cTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACATAAACAGCcggtt                              <  1:302127/93‑1 (MQ=255)
                                           cTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggtt                              <  1:5984553/93‑1 (MQ=255)
                                            tGCAATTTGTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggctc                             >  1:6539055/1‑93 (MQ=255)
                                             gCAATTTTTCTATTGCGGCCTGCGGAGACCTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcg                            <  1:942076/93‑1 (MQ=255)
                                              cAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcgg                           >  2:3923625/1‑93 (MQ=255)
                                                aTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcggtt                         <  1:1836435/93‑1 (MQ=255)
                                                aTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcggtt                         <  2:3547465/93‑1 (MQ=255)
                                                aTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcggtt                         <  2:1169511/93‑1 (MQ=255)
                                                aTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggctcggtt                         >  2:511295/1‑93 (MQ=255)
                                                 tttttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTg                        <  1:4954003/93‑1 (MQ=255)
                                                  ttttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGa                       <  1:709394/93‑1 (MQ=255)
                                                   tttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGaa                      <  2:3368749/93‑1 (MQ=255)
                                                   tttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGaa                      >  2:3386547/1‑93 (MQ=255)
                                                   tttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGaa                      <  1:4431842/93‑1 (MQ=255)
                                                    ttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAag                     >  2:2624851/1‑93 (MQ=255)
                                                    ttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAag                     >  1:2482440/1‑93 (MQ=255)
                                                     tCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAaga                    <  1:4893134/93‑1 (MQ=255)
                                                      cTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAagag                   <  2:5593600/93‑1 (MQ=255)
                                                      cTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACATCACACAAACAGCCGGTTCGGTTGAagag                   <  1:5147620/93‑1 (MQ=255)
                                                        aTTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaa                 >  2:2271406/1‑93 (MQ=255)
                                                        aTTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaa                 >  1:2997398/1‑93 (MQ=255)
                                                        aTTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaa                 >  2:6313953/1‑93 (MQ=255)
                                                             ggCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATc            >  1:1132753/1‑93 (MQ=255)
                                                              gCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATcc           <  2:1452247/93‑1 (MQ=255)
                                                              gCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATcc           >  2:3095288/1‑93 (MQ=255)
                                                               ccTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCt          <  1:1381613/93‑1 (MQ=255)
                                                                 tGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGa        >  2:6244865/1‑93 (MQ=255)
                                                                  gCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGaa       <  1:100665/93‑1 (MQ=255)
                                                                   cGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGaaa      <  2:436613/93‑1 (MQ=37)
                                                                   cGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGaaa      <  1:4534847/93‑1 (MQ=37)
                                                                    ggAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAt     <  2:1512692/93‑1 (MQ=25)
                                                                    ggAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAt     >  1:2298931/1‑93 (MQ=25)
                                                                     gagaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAtt    >  1:5063057/1‑93 (MQ=21)
                                                                     gagaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAtt    >  1:343517/1‑93 (MQ=21)
                                                                     gagaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGCTGAAGAGAAAAATCCTGAAAtt    <  1:4096137/93‑1 (MQ=255)
                                                                       gaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCa  >  2:2422718/1‑93 (MQ=17)
                                                                       gaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCa  >  2:41233/1‑93 (MQ=17)
                                                                       gaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCa  >  2:2725329/1‑93 (MQ=17)
                                                                       gaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCa  >  1:686432/1‑93 (MQ=17)
                                                                                          |                                                                         
TAAAGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCA  >  NC_000913/4205348‑4205511

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