Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,519,062 A→G I176T (ATA→ACA)  yjhF ← putative transport system permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,519,0620AG99.8% 1485.1 / ‑5.5 451I176T (ATA→ACA) yjhFputative transport system permease
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (269/181);  total (269/182)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

Only 100 of 468 total aligned reads displayed.

TTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCAAA  >  NC_000913/4518972‑4519151
                                                                                          |                                                                                         
 tAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCa                                                                                        >  2:8101196/1‑93 (MQ=255)
 tAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCa                                                                                        >  2:693653/1‑93 (MQ=255)
 tAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCa                                                                                        <  1:6631704/93‑1 (MQ=255)
 tAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCa                                                                                        >  1:5285717/1‑93 (MQ=255)
   aaGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGt                                                                                      >  2:1824046/1‑93 (MQ=255)
      gCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTAcc                                                                                   <  1:5132619/93‑1 (MQ=255)
        gTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCaa                                                                                 <  2:8663302/93‑1 (MQ=255)
         tCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAAt                                                                                >  2:9835206/1‑93 (MQ=255)
            ggAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCt                                                                             <  1:1183956/93‑1 (MQ=255)
            ggAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCt                                                                             >  1:84174/1‑93 (MQ=255)
             gAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTa                                                                            >  2:9333393/1‑93 (MQ=255)
             gAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTa                                                                            >  2:4456243/1‑93 (MQ=255)
              aataatGCGCTGGCAGAATTTATTCAATACCGGGCCGTAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTAc                                                                           <  2:6098777/93‑1 (MQ=255)
              aataatGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTAc                                                                           <  1:7981805/93‑1 (MQ=255)
              aataatGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTAc                                                                           <  2:6267347/93‑1 (MQ=255)
               ataatGCGCTGGCAGAATTTATTCAATCCCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACa                                                                          <  1:1709350/93‑1 (MQ=255)
               ataatGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACa                                                                          >  1:9439966/1‑93 (MQ=255)
               ataatGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACa                                                                          >  2:4167771/1‑93 (MQ=255)
               ataatGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACa                                                                          >  1:2288857/1‑93 (MQ=255)
               ataatGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACa                                                                          <  1:2501712/93‑1 (MQ=255)
                aaatGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACAt                                                                         <  2:2776104/92‑1 (MQ=255)
                   tGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAg                                                                      <  2:7352405/93‑1 (MQ=255)
                     cgcTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAgcg                                                                    >  2:1079872/1‑93 (MQ=255)
                      gcTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAgcgc                                                                   >  1:275215/1‑93 (MQ=255)
                      gcTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAgcgc                                                                   >  2:7199721/1‑93 (MQ=255)
                      gcTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAgcgc                                                                   >  2:7906601/1‑93 (MQ=255)
                       cTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAgcgca                                                                  >  1:8491318/1‑92 (MQ=255)
                        tGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCa                                                                 >  2:1255496/1‑93 (MQ=255)
                         ggCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCat                                                                <  1:4666218/93‑1 (MQ=255)
                          gCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCata                                                               <  1:4583612/93‑1 (MQ=255)
                          gCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCata                                                               >  2:7638050/1‑93 (MQ=255)
                          gCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGGTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCata                                                               <  1:6220837/93‑1 (MQ=255)
                           cAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCatat                                                              >  1:7060871/1‑93 (MQ=255)
                           cAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCatat                                                              <  2:4958198/93‑1 (MQ=255)
                           cAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCatat                                                              >  2:7780387/1‑93 (MQ=255)
                             gAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTc                                                            <  1:6440327/93‑1 (MQ=255)
                             gAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTc                                                            <  2:8313186/93‑1 (MQ=255)
                              aaTTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTcc                                                           >  2:223277/1‑93 (MQ=255)
                              aaTTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTcc                                                           <  1:8657535/93‑1 (MQ=255)
                                 ttATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTg                                                        <  1:8864056/93‑1 (MQ=255)
                                 ttATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTg                                                        >  2:8474394/1‑93 (MQ=255)
                                  tATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTgg                                                       <  1:7633439/93‑1 (MQ=255)
                                    ttCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAAACCTACATCAGCGCCATATTCCCTGGCg                                                     <  2:4628326/93‑1 (MQ=255)
                                      cAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGAt                                                   <  1:5833152/93‑1 (MQ=255)
                                        aTACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATaa                                                 >  2:6186727/1‑93 (MQ=255)
                                         tACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAAc                                                >  1:3919560/1‑93 (MQ=255)
                                         tACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAAc                                                >  2:8887355/1‑93 (MQ=255)
                                         tACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAAc                                                <  2:8085650/93‑1 (MQ=255)
                                         aaCCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAAc                                                <  1:8091322/92‑1 (MQ=255)
                                          aCCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACa                                               <  1:1375899/93‑1 (MQ=255)
                                           ccGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAg                                              >  1:3288388/1‑93 (MQ=255)
                                           ccGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAg                                              <  2:1542025/93‑1 (MQ=255)
                                            cGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGt                                             <  2:342827/93‑1 (MQ=255)
                                            cGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGt                                             <  2:951033/93‑1 (MQ=255)
                                              ggCCGCAAAGAATAAAGGTTGGAATGGCGACAAGGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCg                                           <  1:2959161/93‑1 (MQ=255)
                                                ccGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCggg                                         <  2:5745019/93‑1 (MQ=255)
                                                 cGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGc                                        >  1:4870567/1‑93 (MQ=255)
                                                   cAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCg                                      >  2:5893707/1‑93 (MQ=255)
                                                   cAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCg                                      <  2:3050895/93‑1 (MQ=255)
                                                    aaaGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCgg                                     >  2:5150467/1‑93 (MQ=255)
                                                      aGAATAAAGGTTGGAATGCCGACAATGATACCGGATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCgggg                                   <  1:8222634/93‑1 (MQ=255)
                                                       gAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCggggg                                  >  2:3778557/1‑92 (MQ=255)
                                                         aTAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTgc                                >  1:4247295/1‑93 (MQ=255)
                                                         aTAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTgc                                <  1:5915784/93‑1 (MQ=255)
                                                         aTAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGGgc                                >  1:2068/1‑93 (MQ=255)
                                                          tAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTgcg                               >  1:4057118/1‑93 (MQ=255)
                                                          tAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTgcg                               >  2:1022708/1‑93 (MQ=255)
                                                            aaGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTgcggc                             <  2:2692951/93‑1 (MQ=255)
                                                            aaGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCATGGCGATAACAGTCGGGCCGGGGGgcgtc                             >  1:7471620/1‑91 (MQ=255)
                                                            aaGGTTGGAATGCCGACAATGATACCGTATGCCAGCACCAACCCTAAAGCAGCGCCATATTACTCTCCGATAAAAGCCGGGGCGGGGGggggc                             >  1:5148432/1‑93 (MQ=255)
                                                              ggTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGcggcgg                           >  2:6717951/1‑93 (MQ=255)
                                                               gTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCggcggc                          >  2:913598/1‑93 (MQ=255)
                                                                 tGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCaa                        >  1:10020017/1‑93 (MQ=255)
                                                                      tGCCGACAATGATACCGTATGTCAGTACCAATCGTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCaaaaagc                   <  2:7044259/93‑3 (MQ=255)
                                                                      tGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATCACCGTCGGGCCGGGGTGCGGCGGCACAacac                   >  1:7096356/1‑93 (MQ=255)
                                                                      tGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAc                   >  2:2115584/1‑93 (MQ=255)
                                                                      tGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAc                   >  2:7760150/1‑93 (MQ=255)
                                                                       gCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAcc                  >  1:6798649/1‑93 (MQ=255)
                                                                          gACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCgtg               >  1:6116076/1‑93 (MQ=255)
                                                                          gACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCgtg               >  1:1505008/1‑93 (MQ=255)
                                                                          gACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCgtg               <  2:5408758/93‑1 (MQ=255)
                                                                          gACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGCGCGGCGGCAAGAAACCgtg               >  1:2246611/1‑93 (MQ=255)
                                                                           aCAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCgtgt              <  2:3672497/93‑1 (MQ=255)
                                                                           aCAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCgtgt              >  1:9101766/1‑93 (MQ=255)
                                                                            cAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACtgtgtg             <  1:6935762/93‑1 (MQ=255)
                                                                             aaTGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGc            >  1:371356/1‑93 (MQ=255)
                                                                             aaTGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGc            >  2:5657284/1‑93 (MQ=255)
                                                                               tGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGa          >  1:3138363/1‑93 (MQ=255)
                                                                               tGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGa          >  1:5751177/1‑93 (MQ=255)
                                                                               tGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGa          >  1:7813375/1‑93 (MQ=255)
                                                                               tGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGa          >  1:8522625/1‑93 (MQ=255)
                                                                                 aTACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGAcc        <  1:5514475/93‑1 (MQ=255)
                                                                                  tACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCa       >  2:2455209/1‑93 (MQ=255)
                                                                                  tACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCa       <  1:3848890/93‑1 (MQ=255)
                                                                                  tACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCa       >  2:1300532/1‑93 (MQ=255)
                                                                                   aCCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAg      >  1:2450930/1‑93 (MQ=255)
                                                                                   aCCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAg      >  2:4045846/1‑93 (MQ=255)
                                                                                      gTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCaa   >  1:7330729/1‑93 (MQ=255)
                                                                                       tatgTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCaaa  >  1:4172030/1‑93 (MQ=255)
                                                                                       tatgTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCaaa  >  1:9896416/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCAAA  >  NC_000913/4518972‑4519151

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: