Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,259,162 A→G intergenic (+277/‑86) yjbM → / → yjbN hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,259,1620AG99.5% 678.1 / ‑4.9 211intergenic (+277/‑86)yjbM/yjbNhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (105/104);  total (106/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

Only 100 of 228 total aligned reads displayed.

TTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGC  >  NC_000913/4259072‑4259251
                                                                                          |                                                                                         
ttattTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTatggtatggt                                                                                         >  1:3204227/1‑93 (MQ=255)
 tattTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTAtggtatggta                                                                                        >  1:2391866/1‑93 (MQ=255)
  attTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTAc                                                                                       >  1:4071871/1‑93 (MQ=255)
  attTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTAc                                                                                       >  1:1831076/1‑93 (MQ=255)
   tttCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTAcc                                                                                      >  2:1435130/1‑93 (MQ=255)
    ttCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTAccc                                                                                     >  2:3245496/1‑93 (MQ=255)
      ccATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCg                                                                                   <  1:3610887/93‑1 (MQ=255)
        aTCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCt                                                                                 >  2:2895876/1‑93 (MQ=255)
        aTCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCt                                                                                 >  1:1363427/1‑93 (MQ=255)
        aTCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCt                                                                                 <  2:4410604/93‑1 (MQ=255)
        aTCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCt                                                                                 >  2:881378/1‑93 (MQ=255)
          cACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAt                                                                               >  2:924329/1‑93 (MQ=255)
          cACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAt                                                                               <  1:4016280/93‑1 (MQ=255)
           aCTTCGTGAGAGCTTCATGTCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAtt                                                                              >  2:4565781/1‑93 (MQ=255)
           aCTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAtt                                                                              <  2:2993462/93‑1 (MQ=255)
           aCTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAtt                                                                              >  2:662654/1‑93 (MQ=255)
           aCTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAtt                                                                              >  1:3908977/1‑93 (MQ=255)
            cTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTg                                                                             >  2:141397/1‑93 (MQ=255)
             ttCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGa                                                                            <  1:2371278/93‑1 (MQ=255)
             ttCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGa                                                                            >  2:523883/1‑93 (MQ=255)
             ttCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGa                                                                            >  1:910682/1‑93 (MQ=255)
             ttCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGa                                                                            <  1:1003062/93‑1 (MQ=255)
              tCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAg                                                                           >  1:3789450/1‑93 (MQ=255)
                gTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGcc                                                                         >  2:2792898/1‑93 (MQ=255)
                 tGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGcct                                                                        >  2:2082749/1‑93 (MQ=255)
                 tGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGcct                                                                        >  2:2582023/1‑93 (MQ=255)
                   agagCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGcctcc                                                                      <  2:901327/93‑1 (MQ=255)
                      gCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGa                                                                   >  1:2472897/1‑93 (MQ=255)
                       cTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGaa                                                                  <  2:160612/93‑1 (MQ=255)
                       cTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGaa                                                                  <  1:2339248/93‑1 (MQ=255)
                        ttCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAAt                                                                 >  1:1918397/1‑93 (MQ=255)
                        ttCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAAt                                                                 >  1:4537560/1‑93 (MQ=255)
                         tCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATa                                                                >  2:4281822/1‑93 (MQ=255)
                           aTGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGt                                                              >  1:2901602/1‑93 (MQ=255)
                           aTGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGt                                                              <  1:3375388/93‑1 (MQ=255)
                             gCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGa                                                            <  2:3460506/93‑1 (MQ=255)
                              ccTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGAt                                                           <  2:1954699/93‑1 (MQ=255)
                               cTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATg                                                          >  1:4615806/1‑93 (MQ=255)
                                ttGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGc                                                         >  1:1346280/1‑93 (MQ=255)
                                ttGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGc                                                         <  2:770161/93‑1 (MQ=255)
                                   agagGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGa                                                      <  1:197562/93‑1 (MQ=255)
                                     agGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAat                                                    <  1:3510072/93‑1 (MQ=255)
                                      ggATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAata                                                   >  2:2979907/1‑93 (MQ=255)
                                       gATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataa                                                  >  2:448345/1‑93 (MQ=255)
                                       gATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataa                                                  <  2:1119566/93‑1 (MQ=255)
                                        aTCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataac                                                 <  1:113554/93‑1 (MQ=255)
                                         tctcAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataaca                                                <  2:2191645/93‑1 (MQ=255)
                                         tctcAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataaca                                                <  1:1551701/93‑1 (MQ=255)
                                         tctcAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataaca                                                >  1:3231941/1‑93 (MQ=255)
                                         tctcAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataaca                                                <  1:660880/93‑1 (MQ=255)
                                         tctcAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataaca                                                <  2:2290199/93‑1 (MQ=255)
                                          ctcAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacat                                               >  1:477588/1‑93 (MQ=255)
                                            cAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacataa                                             <  2:3561266/93‑1 (MQ=255)
                                            cAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacataa                                             <  1:3646304/93‑1 (MQ=255)
                                             aaTTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacataac                                            <  1:1655340/93‑1 (MQ=255)
                                             aaTTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacataac                                            <  2:600383/93‑1 (MQ=255)
                                             aaTTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacataac                                            <  1:3403506/93‑1 (MQ=255)
                                              aTTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAAcc                                           >  1:1208872/1‑93 (MQ=255)
                                              aTTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAAcc                                           <  2:3143345/93‑1 (MQ=255)
                                              aTTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAAcc                                           <  2:4596324/93‑1 (MQ=255)
                                              aTTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATACTGAATAACATAAcc                                           >  1:1239970/1‑93 (MQ=255)
                                               ttttCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAAccc                                          <  2:1774534/93‑1 (MQ=255)
                                                tttCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAAcccc                                         >  1:4232279/1‑93 (MQ=255)
                                                tttCTTTTACAATGAGACATGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAAcccc                                         >  2:4578025/1‑93 (MQ=255)
                                                    ttttGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGa                                     >  1:1506666/1‑93 (MQ=255)
                                                    ttttGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGa                                     <  1:787662/93‑1 (MQ=255)
                                                     tttGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGat                                    <  1:2688226/93‑1 (MQ=255)
                                                        gCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGatata                                 <  2:3674326/93‑1 (MQ=255)
                                                        gCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGatata                                 <  2:250661/93‑1 (MQ=255)
                                                          aaTGAGACAGGCGCTTCCTGTTGTTATGGTATGGTCCCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATc                               <  1:255333/93‑1 (MQ=255)
                                                          aaTGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATc                               <  2:788860/93‑1 (MQ=255)
                                                          aaTGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATc                               >  1:2451317/1‑93 (MQ=255)
                                                           aTGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCg                              >  1:1061354/1‑93 (MQ=255)
                                                            tGAGACAAGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGa                             >  2:4453739/1‑93 (MQ=255)
                                                             gagaCAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGAt                            <  1:1086460/93‑1 (MQ=255)
                                                             gagaCAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGAt                            <  2:433028/93‑1 (MQ=255)
                                                               gaCAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATaa                          >  2:721822/1‑93 (MQ=255)
                                                               gaCAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATaa                          <  1:698407/93‑1 (MQ=255)
                                                               gaCAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATaa                          <  1:4270024/93‑1 (MQ=255)
                                                                  aGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAaat                       <  2:4391960/93‑1 (MQ=255)
                                                                   ggCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAaata                      <  1:2953882/93‑1 (MQ=255)
                                                                    gcgcTGCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACACCATGATATATCGATAAaataa                     <  1:2992198/93‑1 (MQ=255)
                                                                      gcTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAAtc                   <  1:157968/93‑1 (MQ=255)
                                                                      gcTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAAtc                   >  1:176078/1‑93 (MQ=255)
                                                                      gcTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAAtc                   <  2:4014717/93‑1 (MQ=255)
                                                                       cTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAAtct                  <  1:4105198/93‑1 (MQ=255)
                                                                        ttCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATCATATATCGATAAAATAAtctc                 <  1:4647973/93‑1 (MQ=255)
                                                                        ttCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGAGATGCTGAATAACATAACCCCATGATATATCGATAAAATAAtctc                 >  1:3943497/1‑93 (MQ=255)
                                                                         tCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCGCATGATATATCGATAAAATAATctct                >  2:3796515/1‑93 (MQ=255)
                                                                           cTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTAc              <  1:2209090/93‑1 (MQ=255)
                                                                           cTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTAc              <  2:3846674/93‑1 (MQ=255)
                                                                            tgttgtTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACa             >  2:3993870/1‑93 (MQ=255)
                                                                             gttgttATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACAt            >  2:2321512/1‑93 (MQ=255)
                                                                             gttgttATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACAt            >  1:726785/1‑93 (MQ=255)
                                                                               tgttATCGTATGTCACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATACCATAGCCCCATGATACATCGATAAAACAATCTCTACGttt          <  2:2361012/93‑1 (MQ=255)
                                                                                 ttATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGa        <  2:1258618/93‑1 (MQ=255)
                                                                                 ttATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGa        >  1:816151/1‑93 (MQ=255)
                                                                                  tatggtatggTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGaa       >  1:2025705/1‑93 (MQ=255)
                                                                                   atggtatggtACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGaaa      >  2:4070815/1‑93 (MQ=255)
                                                                                    tggtatggtaCCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGaaaa     <  1:1722027/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGC  >  NC_000913/4259072‑4259251

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: