Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,424,468 A→C 15.2% intergenic (+82/+48) ulaF → / ← yjfY L‑ribulose 5‑phosphate 4‑epimerase/YhcN family protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,424,4680AC15.2% 470.6 / 81.2 189intergenic (+82/+48)ulaF/yjfYL‑ribulose 5‑phosphate 4‑epimerase/YhcN family protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (88/72);  new base C (14/15);  total (102/87)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

Only 100 of 189 total aligned reads displayed.

TTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTT  >  NC_000913/4424371‑4424565
                                                                                                  |                                                                                                 
ttaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAg                                                                                                 >  1:556559/1‑100 (MQ=255)
 taCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAgg                                                                                                <  1:351670/100‑1 (MQ=255)
 cgatctcAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAgg                                                                                                <  2:4058762/94‑1 (MQ=255)
    gggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAggggn                                                                                             <  2:4122563/100‑2 (MQ=255)
    gggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAggggg                                                                                             >  2:934403/1‑100 (MQ=255)
    gggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAggggg                                                                                             >  2:1421867/1‑100 (MQ=255)
     ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAggggga                                                                                            >  1:368917/1‑99 (MQ=255)
     ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGgc                                                                                            >  1:2464296/1‑100 (MQ=255)
     ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGgc                                                                                            >  1:3726394/1‑100 (MQ=255)
     ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTACATCGAGAGATAACTGGTCTAATTCCGCAGAAAAAAAAA‑AAAAAAAAGCCCCCTGCGAAggggg                                                                                             <  2:368917/99‑1 (MQ=255)
       cAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGgcaa                                                                                          >  1:929656/1‑100 (MQ=255)
         gAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCTGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGgcaaa                                                                                         <  2:718363/99‑1 (MQ=255)
         gAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGgcaaag                                                                                        <  1:954059/100‑1 (MQ=255)
          aaGTAGAACACGCGCTGCGTAATTTTCCTTCATCGTGAGATAACTGGCCTAATTCCGCAGCCTTTTTTC‑AAAAAAAAGCCCCCTTCGCAGGGGgct                                                                                           <  2:3774882/96‑2 (MQ=255)
          aaGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCAGGGGgcttttt                                                                                       <  2:1086648/100‑6 (MQ=255)
          aaGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCAGGGGgcttttt                                                                                       <  1:3091020/100‑6 (MQ=255)
            gTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCAGGGGgcttttttt                                                                                     <  1:1915413/100‑8 (MQ=255)
            gTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCAGGGGgcttttttt                                                                                     <  2:1555014/100‑8 (MQ=255)
            gTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGgcaacgcaa                                                                                     >  1:3173780/1‑100 (MQ=255)
             tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCAGGGGgctttttttt                                                                                    <  1:2264139/100‑9 (MQ=38)
             tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCAGGGGgctttttttt                                                                                    <  2:3666319/100‑9 (MQ=38)
             tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCAGGGGgctttttttt                                                                                    <  1:2789581/100‑9 (MQ=38)
              aGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCAGGGGgcttttttttg                                                                                   <  1:2187777/100‑10 (MQ=38)
              aGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTTCGCAGGGGgcttttttttg                                                                                   <  1:1413284/100‑10 (MQ=38)
              aGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAAc                                                                                   >  1:23310/1‑100 (MQ=255)
              aGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAAc                                                                                   >  2:2190969/1‑100 (MQ=255)
                aaCACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAACCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTa                                                                                 >  1:1819134/1‑100 (MQ=255)
                 acacGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑GAAAAAAAGCACCCTGAGAAGGGGGCAAAGCAAACTAt                                                                                >  1:833210/1‑100 (MQ=255)
                 acacGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTAt                                                                                >  1:3945265/1‑100 (MQ=255)
                 acacGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGCGGCAAAGCAAACTAt                                                                                >  1:4125047/1‑100 (MQ=255)
                   acGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATg                                                                               <  1:3540698/100‑1 (MQ=255)
                   acGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATgg                                                                              <  2:2612494/100‑1 (MQ=255)
                    cgcgCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGc                                                                             <  1:3962344/100‑1 (MQ=255)
                    cgcgCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGc                                                                             >  1:1263923/1‑100 (MQ=255)
                        cTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGAAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATg                                                                         <  1:3819726/100‑1 (MQ=255)
                         tGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGt                                                                        >  2:1991842/1‑100 (MQ=255)
                          gCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATg                                                                         <  1:1432375/98‑1 (MQ=255)
                          gCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGtt                                                                       >  1:3052989/1‑100 (MQ=255)
                           cGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGttt                                                                      >  1:1047199/1‑100 (MQ=255)
                             gAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                                                                    <  1:2934670/100‑1 (MQ=255)
                                aTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTg                                                                 >  2:1194541/1‑100 (MQ=255)
                                 tttCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTTGTTgg                                                                <  1:427728/100‑1 (MQ=255)
                                 tttCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTgg                                                                >  1:3297998/1‑100 (MQ=255)
                                 tttCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTgg                                                                >  1:3992589/1‑100 (MQ=255)
                                 tttCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCACTGCGAAGGGGGCAAAGCAAACTATGGACATGTTTCGTTgg                                                                >  2:834033/1‑100 (MQ=255)
                                  ttcctttcTCGGGAGTGAACTGGTCGACATCTTCAGCCATTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt                                                               <  1:327169/100‑1 (MQ=255)
                                  ttccttccTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt                                                               <  1:3019781/100‑1 (MQ=255)
                                   tccttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGtt                                                              <  2:3889594/100‑1 (MQ=255)
                                   tccttcctCGGGAGATAACTGGTCTAATTACGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGtt                                                              <  2:2325356/100‑1 (MQ=255)
                                    ccttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTta                                                             <  2:2471967/100‑1 (MQ=255)
                                    ccttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTta                                                             >  1:164403/1‑100 (MQ=255)
                                      ttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGCTTTTC‑AAAAAAAAGCCCCCTGCGAAAGGGGCAAAGCAAACTATGGcgatggttcgatggtt                                                              >  2:2769567/1‑97 (MQ=255)
                                            gggAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑TAATAAATTCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt                                                     >  2:53165/1‑100 (MQ=255)
                                            gggAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt                                                     >  1:1105760/1‑100 (MQ=255)
                                            gggAGAGAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt                                                     <  2:4037607/100‑1 (MQ=255)
                                             ggAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTgg                                                      <  2:1686117/98‑1 (MQ=255)
                                             ggAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTg                                                    >  1:2435028/1‑100 (MQ=255)
                                              gagaTAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGc                                                   <  1:72048/100‑1 (MQ=255)
                                               agaTAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCt                                                  >  1:3750533/1‑100 (MQ=255)
                                                gaTAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTa                                                 >  1:1513822/1‑100 (MQ=255)
                                                  tAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGc                                               >  2:380902/1‑100 (MQ=255)
                                                   aaCTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAATCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCg                                              >  1:1662149/1‑100 (MQ=255)
                                                   aaCTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCg                                              >  1:1292087/1‑100 (MQ=255)
                                                   aaCTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGAGCAAAGCAAACTATGGCAATGTTCGGTTGGTTATACCTGGTGCTAGCg                                              >  2:2046295/1‑100 (MQ=255)
                                                         tCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAat                                        <  1:1458365/100‑1 (MQ=255)
                                                          cTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAata                                       <  2:3953796/100‑1 (MQ=255)
                                                           tAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAatat                                      >  1:782592/1‑100 (MQ=255)
                                                             aTTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATcc                                    >  1:2231374/1‑100 (MQ=255)
                                                             aTTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATcc                                    >  1:3801634/1‑100 (MQ=255)
                                                              ttCCGCCGCCGTTTTTC‑AACAAAAAGCCCCCTGCGAAGGGGGCAAAGCCAACTATGGCAATTGTTCGTTGGTTATACCTGTTGCTATCGATAAATATCcg                                   >  2:1633942/1‑100 (MQ=255)
                                                              ttCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTAGGGCAATGTTTCGTTGGTTATTCCTGGTGCTAGCGATAAATATCcg                                   <  2:1158317/100‑1 (MQ=255)
                                                                ccGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCcgcg                                 >  2:3302831/1‑100 (MQ=255)
                                                                ccGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCcgcg                                 >  2:2533754/1‑100 (MQ=255)
                                                                 cGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCgcgc                                <  1:3191025/100‑1 (MQ=255)
                                                                  gCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCt                               >  2:3083414/1‑100 (MQ=255)
                                                                   cAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCt                               >  1:392045/1‑99 (MQ=255)
                                                                      ccGTTTTTC‑AAAAAAAAGGCCACTGCGAAGCGGGGAAAGCAAGCTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCg                           <  2:4181744/100‑1 (MQ=255)
                                                                           tttC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCat                      <  1:283958/100‑1 (MQ=255)
                                                                            ttC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCata                     >  2:1824495/1‑100 (MQ=255)
                                                                              c‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATAtt                   >  1:453792/1‑100 (MQ=255)
                                                                                aaaaaaaaGGCCCCTGCGAAGGGGGCAAAACAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCg                           >  2:1607041/1‑91 (MQ=255)
                                                                                aaaaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTg                  >  2:3822559/1‑100 (MQ=255)
                                                                                 aaaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGc                 >  2:2406369/1‑100 (MQ=255)
                                                                                 aaaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGc                 <  1:2920896/100‑1 (MQ=255)
                                                                                 aaaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGc                 >  2:1669031/1‑100 (MQ=255)
                                                                                 aaaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGc                 >  2:1234311/1‑100 (MQ=255)
                                                                                   gctttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGCGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCg               >  2:2277234/6‑100 (MQ=255)
                                                                                    ctttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGt              >  1:2533777/5‑100 (MQ=255)
                                                                                    ctttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGt              >  1:1569279/5‑100 (MQ=255)
                                                                                    ctttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGt              >  2:727498/5‑100 (MQ=255)
                                                                                     tttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGtt             >  1:644681/4‑100 (MQ=255)
                                                                                     aaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGtt             <  1:2298480/100‑1 (MQ=255)
                                                                                     aaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGtt             <  1:2699021/100‑1 (MQ=255)
                                                                                     aaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGtt             <  1:3006013/100‑1 (MQ=255)
                                                                                      ttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTa            >  1:2601048/3‑100 (MQ=255)
                                                                                      ttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCACTGGCGTGCATATTGCCGTTa            >  2:621045/3‑100 (MQ=255)
                                                                                        gCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACt          <  1:3852752/100‑1 (MQ=255)
                                                                                         cccccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTc         >  2:3575837/1‑100 (MQ=255)
                                                                                         cccccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTc         >  1:3657856/1‑100 (MQ=255)
                                                                                            ccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTcccc      <  1:753034/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
TTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTT  >  NC_000913/4424371‑4424565

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑