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,1020TC99.1% 345.2 / ‑4.5 112H55R (CAT→CGT) ydeWputative transcriptional regulator, sorC family
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (67/44);  total (68/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 = 7.11e-01

Only 100 of 113 total aligned reads displayed.

GCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAA  >  NC_000913/1599013‑1599191
                                                                                         |                                                                                         
gCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGcc                                                                                        <  2:1477442/93‑1 (MQ=255)
 cAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGccc                                                                                       <  2:107302/93‑1 (MQ=255)
   gCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCtt                                                                                     <  1:4507443/93‑1 (MQ=255)
   gCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCtt                                                                                     <  1:4310804/93‑1 (MQ=255)
    cGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCttt                                                                                    <  1:1131857/93‑1 (MQ=255)
     gAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtc                                                                                   >  2:1460330/1‑93 (MQ=255)
     gAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtc                                                                                   <  1:559411/93‑1 (MQ=255)
     gAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtc                                                                                   >  1:4619999/1‑93 (MQ=255)
     gAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtc                                                                                   >  2:2567318/1‑93 (MQ=255)
     gAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtc                                                                                   >  1:3289345/1‑93 (MQ=255)
     gAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtc                                                                                   <  2:665342/93‑1 (MQ=255)
      aaaaCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtct                                                                                  <  1:3171758/93‑1 (MQ=255)
      aaaaCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtct                                                                                  >  1:2279335/1‑93 (MQ=255)
       aaaCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtctc                                                                                 >  2:1026842/1‑93 (MQ=255)
       aaaCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtctc                                                                                 <  1:1953946/93‑1 (MQ=255)
       aaaCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTtctc                                                                                 >  2:374762/1‑93 (MQ=255)
        aaCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTcc                                                                                >  2:15951/1‑93 (MQ=255)
         aCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCa                                                                               <  2:823500/93‑1 (MQ=255)
          cTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGCCGCCCTTTCTCCAg                                                                              >  1:4459360/1‑93 (MQ=255)
          cTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAg                                                                              >  2:442822/1‑93 (MQ=255)
           tGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGc                                                                             >  2:260262/1‑93 (MQ=255)
            gacgacGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCa                                                                            >  2:493129/1‑93 (MQ=255)
            gacgacGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCa                                                                            >  2:4428278/1‑93 (MQ=255)
             acgacgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCaa                                                                           <  1:171634/93‑1 (MQ=255)
              cgacgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAAt                                                                          <  2:3937067/93‑1 (MQ=255)
               gacgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATc                                                                         >  2:1865626/1‑93 (MQ=255)
               gacgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATc                                                                         >  1:293419/1‑93 (MQ=255)
                acgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATcg                                                                        <  2:1213246/93‑1 (MQ=255)
                 cgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATcgc                                                                       <  1:3669785/93‑1 (MQ=255)
                 cgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATcgc                                                                       >  2:2988584/1‑93 (MQ=255)
                 cgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGACGGACTGACGCCCTTTCTCCAGCAATcgc                                                                       >  2:2069893/1‑93 (MQ=255)
                   tAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGa                                                                     <  1:1204742/93‑1 (MQ=255)
                   tAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGa                                                                     >  1:2842987/1‑93 (MQ=255)
                   tAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGa                                                                     >  1:3202719/1‑93 (MQ=255)
                    aaTTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACTCGCTTTCTCCAGCAATCGCGac                                                                    <  2:2868357/93‑1 (MQ=255)
                    aaTTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGac                                                                    >  2:3285135/1‑93 (MQ=255)
                    aaTTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGac                                                                    <  2:4309986/93‑1 (MQ=255)
                    aaTTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCATTTCTCCAGCAATCGCGac                                                                    >  1:4631970/1‑93 (MQ=255)
                     aTTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGaca                                                                   >  1:1479179/1‑93 (MQ=255)
                     aTTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGaca                                                                   >  2:2484934/1‑93 (MQ=255)
                       tGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACt                                                                 >  2:1545148/1‑93 (MQ=255)
                          gTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCATTTCTCCAGCAATCGCGACACTTTc                                                              >  1:2312953/1‑93 (MQ=255)
                            ttCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCaa                                                            >  2:4462521/1‑93 (MQ=255)
                             tCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCaaa                                                           <  2:2822776/93‑1 (MQ=255)
                             tCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCaaa                                                           <  1:979522/93‑1 (MQ=255)
                             tCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCaaa                                                           <  2:2992789/93‑1 (MQ=255)
                                tatTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACgt                                                        >  2:3866092/1‑93 (MQ=255)
                                tatTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACgt                                                        <  1:1328481/93‑1 (MQ=255)
                                tatTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACgt                                                        >  2:3314900/1‑93 (MQ=255)
                                  ttCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACgtgt                                                      >  2:82578/1‑93 (MQ=255)
                                    ccAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCa                                                    <  1:4006393/93‑1 (MQ=255)
                                      aGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGACCTTTCTCCAGCAATCGCGACACTTTCACACGTGTCAgc                                                  >  2:4660365/1‑92 (MQ=255)
                                       gaaaGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGc                                                 >  2:2688254/4‑93 (MQ=255)
                                       gACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGc                                                 >  2:2580063/1‑93 (MQ=255)
                                         cAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCg                                               >  2:3914358/1‑93 (MQ=255)
                                           gCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCgag                                             >  1:3575484/1‑93 (MQ=255)
                                           gCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCgag                                             >  2:248377/1‑93 (MQ=255)
                                            ccTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCgaga                                            <  1:3121835/93‑1 (MQ=255)
                                            ccTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCgaga                                            <  2:3713762/93‑1 (MQ=255)
                                            ccTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCgaga                                            <  1:746031/93‑1 (MQ=255)
                                              ttCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACg                                          >  1:4460159/1‑93 (MQ=255)
                                               tCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGa                                         <  1:1196238/93‑1 (MQ=255)
                                               tCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGa                                         >  2:326740/1‑93 (MQ=255)
                                                 aaaGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATc                                       <  2:2777886/93‑1 (MQ=255)
                                                 aaaGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATc                                       >  1:930182/1‑93 (MQ=255)
                                                  aaGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCg                                      >  2:4488162/1‑93 (MQ=255)
                                                   aGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGc                                     >  2:738688/1‑93 (MQ=255)
                                                   aGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGc                                     <  2:3495041/93‑1 (MQ=255)
                                                      gagaATTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGa                                  >  1:157760/1‑93 (MQ=255)
                                                        gaATTAATCTGTACGCGAATAATGCCGGACTGACACCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGAtc                                <  2:2800769/93‑1 (MQ=255)
                                                         aaTTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGAtct                               >  1:4294256/1‑93 (MQ=255)
                                                         aaTTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGAtct                               >  2:3560500/1‑93 (MQ=255)
                                                          aTTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGAtctc                              >  2:39385/1‑93 (MQ=255)
                                                          aTTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGAtctc                              <  2:3153346/93‑1 (MQ=255)
                                                          aTTAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACCATAGCTGAtctc                              >  2:2251033/1‑93 (MQ=255)
                                                           ttAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCg                             >  1:3075416/1‑93 (MQ=255)
                                                            tAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGc                            <  2:3458286/93‑1 (MQ=255)
                                                            tAATCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGc                            >  1:4478119/1‑93 (MQ=255)
                                                             aaTCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGct                           >  1:3079777/1‑93 (MQ=255)
                                                              aTCTGTACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGctc                          <  2:4418716/93‑1 (MQ=255)
                                                                   tACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGt                     <  2:1132267/93‑1 (MQ=255)
                                                                   tACGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGt                     <  2:1112305/93‑1 (MQ=255)
                                                                    aCGCGAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTc                    <  1:13899/93‑1 (MQ=255)
                                                                     cgcgAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCa                   >  2:2396774/1‑93 (MQ=255)
                                                                        gAATAATGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGcc                >  1:424728/1‑93 (MQ=255)
                                                                         aataatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGccc               >  1:2437479/1‑93 (MQ=255)
                                                                         aataatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGccc               <  1:2631323/93‑1 (MQ=255)
                                                                         aataatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGccc               <  2:4082931/93‑1 (MQ=255)
                                                                         aataatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCAGGGTCATccc               >  2:4164560/1‑93 (MQ=255)
                                                                          ataatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCcg              >  1:3294765/1‑93 (MQ=255)
                                                                           taatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCcgt             >  2:2116809/1‑93 (MQ=255)
                                                                            aatGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCcgtc            <  1:2382851/93‑1 (MQ=255)
                                                                              tGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCcgtcgt          >  2:1642393/1‑93 (MQ=255)
                                                                              tGCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCcgtcgt          >  2:758785/1‑93 (MQ=255)
                                                                               gCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTg         >  1:4488631/1‑93 (MQ=255)
                                                                               gCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTg         >  1:398302/1‑93 (MQ=255)
                                                                               gCCGGACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTg         >  1:587059/1‑93 (MQ=255)
                                                                                  ggACTGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATa      <  1:3406378/93‑1 (MQ=255)
                                                                                   gACTGCCGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAg     <  2:128478/93‑1 (MQ=255)
                                                                                      tGACGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTaa  >  2:2464480/1‑93 (MQ=255)
                                                                                         |                                                                                         
GCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAA  >  NC_000913/1599013‑1599191

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