Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,210,781 T→C M68V (ATG→GTG)  b2125 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,210,7810TC99.3% 432.5 / ‑4.8 139M68V (ATG→GTG) b2125hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (80/58);  total (81/58)
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.05e-01

Only 100 of 144 total aligned reads displayed.

ATCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCATCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATCCC  >  NC_000913/2210692‑2210868
                                                                                         |                                                                                       
aTCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCAcccc                                                                                      <  1:4570129/93‑1 (MQ=255)
aTCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCAcccc                                                                                      <  1:853212/93‑1 (MQ=255)
aTCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCAcccc                                                                                      >  1:2138742/1‑93 (MQ=255)
 tCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACcccc                                                                                     >  1:4404055/1‑93 (MQ=255)
  cAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCcca                                                                                    >  2:58580/1‑93 (MQ=255)
    aaGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacc                                                                                  >  1:1295409/1‑93 (MQ=255)
    aaGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacc                                                                                  >  2:1629091/1‑93 (MQ=255)
     aGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacca                                                                                 <  1:579756/93‑1 (MQ=255)
     aGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCccacca                                                                                 >  1:167756/1‑93 (MQ=255)
      ggCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCAt                                                                                >  1:749502/1‑93 (MQ=255)
      ggCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCAt                                                                                >  2:3276938/1‑93 (MQ=255)
        cATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCAttt                                                                              >  1:1386996/1‑93 (MQ=255)
          tGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTcc                                                                            >  1:4019680/1‑93 (MQ=255)
          tGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTcc                                                                            >  2:19332/1‑93 (MQ=255)
          tGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTcc                                                                            >  2:2198890/1‑93 (MQ=255)
           gTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCa                                                                           >  2:629315/1‑93 (MQ=255)
           gTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCa                                                                           <  2:2581352/93‑1 (MQ=255)
               ttcAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGAcc                                                                       >  1:191186/1‑93 (MQ=255)
                  aaaGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACt                                                                    >  2:528422/1‑93 (MQ=255)
                   aaGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTg                                                                   >  2:1434720/1‑93 (MQ=255)
                    aGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGa                                                                  >  1:2241917/1‑93 (MQ=255)
                    aGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGa                                                                  >  1:2116203/1‑93 (MQ=255)
                    aGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGa                                                                  <  1:3737379/93‑1 (MQ=255)
                    aGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGa                                                                  <  1:2907215/93‑1 (MQ=255)
                    aGGCTGTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCCCCCCCACCATTTCCAGACCACTGa                                                                  <  2:2267545/93‑1 (MQ=255)
                     ggCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCATCCCCACCATTTCCAGACCACTGAt                                                                 >  1:1394727/1‑93 (MQ=255)
                     ggCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGAt                                                                 <  1:2409305/93‑1 (MQ=255)
                     ggCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGAt                                                                 >  2:4676621/1‑93 (MQ=255)
                     agCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGAt                                                                 <  2:3773983/92‑1 (MQ=255)
                      gCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATg                                                                <  2:3555209/93‑1 (MQ=255)
                       cTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgc                                                               <  1:2271034/93‑1 (MQ=255)
                        tttAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgcg                                                              <  2:1001688/93‑1 (MQ=255)
                         ttAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgcgc                                                             <  1:2784991/93‑1 (MQ=255)
                         ttAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATgcgc                                                             >  2:1743701/1‑93 (MQ=255)
                          tAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGcgcg                                                            >  2:1328966/1‑93 (MQ=255)
                           aaTTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCgg                                                           >  1:1426321/1‑93 (MQ=255)
                            aTTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGc                                                          >  1:4449829/1‑93 (MQ=255)
                             ttGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCa                                                         <  1:1506228/93‑1 (MQ=255)
                             ttGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCa                                                         >  1:678527/1‑93 (MQ=255)
                               gCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATc                                                       <  1:68993/93‑1 (MQ=255)
                               gCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATc                                                       <  2:3846574/93‑1 (MQ=255)
                                 gTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTg                                                     >  1:125539/1‑93 (MQ=255)
                                 gTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTg                                                     >  1:3580002/1‑93 (MQ=255)
                                 gTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTg                                                     <  1:1592695/93‑1 (MQ=255)
                                   aTTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGa                                                   >  1:3321113/1‑93 (MQ=255)
                                    ttCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGat                                                  >  2:1266919/1‑93 (MQ=255)
                                    ttCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGat                                                  >  1:3002173/1‑93 (MQ=255)
                                     tCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGata                                                 <  2:929045/93‑1 (MQ=255)
                                         cAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCga                                             <  1:3440658/93‑1 (MQ=255)
                                         cAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCga                                             <  1:926290/93‑1 (MQ=255)
                                            aCGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGGTGCGCGGCATCGGGATATCGAGaa                                          >  2:2035587/1‑93 (MQ=255)
                                             cGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGaaa                                         <  2:4386475/93‑1 (MQ=255)
                                             cGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGaaa                                         <  1:1695914/93‑1 (MQ=255)
                                             cGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGaaa                                         >  1:983196/1‑93 (MQ=255)
                                             agCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCCCTGATGCGCGGCATCTGGATATCGAGaaa                                         <  1:4166050/92‑1 (MQ=255)
                                              gCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAAc                                        >  2:1532180/1‑93 (MQ=255)
                                                aGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAg                                      >  2:2442008/1‑93 (MQ=255)
                                                aGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAg                                      <  1:1897636/93‑1 (MQ=255)
                                                 gTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGc                                     >  2:1074285/1‑93 (MQ=255)
                                                  tGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGca                                    >  1:582832/1‑93 (MQ=255)
                                                  tGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGca                                    <  1:611697/93‑1 (MQ=255)
                                                    agaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGcaca                                  <  1:3782276/93‑1 (MQ=255)
                                                    agaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGcaca                                  <  1:4156774/93‑1 (MQ=255)
                                                     gaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACAt                                 >  2:254447/1‑93 (MQ=255)
                                                     gaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACAt                                 >  2:3591986/1‑93 (MQ=255)
                                                     gaAAAACAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACAt                                 >  1:2987408/1‑93 (MQ=255)
                                                       aaaaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATcc                               >  1:4064752/1‑93 (MQ=255)
                                                         aaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCgg                             >  2:3134942/1‑93 (MQ=255)
                                                         aaCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCgg                             <  2:4323483/93‑1 (MQ=255)
                                                          aCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCggg                            >  2:1223576/1‑93 (MQ=255)
                                                          aCAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCggg                            >  2:159339/1‑93 (MQ=255)
                                                           cAATATACGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGgc                           >  1:4597177/1‑93 (MQ=255)
                                                              tataCGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCTGGCGCa                        >  1:1196137/1‑93 (MQ=255)
                                                              tataCGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCa                        >  2:882364/1‑93 (MQ=255)
                                                              tataCGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCa                        <  1:3504469/93‑1 (MQ=255)
                                                               ataCGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAg                       >  2:1131889/1‑93 (MQ=255)
                                                                  cGGGCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGttt                    >  2:4525469/1‑93 (MQ=255)
                                                                   gggCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTa                   <  2:3293326/93‑1 (MQ=255)
                                                                   gggCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTa                   <  1:1280859/93‑1 (MQ=255)
                                                                   gggCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACAGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTa                   <  2:3293452/93‑1 (MQ=255)
                                                                    ggCGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTAt                  <  2:1175389/93‑1 (MQ=255)
                                                                      cGATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGc                >  1:3689634/1‑93 (MQ=255)
                                                                       gATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCa               <  1:3828849/93‑1 (MQ=255)
                                                                       gATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCa               <  1:306419/93‑1 (MQ=255)
                                                                       gATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCa               >  2:2851263/1‑93 (MQ=255)
                                                                       gATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCa               >  2:115200/1‑93 (MQ=255)
                                                                       gATGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCa               >  2:1293203/1‑93 (MQ=255)
                                                                        aTGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCAc              >  1:3576730/1‑93 (MQ=255)
                                                                        aTGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCAc              >  2:4177359/1‑93 (MQ=255)
                                                                         tGTTCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCAcc             <  1:1591638/93‑1 (MQ=255)
                                                                           ttCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACcgc           >  1:1348449/1‑93 (MQ=255)
                                                                           ttCCGGGTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACcgc           <  1:2777299/93‑1 (MQ=255)
                                                                               gggTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATACGGGCGCAGTTTATGCACCGCGCCg       >  1:1514579/1‑93 (MQ=255)
                                                                                ggTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGa      <  2:4064731/93‑1 (MQ=255)
                                                                                 gTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGAt     >  1:3274967/1‑93 (MQ=255)
                                                                                 gTCAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGAt     <  1:1991292/93‑1 (MQ=255)
                                                                                   cAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATcc   <  1:1557111/93‑1 (MQ=255)
                                                                                   cAAGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATcc   <  1:4088617/93‑1 (MQ=255)
                                                                                    aaGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATccc  <  1:1458373/93‑1 (MQ=255)
                                                                                    aaGCACCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATccc  <  1:2452194/93‑1 (MQ=255)
                                                                                         |                                                                                       
ATCAAAGGCATGTTCTTCAAAGGCTTTAATTGCGTATTCGTCAAACGCAGTGAGAAAAACAATATACGGGCGATGTTCCGGGTCAAGCATCCCCACCATTTCCAGACCACTGATGCGCGGCATCTGGATATCGAGAAACAGCACATCCGGGCGCAGTTTATGCACCGCGCCGATCCC  >  NC_000913/2210692‑2210868

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