Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,593,336 +G intergenic (‑23/‑208) yjjM ← / → yjjN hypothetical protein/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,593,3351.G99.6% 769.6 / ‑5.2 235intergenic (‑22/‑209)yjjM/yjjNhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (97/137);  minor base A (0/1);  total (97/138)
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.49e-01

Only 100 of 265 total aligned reads displayed.

TGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATA  >  NC_000913/4593243‑4593421
                                                                                             |                                                                                      
tGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGAt                                                                                         <  2:6043872/93‑1 (MQ=255)
tGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGAt                                                                                         <  2:5029224/93‑1 (MQ=255)
  cGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgg                                                                                      <  2:586562/93‑2 (MQ=255)
   gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga                                                                                     >  1:4340272/1‑91 (MQ=255)
   gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga                                                                                     >  2:5672418/1‑91 (MQ=255)
   gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga                                                                                     <  1:6480922/93‑3 (MQ=255)
   gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga                                                                                     >  2:3618365/1‑91 (MQ=255)
   gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga                                                                                     >  1:2558678/1‑91 (MQ=255)
    tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa                                                                                    >  2:4977524/1‑90 (MQ=255)
    tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa                                                                                    <  1:1808741/93‑4 (MQ=255)
    tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa                                                                                    >  2:2495025/1‑90 (MQ=255)
     cAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaag                                                                                   <  2:3874937/93‑5 (MQ=255)
     cAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaag                                                                                   <  2:257435/93‑5 (MQ=255)
         acgccgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTggg                                                                                <  2:4014354/89‑1 (MQ=255)
         acgccgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTggg                                                                                >  1:5334380/5‑93 (MQ=255)
          cgccgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTgggg                                                                               >  2:865782/4‑93 (MQ=255)
           gccgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGa                                                                              >  1:1163030/3‑93 (MQ=255)
             cgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAg                                                                            <  2:2390645/93‑1 (MQ=255)
             cgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAg                                                                            <  1:1937582/93‑1 (MQ=255)
              gcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAgg                                                                           >  2:1327865/1‑93 (MQ=255)
              gcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAgg                                                                           <  2:1708204/93‑1 (MQ=255)
                aTTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGa                                                                         >  2:2892665/1‑93 (MQ=255)
                 ttCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGat                                                                        >  2:523679/1‑93 (MQ=255)
                 ttCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGat                                                                        >  1:620322/1‑93 (MQ=255)
                  tCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGata                                                                       <  2:6380866/93‑1 (MQ=255)
                   cgcgTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGatat                                                                      >  2:4866578/1‑93 (MQ=255)
                      gTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAc                                                                   >  1:3801951/1‑93 (MQ=255)
                        aaGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAccc                                                                 <  2:671962/93‑1 (MQ=255)
                         aGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAcccc                                                                >  2:5559257/1‑93 (MQ=255)
                           tGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCcgc                                                              <  2:5324231/93‑1 (MQ=255)
                             gCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCgcgc                                                            >  2:4667211/1‑93 (MQ=255)
                              gtGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTATATTTCCGCCGGCGCTCTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGcc                                                           <  2:3014059/92‑1 (MQ=255)
                               tGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCa                                                          >  2:4249854/1‑93 (MQ=255)
                                gagaATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCat                                                         >  1:5088869/1‑93 (MQ=255)
                                  gaATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCatat                                                       <  1:6551906/93‑1 (MQ=255)
                                   aaTATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatc                                                      <  2:3295859/93‑1 (MQ=255)
                                    atatGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatca                                                     <  1:1197275/93‑1 (MQ=255)
                                    atatGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatca                                                     >  2:181388/1‑93 (MQ=255)
                                     tatGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatcat                                                    <  1:5921572/93‑1 (MQ=255)
                                        gACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCaa                                                 <  1:1603693/93‑1 (MQ=255)
                                        gACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCaa                                                 <  1:3352643/93‑1 (MQ=255)
                                        gACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCaa                                                 <  1:4486537/93‑1 (MQ=255)
                                        gACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCaa                                                 <  1:2776645/93‑1 (MQ=255)
                                         aCGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAAt                                                <  1:2379510/93‑1 (MQ=255)
                                         aCGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAAt                                                <  2:2035259/93‑1 (MQ=255)
                                         aCGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAAt                                                <  1:5232313/93‑1 (MQ=255)
                                          cGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATg                                               >  2:3493477/1‑93 (MQ=255)
                                          cGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATg                                               >  1:297499/1‑93 (MQ=255)
                                             aCTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGAcc                                            >  1:586426/1‑93 (MQ=255)
                                             aCTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGAcc                                            >  2:3181390/1‑93 (MQ=255)
                                             aCTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGAcc                                            >  1:2808419/1‑93 (MQ=255)
                                              cTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGTGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCt                                           >  2:5532972/1‑93 (MQ=255)
                                               tgtgGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGTCCCCGCTCCATATCATCAATGACCTg                                          <  2:1884441/93‑1 (MQ=255)
                                               tgtgGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTg                                          <  1:6127111/93‑1 (MQ=255)
                                               tgtgGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTg                                          >  2:3567493/1‑93 (MQ=255)
                                                gtggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTgg                                         >  1:6680774/1‑93 (MQ=255)
                                                  ggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGtt                                       <  2:5416002/93‑1 (MQ=255)
                                                  ggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGtt                                       <  2:2128127/93‑1 (MQ=255)
                                                    tgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTaa                                     <  2:3527161/93‑1 (MQ=255)
                                                     gCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAAt                                    <  1:5015157/93‑1 (MQ=255)
                                                      cGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATc                                   >  1:3435442/1‑93 (MQ=255)
                                                      cGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATc                                   >  2:4186229/1‑93 (MQ=255)
                                                       ggCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATca                                  <  1:3187091/93‑1 (MQ=255)
                                                        gCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATcac                                 <  2:1503589/93‑1 (MQ=255)
                                                        gCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATcac                                 >  1:1039214/1‑93 (MQ=255)
                                                        gCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATcac                                 >  2:2243554/1‑93 (MQ=255)
                                                         cTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATcaca                                <  1:4146193/93‑1 (MQ=255)
                                                         cTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATcaca                                <  1:5559916/93‑1 (MQ=255)
                                                          tGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAt                               <  1:5894091/93‑1 (MQ=255)
                                                          tGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAt                               <  1:2355077/93‑1 (MQ=255)
                                                           gATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAtt                              <  1:1296953/93‑1 (MQ=255)
                                                           gATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAtt                              >  2:3834617/1‑93 (MQ=255)
                                                             ttttgtATATTGCCGCCAGTGCTCTTTGGGATGGAAGTGGGGAAGGTATATGACCCCGCGGCATATCATCAATGACCTGGTTAATCGCATtgt                            <  1:4772088/91‑1 (MQ=255)
                                                             tgttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATtgt                            >  1:3026206/1‑93 (MQ=255)
                                                             tgttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATtgt                            <  2:6308122/93‑1 (MQ=255)
                                                              gttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATtgtg                           <  2:5899360/93‑1 (MQ=255)
                                                              gttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGAGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATtgtg                           <  1:1059975/93‑1 (MQ=255)
                                                                tgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGc                         <  1:4083055/93‑1 (MQ=255)
                                                                tgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGc                         >  2:3051035/1‑93 (MQ=255)
                                                                 gtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCg                        >  1:1442911/1‑93 (MQ=255)
                                                                 gtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCg                        <  2:1145209/93‑1 (MQ=255)
                                                                 gtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCg                        >  2:699874/1‑93 (MQ=255)
                                                                   aCATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTa                      >  2:3363393/1‑93 (MQ=255)
                                                                      tttCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAt                   <  2:80029/93‑1 (MQ=255)
                                                                      tttCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAt                   <  1:5723224/93‑1 (MQ=255)
                                                                      tttCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAt                   >  2:1361140/1‑93 (MQ=255)
                                                                       ttCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAtt                  <  1:3538862/93‑1 (MQ=255)
                                                                          cGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTtg               >  2:2309731/1‑93 (MQ=255)
                                                                            ccAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTtgtg             >  1:3987318/1‑93 (MQ=255)
                                                                              aGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGaa           <  2:4042457/93‑1 (MQ=255)
                                                                               gcgcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGTAc          >  1:2116843/1‑93 (MQ=255)
                                                                               gcgcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAAc          >  2:1954514/1‑93 (MQ=255)
                                                                               gcgcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAAc          <  1:5740572/93‑1 (MQ=255)
                                                                                  cACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGAc       <  2:4770807/93‑1 (MQ=255)
                                                                                   aCTTTGCGATGGCAGTGGGGAACGGATATGAGCGCGCGCCATATCATCAATGACCTGGTTAATCACATTGTAGCGTAAATTTTGTGAACGACt      >  1:891165/1‑93 (MQ=255)
                                                                                   aCTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACt      <  1:5408709/93‑1 (MQ=255)
                                                                                   aCTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACt      <  1:459852/93‑1 (MQ=255)
                                                                                     tttGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCa    >  1:1513089/1‑93 (MQ=255)
                                                                                       tGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCata  <  2:1042302/93‑1 (MQ=255)
                                                                                       tGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCGGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCata  <  1:293338/93‑1 (MQ=255)
                                                                                             |                                                                                      
TGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATA  >  NC_000913/4593243‑4593421

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