Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,923,014 A→G I83T (ATA→ACA)  gidA ← glucose‑inhibited division

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,923,0140AG85.0% 711.7 / 114.1 306I83T (ATA→ACA) gidAglucose‑inhibited division
Reads supporting (aligned to +/- strand):  ref base A (21/25);  new base G (138/122);  total (159/147)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 307 total aligned reads displayed.

TCCAGCGCCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCG  >  NC_000913/3922924‑3923106
                                                                                          |                                                                                            
tcCAGCGCCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTc                                                                                            <  1:6139294/93‑1 (MQ=255)
tcCAGCGCCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTc                                                                                            <  2:3373593/93‑1 (MQ=255)
tcCAGCGCCGTACGCACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTc                                                                                            <  1:3384352/93‑1 (MQ=255)
 ccAGCGCCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTcc                                                                                           >  1:5903737/1‑93 (MQ=255)
   aGCGCCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTa                                                                                         <  2:3848557/93‑1 (MQ=255)
    gcgcCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTaa                                                                                        <  2:6209046/93‑1 (MQ=255)
    gcgcCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTaa                                                                                        <  1:6444848/93‑1 (MQ=255)
    gcgcCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTATCCTaa                                                                                        >  1:3910730/1‑93 (MQ=255)
      gcCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAAc                                                                                      <  2:32224/93‑1 (MQ=255)
      gcCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAAc                                                                                      >  2:1258079/1‑93 (MQ=255)
      gcCGCACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAAc                                                                                      >  2:7456765/1‑93 (MQ=255)
        cgtacgtaCCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTg                                                                                    <  2:731874/93‑1 (MQ=255)
         gtacgtacCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTgg                                                                                   >  2:1061608/1‑93 (MQ=255)
          tacgtacCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGa                                                                                  >  2:328711/1‑93 (MQ=255)
          tacgtacCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGa                                                                                  <  1:5047293/93‑1 (MQ=255)
            cgtacCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATa                                                                                >  2:5005827/1‑93 (MQ=255)
            cgtacCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATa                                                                                <  2:5989050/93‑1 (MQ=255)
             gtacCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATAc                                                                               <  2:4769964/93‑1 (MQ=255)
               acCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATAccc                                                                             >  2:302173/1‑93 (MQ=255)
                ccGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCg                                                                            <  1:4744399/93‑1 (MQ=255)
                 cGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGc                                                                           >  1:2291312/1‑93 (MQ=255)
                  gCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGcc                                                                          >  1:839188/1‑93 (MQ=255)
                  gCCGGACGGTAGAGCACCCGATCCGCCTGAGCTCGGGTAGCGCGACCCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATCCCCGcc                                                                          <  1:1281846/93‑1 (MQ=255)
                   ccTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCt                                                                         <  2:4369644/93‑1 (MQ=255)
                    cTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTg                                                                        >  1:1934757/1‑93 (MQ=255)
                    cTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTg                                                                        >  1:2590608/1‑93 (MQ=255)
                     tGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTga                                                                       <  1:7054037/93‑1 (MQ=255)
                     tGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTga                                                                       <  1:2051710/93‑1 (MQ=255)
                     tGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTga                                                                       >  2:1991696/1‑93 (MQ=255)
                       aCGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTgatc                                                                     >  2:6098463/1‑93 (MQ=255)
                       aCGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTgatc                                                                     <  1:4048441/93‑1 (MQ=255)
                         ggTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTgatcga                                                                   >  2:3729586/1‑93 (MQ=255)
                          gTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTgatcgat                                                                  <  2:1714252/93‑1 (MQ=255)
                                 acacGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTc                                                           >  1:6671678/1‑93 (MQ=255)
                                  cacGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCg                                                          >  2:2956197/1‑93 (MQ=255)
                                    cGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGcc                                                        >  2:5475369/1‑93 (MQ=255)
                                    cGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGcc                                                        <  1:5809101/93‑1 (MQ=255)
                                    cGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGcc                                                        >  1:1854088/1‑93 (MQ=255)
                                       tCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATc                                                     >  1:3569705/1‑93 (MQ=255)
                                       tCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATc                                                     <  1:110499/93‑1 (MQ=255)
                                       tCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATc                                                     >  1:1352565/1‑93 (MQ=255)
                                       tCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATc                                                     >  2:7272636/1‑93 (MQ=255)
                                        ccGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCa                                                    <  1:6264798/93‑1 (MQ=255)
                                          gCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGa                                                  >  1:5151803/1‑93 (MQ=255)
                                          gCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGa                                                  >  2:3611528/1‑93 (MQ=255)
                                            cTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGAcc                                                >  2:6669219/1‑93 (MQ=255)
                                            cTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGAcc                                                >  1:6815781/1‑93 (MQ=255)
                                              gAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGAccgc                                              >  1:3667260/1‑93 (MQ=255)
                                               aGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGAccgcc                                             <  1:6513119/93‑1 (MQ=255)
                                               aGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGAccgcc                                             <  1:5604113/93‑1 (MQ=255)
                                                gCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGAccgccg                                            >  2:2920369/1‑93 (MQ=255)
                                                 cTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGa                                           <  1:3072227/93‑1 (MQ=255)
                                                  tCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAg                                          >  2:7358271/1‑93 (MQ=255)
                                                  tCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAg                                          >  2:422443/1‑93 (MQ=255)
                                                    gggTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTg                                        >  1:6040023/1‑93 (MQ=255)
                                                      ggaGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCa                                      >  1:1161829/3‑93 (MQ=255)
                                                        aGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATc                                    <  1:3768591/93‑1 (MQ=255)
                                                         gcgcGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATcc                                   >  2:4647898/1‑93 (MQ=255)
                                                         gcgcGAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATcc                                   <  2:2716280/93‑1 (MQ=255)
                                                          cgcgAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCa                                  <  2:6350593/93‑1 (MQ=255)
                                                           gcgAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCAc                                 <  2:4574225/93‑1 (MQ=255)
                                                           gcgAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCAc                                 >  2:7467011/1‑93 (MQ=255)
                                                            cgAACCGCCGGTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCAct                                <  2:5869684/93‑1 (MQ=255)
                                                             gAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCActt                               >  1:6885319/1‑93 (MQ=255)
                                                             gAACCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCActt                               >  2:4816080/1‑93 (MQ=255)
                                                               aCCGCCGGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCActtct                             >  2:3887421/1‑93 (MQ=255)
                                                                 cgccgGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCttt                           >  1:4527420/1‑93 (MQ=255)
                                                                 cgccgGTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCttt                           >  2:5627494/1‑93 (MQ=255)
                                                                  gccgGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCtttt                          >  1:945166/1‑93 (MQ=255)
                                                                  gccgGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCtttt                          >  1:2795694/1‑93 (MQ=255)
                                                                   ccgGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTa                         >  1:1608263/1‑93 (MQ=255)
                                                                    cgGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATCCCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTAc                        <  2:4010105/93‑1 (MQ=255)
                                                                    cgGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTAc                        >  1:1443097/1‑93 (MQ=255)
                                                                    cgGTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTAc                        <  1:7186294/93‑1 (MQ=255)
                                                                      gTCCTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCa                      >  2:5317405/1‑93 (MQ=255)
                                                                      gTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCa                      >  1:1770579/1‑93 (MQ=255)
                                                                      gTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCa                      >  2:708304/1‑93 (MQ=255)
                                                                        ccTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGa                    >  1:1326511/1‑93 (MQ=255)
                                                                        ccTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGa                    <  1:4846741/93‑1 (MQ=255)
                                                                        ccTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGa                    <  2:7189907/93‑1 (MQ=255)
                                                                        ccTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGCGCATCCACTTCTTTTACCAGa                    <  1:2325334/93‑1 (MQ=255)
                                                                         cTTTGCTTGCGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGAt                   >  1:6594425/1‑93 (MQ=255)
                                                                           ttgcttgcGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGt                 >  1:2200218/1‑93 (MQ=255)
                                                                           ttgcttgcGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGt                 >  1:1732644/1‑93 (MQ=255)
                                                                            tgcttgcGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTc                >  2:4324958/1‑93 (MQ=255)
                                                                            tgcttgcGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTc                <  2:2521728/93‑1 (MQ=255)
                                                                             gcttgcGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTcc               >  1:5230921/1‑93 (MQ=255)
                                                                             gcttgcGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTcc               <  1:5088595/93‑1 (MQ=255)
                                                                              cttgcGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTccc              <  2:4422939/93‑1 (MQ=255)
                                                                               ttgcGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCt             <  2:5283261/93‑1 (MQ=255)
                                                                               ttgcGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCt             >  1:7432690/1‑93 (MQ=255)
                                                                               ttgcGTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCt             <  2:2666865/93‑1 (MQ=255)
                                                                                tgcGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCtt            <  2:6042225/93‑1 (MQ=255)
                                                                                  gttttGGTGTCCTAAACGGGTTTCCCGCCTCATCGATCGATTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTcc          <  1:291659/91‑1 (MQ=255)
                                                                                   gTTTAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTccc         <  1:3264036/93‑1 (MQ=255)
                                                                                   gTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTccc         <  2:1401940/93‑1 (MQ=255)
                                                                                    tttAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCa        >  1:297124/1‑93 (MQ=255)
                                                                                      tAGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAAt      <  2:2320130/93‑1 (MQ=255)
                                                                                       aGTGTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATa     >  1:6377590/1‑93 (MQ=255)
                                                                                          aTCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATAccg  >  1:4376944/1‑93 (MQ=255)
                                                                                          |                                                                                            
TCCAGCGCCGTACGTACCGCCTGACGGTAGAGCACACGATCCGCCTGAGCTCGGGTAGCGCGAACCGCCGGTCCTTTGCTTGCGTTTAGTATCCTAAACTGGATACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCG  >  NC_000913/3922924‑3923106

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