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.7% 1203.2 / ‑4.4 346intergenic (‑22/‑209)yjjM/yjjNhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (175/170);  total (175/171)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 384 total aligned reads displayed.

TGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATAT  >  NC_000913/4593243‑4593422
                                                                                             |                                                                                       
tGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGAt                                                                                          >  2:2539127/1‑93 (MQ=255)
 gCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATg                                                                                        >  2:3934901/1‑93 (MQ=255)
  cGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgg                                                                                       <  2:3671650/93‑2 (MQ=255)
   gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga                                                                                      <  2:2824722/93‑3 (MQ=255)
   gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga                                                                                      <  2:5050556/93‑3 (MQ=255)
     cAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaag                                                                                    <  2:5859705/93‑5 (MQ=255)
            ccgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGaa                                                                              <  2:5448238/92‑1 (MQ=255)
             cgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAg                                                                             <  1:3467491/93‑1 (MQ=255)
              gcATTCGCGTAAGTGGCTGCGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAgg                                                                            >  1:8533456/1‑93 (MQ=255)
              gcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAgg                                                                            <  2:3395278/93‑1 (MQ=255)
               cATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAggg                                                                           <  2:5515492/93‑1 (MQ=255)
                aTTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGa                                                                          <  1:8110215/93‑1 (MQ=255)
                aTTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGa                                                                          >  2:419527/1‑93 (MQ=255)
                 ttCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGat                                                                         >  2:3001092/1‑93 (MQ=255)
                  tCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTgggg                                                                                <  1:4762208/85‑1 (MQ=255)
                  tCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTgggg                                                                                >  2:2780811/1‑85 (MQ=255)
                  tCGCGTAAGTGGCGGAGAATATGACTCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGata                                                                        <  1:1089594/93‑1 (MQ=255)
                     cgTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGa                                                                     >  2:6326584/1‑93 (MQ=255)
                      gTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAc                                                                    >  2:752003/1‑93 (MQ=255)
                       tAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAcc                                                                   >  1:3373946/1‑93 (MQ=255)
                       tAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAcc                                                                   >  2:6021485/1‑93 (MQ=255)
                       tAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGGGGGGAAGGGATATGAcc                                                                   >  2:1608348/1‑93 (MQ=255)
                         aGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAcccc                                                                 <  2:6473999/93‑1 (MQ=255)
                           tGGCTGAGAATATGACGCACTGTAGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCcgc                                                               <  2:4761973/93‑1 (MQ=255)
                            ggCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCcgcg                                                              >  1:542497/1‑93 (MQ=255)
                             gCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCgcgc                                                             >  1:1444885/1‑93 (MQ=255)
                              cTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGcc                                                            >  1:1046206/1‑93 (MQ=255)
                                  gaATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCatat                                                        >  1:519527/1‑93 (MQ=255)
                                  gaATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCatat                                                        <  2:8482754/93‑1 (MQ=255)
                                   aaTATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatc                                                       <  2:6982870/93‑1 (MQ=255)
                                     tatGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatcat                                                     >  2:1705806/1‑93 (MQ=255)
                                     tatGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatcat                                                     >  2:8048686/1‑93 (MQ=255)
                                      atGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatcatc                                                    >  2:3087798/1‑93 (MQ=255)
                                      atGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatcatc                                                    <  2:2225431/93‑1 (MQ=255)
                                       tGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATAtcatca                                                   >  2:6078220/1‑93 (MQ=255)
                                       tGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATAtcatca                                                   >  2:6222116/1‑93 (MQ=255)
                                        gACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCaa                                                  >  1:2215610/1‑93 (MQ=255)
                                        gACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCaa                                                  <  2:8413673/93‑1 (MQ=255)
                                        gACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCaa                                                  <  1:5655085/93‑1 (MQ=255)
                                        gACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCaa                                                  <  1:8422651/93‑1 (MQ=255)
                                         aCGCACTGTGGTGCGGCTGATGGTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCGAt                                                 >  1:2915905/1‑93 (MQ=255)
                                           gCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGa                                               <  2:3216493/93‑1 (MQ=255)
                                           gCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGa                                               >  1:3035265/1‑93 (MQ=255)
                                            cACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGAc                                              >  1:2428596/1‑93 (MQ=255)
                                            cACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGAc                                              <  2:4389107/93‑1 (MQ=255)
                                            cACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGAc                                              >  1:8550511/1‑93 (MQ=255)
                                             aCTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGAcc                                             >  1:5031570/1‑93 (MQ=255)
                                              cTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCt                                            <  2:3457134/93‑1 (MQ=255)
                                              cTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCt                                            >  2:8533272/1‑93 (MQ=255)
                                               tgtgGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTg                                           <  2:1383499/93‑1 (MQ=255)
                                               tgtgGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTg                                           >  1:4958429/1‑93 (MQ=255)
                                                gtggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTgg                                          >  2:7566868/1‑93 (MQ=255)
                                                 tggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGt                                         <  1:4451428/93‑1 (MQ=255)
                                                 tggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGt                                         >  1:6650351/1‑93 (MQ=255)
                                                 tggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGt                                         <  1:6072196/93‑1 (MQ=255)
                                                 ggctgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGt                                         <  2:8668273/90‑1 (MQ=255)
                                                   gtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTa                                       >  1:4835611/1‑93 (MQ=255)
                                                      gggCTGATGTTGTACATTTCCGCCAGCGCACTTTGAGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATc                                    <  2:7206541/92‑1 (MQ=255)
                                                      cGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATc                                    <  2:3982904/93‑1 (MQ=255)
                                                       ggCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATca                                   >  2:2106606/1‑93 (MQ=255)
                                                       ggCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATca                                   <  1:7532954/93‑1 (MQ=255)
                                                         cTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATcaca                                 <  2:1928311/93‑1 (MQ=255)
                                                         cTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATcaca                                 <  1:7264113/93‑1 (MQ=255)
                                                          ttaTGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAt                                <  1:6366792/91‑1 (MQ=255)
                                                          tGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAt                                <  1:2543582/93‑1 (MQ=255)
                                                           gATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAtt                               >  2:5406088/1‑93 (MQ=255)
                                                            aTGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATtg                              >  1:3695990/1‑93 (MQ=255)
                                                             tgttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATtgt                             <  2:6277159/93‑1 (MQ=255)
                                                              gttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATtgtg                            >  2:7086573/1‑93 (MQ=255)
                                                              gttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATtgtg                            >  1:3964075/1‑93 (MQ=255)
                                                               ttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTgg                           >  1:5981354/1‑93 (MQ=255)
                                                               ttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTgg                           <  1:1749331/93‑1 (MQ=255)
                                                               ttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTgg                           >  2:792949/1‑93 (MQ=255)
                                                               ttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTgg                           >  2:5904853/1‑93 (MQ=255)
                                                                tgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGc                          <  1:1594713/93‑1 (MQ=255)
                                                                 gtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCg                         >  2:3787999/1‑93 (MQ=255)
                                                                  tACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGt                        <  2:4115428/93‑1 (MQ=255)
                                                                  tACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGt                        >  1:1236767/1‑93 (MQ=255)
                                                                     ctttCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGGCCCATATCATCAATGCCCTGGTTAATCACATTGTGGCGTaaa                     <  1:2860915/92‑1 (MQ=255)
                                                                     aTTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTaaa                     <  2:6815851/93‑1 (MQ=255)
                                                                     aTTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTaaa                     >  1:2720489/1‑93 (MQ=255)
                                                                     aTTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTaaa                     <  2:3322207/93‑1 (MQ=255)
                                                                          cGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTtg                <  1:687181/93‑1 (MQ=255)
                                                                            ccAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTtgtg              >  2:4932911/1‑93 (MQ=255)
                                                                             cAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGa             >  1:2081901/1‑93 (MQ=255)
                                                                             cAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGa             >  1:4538630/1‑93 (MQ=255)
                                                                              aGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGaa            <  1:7722409/93‑1 (MQ=255)
                                                                               gcgcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAAc           <  1:6267602/93‑1 (MQ=255)
                                                                               gcgcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAAc           <  1:3788547/93‑1 (MQ=255)
                                                                                 gcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGa         >  1:2958906/1‑93 (MQ=255)
                                                                                 gcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGa         <  1:4186498/93‑1 (MQ=255)
                                                                                  cACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGAc        >  1:524555/1‑93 (MQ=255)
                                                                                   aCTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACt       >  1:7078949/1‑93 (MQ=255)
                                                                                   aCTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACt       <  1:7036759/93‑1 (MQ=255)
                                                                                    cTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTc      >  2:669432/1‑93 (MQ=255)
                                                                                     tttGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCa     <  1:7383357/93‑1 (MQ=255)
                                                                                     tttGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCa     >  2:775349/1‑93 (MQ=255)
                                                                                     tttGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCa     <  1:7911591/93‑1 (MQ=255)
                                                                                      ttGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCat    <  1:7136352/93‑1 (MQ=255)
                                                                                       tGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCata   <  2:111379/93‑1 (MQ=255)
                                                                                             |                                                                                       
TGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATAT  >  NC_000913/4593243‑4593422

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: