Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,348,109 C→T A378T (GCC→ACC)  yjdH ← putative 2‑component sensor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,348,1090CT97.5% 826.6 / 2.0 239A378T (GCC→ACC) yjdHputative 2‑component sensor protein
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base T (121/112);  total (122/117)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

Only 100 of 245 total aligned reads displayed.

CGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGCTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTCGTGGGAT  >  NC_000913/4348019‑4348200
                                                                                          |                                                                                           
cgcgGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTg                                                                                           <  1:4079788/93‑1 (MQ=255)
  cgGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGgttgtt                                                                                         >  2:3479633/1‑93 (MQ=255)
  cgGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGgttgtt                                                                                         <  1:455720/93‑1 (MQ=255)
  cgGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGgttgtt                                                                                         <  1:182710/93‑1 (MQ=255)
   ggTTAATCTTGCTGATTAAAAAACCCGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGttt                                                                                        >  2:847437/1‑93 (MQ=255)
   ggTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGttt                                                                                        <  1:1961900/93‑1 (MQ=255)
    gTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGtttt                                                                                       <  1:2576055/93‑1 (MQ=255)
     ttAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTg                                                                                      <  2:1839328/93‑1 (MQ=255)
     ttAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTg                                                                                      >  1:675697/1‑93 (MQ=255)
       aaTCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTgag                                                                                    >  2:1881263/1‑93 (MQ=255)
        aTCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTCGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTgaga                                                                                   <  2:3674757/93‑1 (MQ=255)
         tCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGaa                                                                                  <  2:1439079/93‑1 (MQ=255)
         tCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGaa                                                                                  <  1:572762/93‑1 (MQ=255)
         tCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGaa                                                                                  <  1:4357972/93‑1 (MQ=255)
           ttGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATg                                                                                >  1:3720124/1‑93 (MQ=255)
            tGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGt                                                                               <  2:1053443/93‑1 (MQ=255)
             gCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTa                                                                              >  2:1132309/1‑93 (MQ=255)
                gATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATc                                                                           <  1:4172477/93‑1 (MQ=255)
                gATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATc                                                                           >  1:4636597/1‑93 (MQ=255)
                   tAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTc                                                                        <  2:374676/93‑1 (MQ=255)
                   tAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTc                                                                        <  2:3932539/93‑1 (MQ=255)
                   tAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTc                                                                        <  2:1602335/93‑1 (MQ=255)
                   tAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTc                                                                        <  2:814881/93‑1 (MQ=255)
                     aaaaaCCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCa                                                                      >  1:1352276/1‑93 (MQ=255)
                     aaaaaCCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGCTGTTTTGAGAATGTAATCTTCCa                                                                      <  2:4115824/93‑1 (MQ=255)
                      aaaaCCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCaa                                                                     <  2:3271020/93‑1 (MQ=255)
                      aaaaCCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCaa                                                                     <  1:3968853/93‑1 (MQ=255)
                       aaaCCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAAc                                                                    <  2:1269040/93‑1 (MQ=255)
                       aaaCCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAAc                                                                    <  1:1684466/93‑1 (MQ=255)
                        aaCCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACt                                                                   <  1:247185/93‑1 (MQ=255)
                        aaCCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACt                                                                   >  1:2055338/1‑93 (MQ=255)
                         tccAGCGATAACCGGAGATGTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTGATTGGTTGTTGTGAGAATGTAATCTTCCAACTg                                                                  >  1:4209214/2‑93 (MQ=255)
                         aCCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTg                                                                  >  1:2182134/1‑93 (MQ=255)
                          ccAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGc                                                                 <  1:4581632/93‑1 (MQ=255)
                          ccAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGc                                                                 >  2:4649618/1‑93 (MQ=255)
                             gCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTta                                                              <  2:1029255/93‑1 (MQ=255)
                             gCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTta                                                              >  2:1996437/1‑93 (MQ=255)
                             gCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTta                                                              >  1:2985109/1‑93 (MQ=255)
                             gCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTta                                                              >  1:2082023/1‑93 (MQ=255)
                             gCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTta                                                              <  2:510702/93‑1 (MQ=255)
                             gCGATAACCGCAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTta                                                              <  2:2743304/93‑1 (MQ=255)
                               gATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTtata                                                            >  1:1252208/1‑93 (MQ=255)
                               gATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGAGAGTTATTGGCTGTTTTGAGAATGTAATCTTCCAACTGCTtata                                                            <  2:3332190/93‑1 (MQ=255)
                                 tAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAAc                                                          >  1:2273338/1‑93 (MQ=255)
                                 tAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAAc                                                          <  2:2898529/93‑1 (MQ=255)
                                 tAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTGATAGGTTGTTTTGAGAATGTAATCTTCCAACAGCTTATAAc                                                          >  2:4221537/1‑93 (MQ=255)
                                  aaCCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATCAct                                                         >  1:1183138/1‑93 (MQ=255)
                                   aCCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAActc                                                        <  2:3416685/93‑1 (MQ=255)
                                    ccGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAActct                                                       >  2:1183119/1‑93 (MQ=255)
                                      ggAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTc                                                     >  1:173654/1‑93 (MQ=255)
                                      ggAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTc                                                     <  2:2601095/93‑1 (MQ=255)
                                       gagaTTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCACCTGCTTATAACTCTTCa                                                    <  1:4494329/93‑1 (MQ=255)
                                       gagaTTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCa                                                    >  2:2722663/1‑93 (MQ=255)
                                        agaTTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAg                                                   >  1:4425970/1‑93 (MQ=255)
                                         ggtttGATCTTTCCCAGGAGAGAGCCCATCTCTTTCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGa                                                  <  1:774431/91‑1 (MQ=255)
                                         gaTTTGATCTTACCCAGCCGAGAGCCAATCTCTTCCTGATAGTTATTTGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGa                                                  >  2:3158982/1‑93 (MQ=255)
                                         gaTTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGa                                                  <  1:332024/93‑1 (MQ=255)
                                         gaTTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGa                                                  <  1:3267358/93‑1 (MQ=255)
                                               aTCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCaat                                            >  2:1421203/1‑93 (MQ=255)
                                                tCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCaata                                           >  2:4449749/1‑93 (MQ=255)
                                                tCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCaata                                           >  1:1882401/1‑93 (MQ=255)
                                                tCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCaata                                           >  2:4683312/1‑93 (MQ=255)
                                                 cTTGCCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCaataa                                          <  2:4579914/93‑1 (MQ=255)
                                                   tACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATc                                        <  1:2884713/93‑1 (MQ=255)
                                                    aCCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATcc                                       <  2:3009507/93‑1 (MQ=255)
                                                      ccAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCga                                     <  1:3706518/93‑1 (MQ=255)
                                                      ccAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCga                                     >  1:4006304/1‑93 (MQ=255)
                                                         gcagAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGaa                                  >  2:4381682/1‑93 (MQ=255)
                                                          cagAGAGCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAAt                                 >  1:3089149/1‑93 (MQ=255)
                                                           agagagCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATc                                <  1:2003321/93‑1 (MQ=255)
                                                            gagagCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATca                               <  1:3659060/93‑1 (MQ=255)
                                                            gagagCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATca                               >  1:1487749/1‑93 (MQ=255)
                                                            gagagCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATca                               >  1:24101/1‑93 (MQ=255)
                                                              gagCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATcaca                             <  2:2979105/93‑1 (MQ=255)
                                                              gagCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATcaca                             >  1:4735372/1‑93 (MQ=255)
                                                               agCCAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACAt                            >  1:861376/1‑93 (MQ=255)
                                                                gCCAATCTCTTCCTGATAGTTATTGGCTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATg                           >  1:724082/1‑93 (MQ=255)
                                                                 ccAATCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGc                          <  2:1743088/93‑1 (MQ=255)
                                                                    aTCTCTTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAAt                       >  2:2074424/1‑93 (MQ=255)
                                                                      ctctTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAAttt                     <  1:3926680/93‑1 (MQ=255)
                                                                      ctctTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAAttt                     <  2:2137147/93‑1 (MQ=255)
                                                                       tctTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTa                    >  2:1084204/1‑93 (MQ=255)
                                                                        ctTCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTAt                   >  2:569564/1‑93 (MQ=255)
                                                                         ttCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTAtt                  >  2:2330145/1‑93 (MQ=255)
                                                                          tCCTGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTc                 >  2:16813/1‑93 (MQ=255)
                                                                             tGATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATa              >  2:2412626/1‑93 (MQ=255)
                                                                              gATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATaa             <  2:3760638/93‑1 (MQ=255)
                                                                              gATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATaa             <  1:1778664/93‑1 (MQ=255)
                                                                              gATAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATaa             <  1:942745/93‑1 (MQ=255)
                                                                                tAGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAAt           >  1:1512320/1‑93 (MQ=255)
                                                                                 aGTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAAtt          <  1:4701926/93‑1 (MQ=255)
                                                                                  gTTATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTc         <  2:158378/93‑1 (MQ=255)
                                                                                    tATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTCGt       >  2:3646934/1‑93 (MQ=255)
                                                                                    tATTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTCGt       >  2:4493739/1‑93 (MQ=255)
                                                                                     aTTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTCGTg      <  1:229295/93‑1 (MQ=255)
                                                                                     aTTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTCGTg      >  1:3475009/1‑93 (MQ=255)
                                                                                     aTTGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTCGTg      <  1:4342186/93‑1 (MQ=255)
                                                                                     aTTGGTTGTTTTGAGAATGTAAGCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTCGTg      >  1:845804/1‑93 (MQ=255)
                                                                                      ttGGTTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTCGTgg     <  2:208947/93‑1 (MQ=255)
                                                                                         gCTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTCGTGGGAt  <  1:4103903/93‑1 (MQ=255)
                                                                                          |                                                                                           
CGCGGTTAATCTTGCTGATTAAAAAACCAGCGATAACCGGAGATTTGATCTTACCCAGCAGAGAGCCAATCTCTTCCTGATAGTTATTGGCTGTTTTGAGAATGTAATCTTCCAACTGCTTATAACTCTTCAGATGCAATAATCCGAGAATCACATGCAATTTATTCATAAATTCGTGGGAT  >  NC_000913/4348019‑4348200

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