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.5% 668.5 / ‑3.9 203intergenic (‑22/‑209)yjjM/yjjNhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (113/89);  total (114/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 225 total aligned reads displayed.

TGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCATAT  >  NC_000913/4593243‑4593422
                                                                                              |                                                                                       
tGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGAt                                                                                          >  2:2384640/1‑93 (MQ=255)
 gCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATg                                                                                        <  2:1949699/93‑1 (MQ=255)
 gCGTCAAGACGCCGCATTCGCGTAAGAGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTCCATTTCCGCCAGCGCACTTTGCGATg                                                                                        >  2:4080871/1‑93 (MQ=255)
    tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa                                                                                     <  1:504959/93‑4 (MQ=255)
    tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa                                                                                     >  1:4398539/1‑90 (MQ=255)
          cgccgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTgggg                                                                                <  2:1596429/90‑1 (MQ=255)
            ccgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGaa                                                                              >  1:390149/2‑93 (MQ=255)
             cgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAg                                                                             <  2:191609/93‑1 (MQ=255)
             cgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAg                                                                             <  2:204158/93‑1 (MQ=255)
              gcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAgg                                                                            <  2:611271/93‑1 (MQ=255)
               cATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAggg                                                                           >  1:804757/1‑93 (MQ=255)
                agtCGCGTAAGTGGCTGAGAATATGACGCACTGGGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGa                                                                          <  2:2442481/91‑1 (MQ=255)
                aTTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGa                                                                          <  1:143549/93‑1 (MQ=255)
                aTTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGa                                                                          >  2:43322/1‑93 (MQ=255)
                aTTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGa                                                                          >  1:4744397/1‑93 (MQ=255)
                 ttCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGat                                                                         >  2:2828035/1‑93 (MQ=255)
                 ttCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGat                                                                         >  1:2803790/1‑93 (MQ=255)
                 ttCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGat                                                                         <  2:697760/93‑1 (MQ=255)
                  tCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGata                                                                        <  1:2594233/93‑1 (MQ=255)
                     cgTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGa                                                                     >  2:3242556/1‑93 (MQ=255)
                     cgTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGa                                                                     <  1:1880583/93‑1 (MQ=255)
                     cgTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGa                                                                     >  1:3991649/1‑93 (MQ=255)
                     cgTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGa                                                                     >  1:2859594/1‑93 (MQ=255)
                      gTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCGCTTTGCGATGGAAGTGGGGAAGGGATATGAc                                                                    >  1:3059765/1‑93 (MQ=255)
                      gTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAc                                                                    >  2:3047897/1‑93 (MQ=255)
                      gTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAc                                                                    >  2:3700353/1‑93 (MQ=255)
                       tAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAcc                                                                   >  1:3708502/1‑93 (MQ=255)
                       tAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAcc                                                                   >  2:3951511/1‑93 (MQ=255)
                        aaGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAccc                                                                  >  1:348563/1‑93 (MQ=255)
                        aaGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAccc                                                                  <  1:2670035/93‑1 (MQ=255)
                         aGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGAcccc                                                                 <  1:1141756/93‑1 (MQ=255)
                          gTGGCTGAGAATATGACGCACTGGGGCGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCcg                                                                <  2:2772011/93‑1 (MQ=255)
                            ggCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCcgcg                                                              <  1:2064134/93‑1 (MQ=255)
                              cTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGcc                                                            >  2:1236300/1‑93 (MQ=255)
                               tGAGAATATGACGCACTGCGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCa                                                           >  1:4555721/1‑93 (MQ=255)
                                gagaATATGACTCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCat                                                          >  1:152972/1‑93 (MQ=255)
                                gagaATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCat                                                          >  2:2738550/1‑93 (MQ=255)
                                gagaATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCat                                                          <  2:2999355/93‑1 (MQ=255)
                                 agaATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCata                                                         >  1:3518199/1‑93 (MQ=255)
                                   aaTATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatc                                                       <  2:919034/93‑1 (MQ=255)
                                   aaTATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatc                                                       <  2:4130101/93‑1 (MQ=255)
                                     tatGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATatcat                                                     >  1:4657886/1‑93 (MQ=255)
                                       tGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATAtcatca                                                   <  1:3357503/93‑1 (MQ=255)
                                        gACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCaa                                                  <  2:3301476/93‑1 (MQ=255)
                                         aCGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAAt                                                 >  2:3164426/1‑93 (MQ=255)
                                           gCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATTGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGa                                               >  1:85483/1‑93 (MQ=255)
                                           gCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGa                                               >  2:1436583/1‑93 (MQ=255)
                                           gCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGa                                               >  2:1074896/1‑93 (MQ=255)
                                             agtgtgGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGAcc                                             >  1:4307919/3‑93 (MQ=255)
                                             aCTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGAcc                                             <  1:3817055/93‑1 (MQ=255)
                                              cTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCt                                            >  2:1723750/1‑93 (MQ=255)
                                               tgtgGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTg                                           >  1:1335039/1‑93 (MQ=255)
                                               tgtgGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTg                                           >  2:1136234/1‑93 (MQ=255)
                                                gtggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAATAGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTgg                                          >  2:2089634/1‑93 (MQ=255)
                                                gtggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTgg                                          <  1:2308995/93‑1 (MQ=255)
                                                  ggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGtt                                        <  2:1302759/93‑1 (MQ=255)
                                                  ggtgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGtt                                        <  1:537995/93‑1 (MQ=255)
                                                  ggtgCGGCTGATGTGGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGtt                                        <  2:1387299/93‑1 (MQ=255)
                                                    tgCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTaa                                      <  2:1843193/93‑1 (MQ=255)
                                                     gCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAAt                                     >  2:4577873/1‑93 (MQ=255)
                                                      cGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATc                                    <  1:1934453/93‑1 (MQ=255)
                                                      cGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATc                                    <  1:4414369/93‑1 (MQ=255)
                                                       ggCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGAAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATca                                   <  2:4699512/93‑1 (MQ=255)
                                                         cTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATcaca                                 <  1:2337041/93‑1 (MQ=255)
                                                         cTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATcaca                                 >  2:1781584/1‑93 (MQ=255)
                                                          tGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAt                                >  2:3326533/1‑93 (MQ=255)
                                                          tGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAt                                >  1:2440139/1‑93 (MQ=255)
                                                           gATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAtt                               >  2:3444968/1‑93 (MQ=255)
                                                           gATGTTGTACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACAtt                               <  1:634874/93‑1 (MQ=255)
                                                             tgttgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATtgt                             <  1:278587/93‑1 (MQ=255)
                                                                tgtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGc                          <  2:901739/93‑1 (MQ=255)
                                                                 gtACATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCg                         <  2:3780668/93‑1 (MQ=255)
                                                                   aCATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTa                       >  1:1803899/1‑93 (MQ=255)
                                                                    cATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTaa                      >  2:3503994/1‑93 (MQ=255)
                                                                    cATTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTCGCGTaa                      >  2:4180834/1‑93 (MQ=255)
                                                                     aTTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTaaa                     <  2:2341432/93‑1 (MQ=255)
                                                                     aTTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTaaa                     <  1:1662066/93‑1 (MQ=255)
                                                                     aTTTCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTaaa                     <  2:4553482/93‑1 (MQ=255)
                                                                      tttCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAt                    >  1:4189386/1‑93 (MQ=255)
                                                                        tCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAttt                  <  1:2925051/93‑1 (MQ=255)
                                                                        tCCGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAttt                  >  1:969974/1‑93 (MQ=255)
                                                                         ccGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAtttt                 >  1:2243404/1‑93 (MQ=255)
                                                                         ccGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAAtttt                 >  1:1018021/1‑93 (MQ=255)
                                                                          cGCCAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTtg                <  1:4426537/93‑1 (MQ=255)
                                                                             cAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGa             >  1:4504694/1‑93 (MQ=255)
                                                                             cAGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGa             <  2:4665793/93‑1 (MQ=255)
                                                                              aGCGCACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGaa            <  1:330058/93‑1 (MQ=255)
                                                                               gcgcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGCCCTGGTTAATCACATTGTGGCGTAAATTTTGTGAAc           >  2:1003403/1‑93 (MQ=255)
                                                                               gcgcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCAGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAAc           <  1:908334/93‑1 (MQ=255)
                                                                                cgcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACg          >  1:1723560/1‑93 (MQ=255)
                                                                                 gcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGa         >  1:2531184/1‑93 (MQ=255)
                                                                                 gcACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGa         <  1:84655/93‑1 (MQ=255)
                                                                                  cACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGAc        >  1:1073124/1‑93 (MQ=255)
                                                                                  cACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGAc        >  2:3978129/1‑93 (MQ=255)
                                                                                  cACTTTGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGAc        <  1:3097372/93‑1 (MQ=255)
                                                                                     tttGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCa     >  2:313476/1‑93 (MQ=255)
                                                                                     tttGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCa     >  2:3139338/1‑93 (MQ=255)
                                                                                       tGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCata   <  1:4323420/93‑1 (MQ=255)
                                                                                       tGCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCata   <  1:1249751/93‑1 (MQ=255)
                                                                                        gCGATGGAAGTGGGGAAGGGATATGACCCCGCGCCATATCATCAATGACCTGGTTAATCACATTGTGGCGTAAATTTTGTGAACGACTCatat  <  1:40580/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 < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: