Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 961,334 C→T D39D (GAC→GAT rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913961,3340CT98.4% 421.5 / ‑3.2 123D39D (GAC→GATrpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (75/46);  total (76/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

Only 100 of 126 total aligned reads displayed.

GAAGAGTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGA  >  NC_000913/961245‑961423
                                                                                         |                                                                                         
gAAGAGTCCTTAAAAGAAATCGAAACCCGCCCGGGATCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGGAGTACTGGTTGATGCt                                                                                        >  1:1156232/1‑93 (MQ=255)
 aaGAGTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTg                                                                                       <  1:3951961/93‑1 (MQ=255)
 aaGAGTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTg                                                                                       <  2:1628200/93‑1 (MQ=255)
    agTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTc                                                                                    <  2:1891374/93‑1 (MQ=255)
    agTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTAGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTc                                                                                    <  1:2248299/93‑1 (MQ=255)
    agTCCTTAAAAGAAATCGAAACCCGCCAGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTc                                                                                    <  1:2885914/93‑1 (MQ=255)
     gTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCt                                                                                   >  2:4315536/1‑93 (MQ=255)
     gTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCt                                                                                   >  2:486386/1‑93 (MQ=255)
        cTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGaa                                                                                <  1:629514/93‑1 (MQ=255)
        cTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGaa                                                                                <  2:1351678/93‑1 (MQ=255)
        cTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGaa                                                                                >  1:1624570/1‑93 (MQ=255)
         ttAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGaaa                                                                               >  2:2907631/1‑93 (MQ=255)
          tAAAAGAAATCGAACCCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAAt                                                                              <  2:4259743/93‑1 (MQ=255)
          tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAAt                                                                              >  1:289910/1‑93 (MQ=255)
          tAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAAt                                                                              <  1:3800136/93‑1 (MQ=255)
           aaaaGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATc                                                                             <  1:4646602/93‑1 (MQ=255)
           aaaaGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATc                                                                             >  1:3169761/1‑93 (MQ=255)
            aaaGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCt                                                                            >  2:2490652/1‑93 (MQ=255)
             aaGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTg                                                                           <  1:2428626/93‑1 (MQ=255)
             aaGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTg                                                                           >  1:3636733/1‑93 (MQ=255)
             aaGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTGGTTGCTGGTCGGAAATCTg                                                                           >  2:1067902/1‑93 (MQ=255)
              aGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGa                                                                          <  2:3728467/93‑1 (MQ=255)
               gAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAg                                                                         >  2:758493/1‑93 (MQ=255)
                 aaTCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTc                                                                       >  1:1221472/1‑93 (MQ=255)
                  aTCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTcc                                                                      <  1:455026/93‑1 (MQ=255)
                   tCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCg                                                                     >  1:537020/1‑93 (MQ=255)
                   tCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCg                                                                     >  2:4583467/1‑93 (MQ=255)
                   tCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCg                                                                     >  1:3145906/1‑93 (MQ=255)
                   tCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTCACTCCg                                                                     >  1:537021/1‑93 (MQ=255)
                   tCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTAGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCg                                                                     >  2:1722761/1‑93 (MQ=255)
                    cGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGc                                                                    >  1:1886148/1‑93 (MQ=255)
                     gAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGcc                                                                   >  1:2187192/1‑93 (MQ=255)
                      aaaCCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCa                                                                  >  2:654309/1‑93 (MQ=255)
                       aaCACGCCCGGGATCTATCGTTCGGGGCGTTGATGTTGCTATCGACAAAGACGAAGTACTGGGTGATGCTGGTCTGAAATCTGAGCCCGCCAt                                                                 >  2:852711/1‑93 (MQ=255)
                         cccgcccgGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATcc                                                               <  2:3722170/93‑1 (MQ=255)
                           cgcccgGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCg                                                             <  1:1399000/93‑1 (MQ=255)
                           cgcccgGGTTCTATCGTTCGTGGCGGTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCg                                                             <  1:1398999/93‑1 (MQ=255)
                             cccgGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGc                                                           >  1:4184504/1‑93 (MQ=255)
                               cgGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTg                                                         <  2:2728156/93‑1 (MQ=255)
                                gggTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCACGGCTGa                                                        >  2:1930791/1‑93 (MQ=255)
                                 ggTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAg                                                       >  2:2627480/1‑93 (MQ=255)
                                  gTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGc                                                      >  1:3301163/1‑93 (MQ=255)
                                   ttCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCa                                                     >  1:3791896/1‑93 (MQ=255)
                                     cTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGt                                                   <  2:4179500/93‑1 (MQ=255)
                                     cTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGt                                                   <  2:1061896/93‑1 (MQ=255)
                                     cTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGt                                                   <  2:1061802/93‑1 (MQ=255)
                                      tATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGtt                                                  >  1:2890832/1‑93 (MQ=255)
                                       aTCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTc                                                 >  1:1196524/1‑93 (MQ=255)
                                       aTCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTc                                                 <  2:2374585/93‑1 (MQ=255)
                                        tcgttcgtGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGGCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCa                                                >  2:1541640/1‑93 (MQ=255)
                                         cgttcgtGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCaa                                               >  2:4417965/1‑93 (MQ=255)
                                           ttcgtGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCaaaa                                             >  1:670057/1‑93 (MQ=255)
                                            tcgtGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCaaaaa                                            <  2:2782457/93‑1 (MQ=255)
                                            tcgtGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCaaaaa                                            <  1:3409625/93‑1 (MQ=255)
                                             cgtGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTTGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAAc                                           <  1:2517893/93‑1 (MQ=255)
                                              gtGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACg                                          >  2:4736952/1‑93 (MQ=255)
                                              gtGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACg                                          >  1:3842078/1‑93 (MQ=255)
                                                 gCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGccc                                       >  2:1320542/1‑93 (MQ=255)
                                                  cGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCa                                      <  2:2101059/93‑1 (MQ=255)
                                                   gttgttgttGCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAg                                     >  1:1438505/1‑93 (MQ=255)
                                                    ttgttgttgCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAgg                                    >  1:3277365/1‑93 (MQ=255)
                                                    ttgttgttgCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAgg                                    >  1:2235721/1‑93 (MQ=255)
                                                    ttgttgttgCTATCGACAAAGACGTAGCACTGGGTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGTCTGAGCAGTTCAAAAACGCCCaag                                    >  2:3339722/1‑91 (MQ=255)
                                                     tgttgttgCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAggg                                   >  1:1072377/1‑93 (MQ=255)
                                                     tgttgttgCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAggg                                   <  1:2132384/93‑1 (MQ=255)
                                                      gttgttgCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGc                                  >  2:2921071/1‑93 (MQ=255)
                                                      gttgttgCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGc                                  <  1:2998940/93‑1 (MQ=255)
                                                       ttgttgCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCg                                 >  1:2004208/1‑93 (MQ=255)
                                                          ttgCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGc                              <  2:3382148/93‑1 (MQ=255)
                                                           tgCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCt                             <  1:771049/93‑1 (MQ=255)
                                                           tgCTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCt                             <  2:3079602/93‑1 (MQ=255)
                                                             cTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTgg                           >  2:2852304/1‑93 (MQ=255)
                                                             cTATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTgg                           >  2:3654128/1‑93 (MQ=255)
                                                              tATCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGa                          >  1:1766934/1‑93 (MQ=255)
                                                              tATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGa                          >  1:3994951/1‑93 (MQ=255)
                                                               aTCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGaa                         >  2:1870138/1‑93 (MQ=255)
                                                               aTCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGaa                         >  2:1870137/1‑93 (MQ=255)
                                                                tCGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGaaa                        >  1:3174399/1‑93 (MQ=255)
                                                                 cGACAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAAt                       >  1:4726300/1‑93 (MQ=255)
                                                                   aCAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATcc                     >  2:2469123/1‑93 (MQ=255)
                                                                    cAAAGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCa                    >  2:855702/1‑93 (MQ=255)
                                                                     aaaGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCag                   <  2:1624570/93‑1 (MQ=255)
                                                                       aGACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCCGGAAATCCagga                 >  1:3749048/1‑92 (MQ=255)
                                                                        gACGTAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCaggta                <  1:4157135/93‑1 (MQ=255)
                                                                        gACGCAGTACTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCaggta                >  1:4148047/1‑93 (MQ=255)
                                                                           gtagtaCTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCaggtaggc             <  1:971092/93‑2 (MQ=255)
                                                                            tagtaCTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTg            >  2:1257300/1‑93 (MQ=255)
                                                                             agtaCTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGa           >  1:169605/1‑93 (MQ=255)
                                                                             agtaCTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGa           >  2:4026296/1‑93 (MQ=255)
                                                                              gtaCTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGAc          <  2:1281561/93‑1 (MQ=255)
                                                                              gtaCTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGAc          <  1:3513809/93‑1 (MQ=255)
                                                                                 cTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGaa       <  1:1694199/93‑1 (MQ=255)
                                                                                 cTGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGaa       >  2:346734/1‑93 (MQ=255)
                                                                                  tGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAg      >  2:4394159/1‑93 (MQ=255)
                                                                                  tGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAg      <  1:3902984/93‑1 (MQ=255)
                                                                                  tGGTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAg      >  2:4394158/1‑93 (MQ=255)
                                                                                   ggTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGt     >  2:3678811/1‑93 (MQ=255)
                                                                                    gTTGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGtt    >  1:2202082/1‑93 (MQ=255)
                                                                                      tGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGa  <  2:666018/93‑1 (MQ=255)
                                                                                      tGATGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGa  >  2:3124385/1‑93 (MQ=255)
                                                                                         |                                                                                         
GAAGAGTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGA  >  NC_000913/961245‑961423

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: