Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,588,821 A→G I15T (ATA→ACA)  yjiY ← putative carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,588,8210AG96.6% 676.5 / 10.0 204I15T (ATA→ACA) yjiYputative carbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (1/6);  new base G (76/121);  total (77/127)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

Only 100 of 224 total aligned reads displayed.

ACGATCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTATGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGTTTGCGTTGCGTTTG  >  NC_000913/4588731‑4588913
                                                                                          |                                                                                            
ccGATCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGtgtg                                                                                            <  2:2957782/92‑1 (MQ=255)
aCGATCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGtgtg                                                                                            >  1:2697834/1‑93 (MQ=255)
 cGATCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGgtgtgt                                                                                           <  1:5175944/93‑1 (MQ=255)
 cGATCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGgtgtgt                                                                                           <  2:1100644/93‑1 (MQ=255)
 cGATCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGgtgtgt                                                                                           <  2:5003470/93‑1 (MQ=255)
 cGATCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGgtgtgt                                                                                           >  2:2558436/1‑93 (MQ=255)
  gATCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGtgtgtg                                                                                          >  2:2520146/1‑93 (MQ=255)
   aTCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGc                                                                                         >  1:1836093/1‑93 (MQ=255)
   aTCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGc                                                                                         >  2:2361578/1‑93 (MQ=255)
     ccACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCtt                                                                                       >  1:82650/1‑93 (MQ=255)
      cacaGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTg                                                                                      <  2:2822254/93‑1 (MQ=255)
        caGGGCTCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGaa                                                                                    <  2:5047140/93‑1 (MQ=255)
          gggCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAata                                                                                  >  1:1557689/1‑93 (MQ=255)
          gggCGCTGATGTGCTCCCCCCGACATAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAata                                                                                  <  1:5098332/93‑1 (MQ=255)
            gcgcTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATc                                                                                <  1:1668149/93‑1 (MQ=255)
             ggcggATGCGCGCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCt                                                                               <  1:5992707/89‑1 (MQ=255)
              gcTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCtt                                                                              <  2:5964814/93‑1 (MQ=255)
                tcaTGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATGCCGAGAATCACCCAGGGTGTGTGCTTGAATATCtttt                                                                            <  1:766490/91‑1 (MQ=255)
                 gATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCttttt                                                                           <  2:6006140/93‑1 (MQ=255)
                 gATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCttttt                                                                           <  2:3300477/93‑1 (MQ=255)
                  aTGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTa                                                                          <  1:32860/93‑1 (MQ=255)
                  aTGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTa                                                                          <  1:702214/93‑1 (MQ=255)
                  aTGTGCGCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTa                                                                          <  1:2808560/93‑1 (MQ=255)
                   tgtgCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAg                                                                         >  1:2745742/1‑93 (MQ=255)
                    gtgCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGt                                                                        <  2:1688684/93‑1 (MQ=255)
                    gtgCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGt                                                                        <  1:201165/93‑1 (MQ=255)
                     tgCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTa                                                                       >  2:2305942/1‑93 (MQ=255)
                     gctcccccccGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTa                                                                       <  2:3468598/89‑1 (MQ=255)
                     gccTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTa                                                                       <  1:1711531/91‑1 (MQ=255)
                      gCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTAt                                                                      <  1:5279672/93‑1 (MQ=255)
                         ccccccGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCa                                                                   <  1:957098/93‑1 (MQ=255)
                          cccccGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCAt                                                                  <  2:4131299/93‑1 (MQ=255)
                          cccccGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCCAGAATCACCCAGGGTGTGTGCTTGCATATCTTTTTAGTATGCAt                                                                  >  1:314502/1‑93 (MQ=255)
                             ccGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGt                                                               <  1:3305345/93‑1 (MQ=255)
                              cGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTa                                                              >  2:1308221/1‑93 (MQ=255)
                              cGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTa                                                              <  1:4044962/93‑1 (MQ=255)
                                aCGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTaaa                                                            >  1:1450208/1‑93 (MQ=255)
                                aCGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTaaa                                                            >  1:5468717/1‑93 (MQ=255)
                                  gTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTATGTGCTTGAATATCTTTTTAGTATCCATAGTAAAAc                                                          <  2:4441286/93‑1 (MQ=255)
                                   tAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAAcc                                                         <  1:2097251/93‑1 (MQ=255)
                                   tAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAAcc                                                         <  2:2024993/93‑1 (MQ=255)
                                   tAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAAcc                                                         >  1:4529232/1‑93 (MQ=255)
                                    aaTGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCt                                                        <  2:1406286/93‑1 (MQ=255)
                                     aTGCAACTCCCGCGAGGCAGAATGCACCGATTATTCCGAGAATCCCCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTg                                                       <  2:1498582/93‑1 (MQ=255)
                                     aTGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTg                                                       >  2:3795108/1‑93 (MQ=255)
                                     aTGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTg                                                       <  2:4038239/93‑1 (MQ=255)
                                     aTGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTg                                                       <  2:5840021/93‑1 (MQ=255)
                                      tGCAACTACCGCGAGACAGAATGCACCGATGATTCCGAGAATCACCCAGGGGGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTgg                                                      >  2:3678447/1‑93 (MQ=255)
                                       gCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGc                                                     >  2:311180/1‑93 (MQ=255)
                                         aaCTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCAt                                                   >  1:1209075/1‑93 (MQ=255)
                                           cTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGt                                                 >  1:4284240/1‑93 (MQ=255)
                                           cTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAAACTGGCATGt                                                 >  1:5447298/1‑93 (MQ=255)
                                            tACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTa                                                <  2:2802768/93‑1 (MQ=255)
                                              ccGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTAtt                                              >  1:2259183/1‑93 (MQ=255)
                                              ccGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTAtt                                              >  2:5895370/1‑93 (MQ=255)
                                                   aGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTa                                         <  1:5648414/93‑1 (MQ=255)
                                                   aGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTa                                         >  1:4247981/1‑93 (MQ=255)
                                                      gaGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATa                                      <  1:2478223/92‑1 (MQ=255)
                                                       aGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAg                                     <  1:699234/93‑1 (MQ=255)
                                                       aGAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAg                                     <  2:3375082/93‑1 (MQ=255)
                                                        gAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGt                                    <  2:56291/93‑1 (MQ=255)
                                                        gAATGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGt                                    >  1:3966038/1‑93 (MQ=255)
                                                         aaTGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGtt                                   <  2:348304/93‑1 (MQ=255)
                                                         aaTGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGtt                                   <  2:678213/93‑1 (MQ=255)
                                                         aaTGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGtt                                   <  2:732778/93‑1 (MQ=255)
                                                         aaTGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGtt                                   <  2:746967/93‑1 (MQ=255)
                                                          aTGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTg                                  <  1:5820091/93‑1 (MQ=255)
                                                          aTGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTg                                  <  2:1082961/93‑1 (MQ=255)
                                                          aTGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTg                                  <  2:3134305/93‑1 (MQ=255)
                                                           tGCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTgg                                 >  2:5165329/1‑93 (MQ=255)
                                                            gCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGc                                >  2:1366195/1‑93 (MQ=255)
                                                            gCACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGc                                <  2:4862902/93‑1 (MQ=255)
                                                            gCACCGATGATTCCGAGAAGTCCACAGGGTGTGTTCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGc                                <  2:2138104/93‑1 (MQ=255)
                                                             cACCGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGcc                               <  1:3209166/93‑1 (MQ=255)
                                                              acggatgatTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCg                              >  1:3862357/4‑93 (MQ=255)
                                                               ccGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGa                             >  1:5730448/1‑93 (MQ=255)
                                                                cGATGATTCCGAGAATCACCCAGGGTGTGTGCTTTAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTAGCCGaa                            >  1:3504780/1‑93 (MQ=255)
                                                                cGATGATTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGaa                            <  2:5571281/93‑1 (MQ=255)
                                                                cGATGATTCCGAGAATCACCCAGGGTGGGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGaa                            <  2:1721680/93‑1 (MQ=255)
                                                                  atgatTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGc                          >  1:3070236/1‑93 (MQ=255)
                                                                  atgatTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGc                          <  2:4610520/93‑1 (MQ=255)
                                                                  atgatTCCAAGAATAACCCAGGGTGTTTGCTTAAATATCTTATTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGc                          <  1:24581/93‑1 (MQ=255)
                                                                   tgatTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGcc                         >  2:4741484/1‑93 (MQ=255)
                                                                     atTCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGt                       >  1:3939168/1‑93 (MQ=255)
                                                                      ttCCGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGtt                      >  2:526742/1‑93 (MQ=255)
                                                                      atCCCAGAATCACCCAGGGTTTTTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAATCCGtt                      <  1:2242632/92‑1 (MQ=255)
                                                                         cGAGAATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTa                   >  2:1262996/1‑93 (MQ=255)
                                                                          gagaATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAg                  >  2:822808/1‑93 (MQ=255)
                                                                          gagaATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAg                  <  1:738288/93‑1 (MQ=255)
                                                                          gagaATCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAg                  >  2:1482279/1‑93 (MQ=255)
                                                                             aaTCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGtt               <  2:5357138/93‑1 (MQ=255)
                                                                             aaTCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGtt               <  2:3889941/93‑1 (MQ=255)
                                                                               tCACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGTttg             <  2:323568/93‑1 (MQ=255)
                                                                                cACCCAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGTttgc            <  2:3523653/93‑1 (MQ=255)
                                                                                   ccAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGTttgcgtt         <  2:2698123/93‑1 (MQ=255)
                                                                                   ccAGGGTATGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGTttgcgtt         <  1:3890475/93‑1 (MQ=255)
                                                                                    cAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGTttgcgttg        >  2:510920/1‑93 (MQ=255)
                                                                                    cAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGTttgcgttg        >  1:2137455/1‑93 (MQ=255)
                                                                                    cAGGGTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGTttgcgttg        >  1:5036996/1‑93 (MQ=255)
                                                                                      gggTGTGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGTttgcgttgcg      >  1:3532025/1‑93 (MQ=255)
                                                                                          |                                                                                            
ACGATCCACAGGGCGCTGATGTGCTCCCCCCGACGTAATGCAACTACCGCGAGGCAGAATGCACCGATGATTCCGAGAATCACCCAGGGTATGTGCTTGAATATCTTTTTAGTATCCATAGTAAAACCTGGCATGTATTGATTAATAGTTGGCCGAAGCCGTTCTAGGTTTGCGTTGCGTTTG  >  NC_000913/4588731‑4588913

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