Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,599,102 T→C H55R (CAT→CGT)  ydeW ← putative transcriptional regulator, sorC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,599,1020TC98.9% 610.4 / ‑3.5 188H55R (CAT→CGT) ydeWputative transcriptional regulator, sorC family
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (125/61);  total (127/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 193 total aligned reads displayed.

TGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAA  >  NC_000913/1599012‑1599191
                                                                                          |                                                                                          
tGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGc                                                                                         <  1:3648683/93‑1 (MQ=255)
  cAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGccc                                                                                       <  2:1749395/93‑1 (MQ=255)
    gCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCtt                                                                                     >  1:4584368/1‑93 (MQ=255)
      gAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtc                                                                                   >  2:3087181/1‑93 (MQ=255)
      gAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtc                                                                                   >  1:3001218/1‑93 (MQ=255)
       aaaaCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtct                                                                                  <  1:1266935/93‑1 (MQ=255)
        aaaCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtctc                                                                                 >  1:3312945/1‑93 (MQ=255)
         aaCTGGCGACGTAATTGAGTTTCATAATCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTcc                                                                                <  2:3744056/93‑1 (MQ=255)
         aaCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTcc                                                                                >  2:170504/1‑93 (MQ=255)
         aaCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTcc                                                                                <  2:3420424/93‑1 (MQ=255)
            tGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCCTTTCTCCAg                                                                              <  2:2383775/93‑1 (MQ=255)
            tGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGc                                                                             >  2:822357/1‑93 (MQ=255)
              acgacgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCaa                                                                           <  2:4162672/93‑1 (MQ=255)
                gacgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATc                                                                         >  2:4313655/1‑93 (MQ=255)
                 acgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATcg                                                                        >  2:4572153/1‑93 (MQ=255)
                  cgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATcgc                                                                       <  1:1470324/93‑1 (MQ=255)
                  cgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATcgc                                                                       <  2:2698304/93‑1 (MQ=255)
                  cgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATcgc                                                                       >  2:3613682/1‑93 (MQ=255)
                   gTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATcgcg                                                                      >  1:2076324/1‑93 (MQ=255)
                   gTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATcgcg                                                                      >  1:1217301/1‑93 (MQ=255)
                    tAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGa                                                                     >  1:594315/1‑93 (MQ=255)
                     aaTTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGac                                                                    <  1:4181076/93‑1 (MQ=255)
                      aTTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGaca                                                                   >  2:4046084/1‑93 (MQ=255)
                      aTTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGaca                                                                   >  2:2676365/1‑93 (MQ=255)
                       ttGAGTTTCATATTCCAGAGAGACTTCAAAGCGAGAATTAATCTGTACGCGAATAATGACGGACTGACGTCCTTTCTCCCGCAATCGCGacac                                                                  >  2:3211875/1‑93 (MQ=255)
                       ttGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGacac                                                                  >  1:1875383/1‑93 (MQ=255)
                             ttCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGCGACACTTTCaa                                                            >  1:2324386/1‑93 (MQ=255)
                             ttCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCaa                                                            >  2:4378012/1‑93 (MQ=255)
                             ttCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCATTTCTCCAGCAATCGCGACACTTTCaa                                                            >  2:2939971/1‑93 (MQ=255)
                              tCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCaaa                                                           >  1:2603967/1‑93 (MQ=255)
                              tCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCaaa                                                           >  1:127497/1‑93 (MQ=255)
                                atatTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTGTCTCCAGCAATCGCGACCCTTTCAAACg                                                         >  1:4486894/1‑93 (MQ=255)
                                 tatTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACgt                                                        >  2:2249386/1‑93 (MQ=255)
                                  atTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACgtg                                                       <  1:3437246/93‑1 (MQ=255)
                                  atTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACgtg                                                       >  1:1966745/1‑93 (MQ=255)
                                  atTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACgtg                                                       <  2:3942669/93‑1 (MQ=255)
                                    tCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTc                                                     <  1:4707213/93‑1 (MQ=255)
                                     ccAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGCGTCa                                                    >  2:739486/1‑93 (MQ=255)
                                      cAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAg                                                   <  2:1115863/93‑1 (MQ=255)
                                      cAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAg                                                   <  1:86904/93‑1 (MQ=255)
                                      cAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAg                                                   >  2:1416281/1‑93 (MQ=255)
                                        gACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGc                                                 >  2:4252756/1‑93 (MQ=255)
                                        gACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGc                                                 >  2:1175692/1‑93 (MQ=255)
                                         aCAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGcc                                                >  1:1404823/1‑93 (MQ=255)
                                         aCAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGcc                                                >  1:1791578/1‑93 (MQ=255)
                                         aCAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGcc                                                >  1:4422983/1‑93 (MQ=255)
                                         aCAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGcc                                                <  1:432336/93‑1 (MQ=255)
                                          cAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCg                                               >  1:574111/1‑93 (MQ=255)
                                            gCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCgag                                             >  1:1656834/1‑93 (MQ=255)
                                            gCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCgag                                             >  1:903669/1‑93 (MQ=255)
                                             gcTTCAAAGCGAGAATTTATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCgaga                                            <  1:1367706/92‑1 (MQ=255)
                                               ttCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACg                                          >  1:970343/1‑93 (MQ=255)
                                               ttCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACg                                          <  1:687721/93‑1 (MQ=255)
                                               ttCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACg                                          >  1:3873624/1‑93 (MQ=255)
                                               ttCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACg                                          >  1:2685792/1‑93 (MQ=255)
                                                tCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGa                                         >  1:3560172/1‑93 (MQ=255)
                                                tCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGa                                         >  2:3319099/1‑93 (MQ=255)
                                                tCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGa                                         >  2:550452/1‑93 (MQ=255)
                                                tCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGAAACTTTCAAACGTGTCAGGCCGAGACGa                                         <  1:3988512/93‑1 (MQ=255)
                                                 cAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGAt                                        <  1:613072/93‑1 (MQ=255)
                                                 cAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGAt                                        <  1:4517936/93‑1 (MQ=255)
                                                   aaGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGGCGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCg                                      <  1:2787026/93‑1 (MQ=255)
                                                   aaGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCg                                      >  1:4256587/1‑93 (MQ=255)
                                                    aGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGc                                     >  2:477790/1‑93 (MQ=255)
                                                     gCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCt                                    >  1:3190851/1‑93 (MQ=255)
                                                       gagaATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTATCAGGCCGAGACGATCGCTGa                                  >  2:4359261/1‑93 (MQ=255)
                                                        agaATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGAt                                 <  2:1279201/93‑1 (MQ=255)
                                                          aaTTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGAtct                               <  2:4444752/93‑1 (MQ=255)
                                                           aTTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATAGCTGAtctc                              >  2:3358621/1‑93 (MQ=255)
                                                             tAATCTGTACGCAAATAATGCCGGACTGACGCCCTTGCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGc                            <  2:2662528/93‑1 (MQ=255)
                                                               aTCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCGAGCAATAGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGctc                          >  1:2236019/1‑93 (MQ=255)
                                                               aTCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGCTCTCGctc                          >  2:3776172/1‑93 (MQ=255)
                                                                tCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGctct                         >  1:117489/1‑93 (MQ=255)
                                                                tCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGctct                         >  1:2224393/1‑93 (MQ=255)
                                                                tCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGctct                         <  2:1582854/93‑1 (MQ=255)
                                                                tCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGctct                         <  2:166310/93‑1 (MQ=255)
                                                                 cTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTg                        <  2:4641825/93‑1 (MQ=255)
                                                                  tGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTgg                       >  1:3140977/1‑93 (MQ=255)
                                                                   gTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTggg                      <  1:3829860/93‑1 (MQ=255)
                                                                    tACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTgggc                     <  1:103140/93‑2 (MQ=255)
                                                                    tACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGt                     <  2:2173972/93‑1 (MQ=255)
                                                                     aCGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTc                    <  1:2531875/93‑1 (MQ=255)
                                                                     aCGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTc                    <  2:556979/93‑1 (MQ=255)
                                                                        cgAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGc                 >  2:2091398/1‑93 (MQ=255)
                                                                          aataatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGccc               <  2:4282936/93‑1 (MQ=255)
                                                                            taatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCcgt             <  2:4687446/93‑1 (MQ=255)
                                                                             aatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCcgtc            >  2:3663045/1‑93 (MQ=255)
                                                                             aatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCcgtc            >  1:4605056/1‑93 (MQ=255)
                                                                              atGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCcgtcg           >  2:2423817/1‑93 (MQ=255)
                                                                               tGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTCGGTCAGCCcgtcgt          >  2:1858136/1‑93 (MQ=255)
                                                                                gCCGGACTGCCGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTg         <  2:3522923/93‑1 (MQ=255)
                                                                                gCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTg         >  2:407056/1‑93 (MQ=255)
                                                                                gCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTg         >  1:286586/1‑93 (MQ=255)
                                                                                 ccGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCACCTCGTGa        >  2:4011375/1‑93 (MQ=255)
                                                                                   ggACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATa      <  1:3579749/93‑1 (MQ=255)
                                                                                    gACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAg     >  2:4199939/1‑93 (MQ=255)
                                                                                    gACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAg     >  1:3233201/1‑93 (MQ=255)
                                                                                      cTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTa   >  2:4456847/1‑93 (MQ=255)
                                                                                       tGACGCCCTTTCTCCAGCAATCGCGACACTTTCAACCGTGTCAGGCCGAGCCGATCGCTGATCTCGCTCTGGGTCAGCCCGTGGTGATAGTaa  >  1:2006360/1‑93 (MQ=255)
                                                                                       tGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTaa  <  1:4458862/93‑1 (MQ=255)
                                                                                          |                                                                                          
TGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAA  >  NC_000913/1599012‑1599191

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