Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,424,465 G→T 18.4% intergenic (+79/+51) 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,4650GT18.4% 307.8 / 56.7 125intergenic (+79/+51)ulaF/yjfYL‑ribulose 5‑phosphate 4‑epimerase/YhcN family protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (55/46);  new base T (10/13);  total (66/59)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

Only 100 of 126 total aligned reads displayed.

GCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCG  >  NC_000913/4424366‑4424562
                                                                                                   |                                                                                                 
   tattaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGa                                                                                                >  2:241706/1‑100 (MQ=255)
      taCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAgg                                                                                             >  2:418005/1‑100 (MQ=255)
      taCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAgg                                                                                             >  1:1322219/1‑100 (MQ=255)
       aCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCATATTTTTTCCCCCTGCGAAggg                                                                                            >  1:119457/1‑100 (MQ=255)
         gggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAggggg                                                                                          >  2:1189153/1‑100 (MQ=255)
          ggtaGAAGTAGAACACGCGCTGCGGAAATTTACTTCCTCGGGAGATAACTGGTCTAATTCCGTAGCCGTTTTTTAAAAAAAAGCCCCCTTCGCAgggg                                                                                           <  1:820714/95‑1 (MQ=255)
          ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgc                                                                                         >  2:707929/1‑100 (MQ=255)
          ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgc                                                                                         >  1:2285288/1‑100 (MQ=255)
          ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGAAGCCGATATTAAAAAAAAAGCCCCCTGCGAAGGGGgc                                                                                         <  2:2340924/100‑1 (MQ=255)
           gCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgca                                                                                        <  2:211048/100‑1 (MQ=255)
            cAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaa                                                                                       >  1:1459391/1‑100 (MQ=255)
            cAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgc                                                                                         >  1:808276/1‑98 (MQ=255)
             aGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaa                                                                                      <  2:462699/100‑1 (MQ=255)
              gcaGTAGAACACGCGCTGCGGAAATTTCCTTACTCGCAAGATACCTTCTCTAACTCCGCAAACGTCTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaag                                                                                     <  1:344717/98‑1 (MQ=255)
              gAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaag                                                                                     <  1:2348245/100‑1 (MQ=255)
              gAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaag                                                                                     <  2:1164749/100‑1 (MQ=255)
                 gTAGACCTCGCGCTGCGCAAATTTCCTTCCTCGGGCGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAACAATGCCCCCTGCGAAGGGGgcaaag                                                                                     <  1:2482222/97‑1 (MQ=255)
                 gTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagcaa                                                                                  >  2:1522775/1‑100 (MQ=255)
                 gTAGAACACGCGCGGCGGAAATTTCCTTCCTCGGGAGATAACTGCTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgcttttttt                                                                                  <  1:1585193/100‑8 (MQ=255)
                  tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgctttttttt                                                                                 <  1:287678/100‑9 (MQ=38)
                  tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgctttttttt                                                                                 <  2:1387196/100‑9 (MQ=38)
                  tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgctttttttt                                                                                 <  1:2115511/100‑9 (MQ=38)
                  tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgctttttttt                                                                                 <  1:1400732/100‑9 (MQ=38)
                  tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgctttttttt                                                                                 <  1:1108891/100‑9 (MQ=38)
                  tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgctttttttt                                                                                 <  2:784568/100‑9 (MQ=38)
                   aGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgcttttttttg                                                                                <  1:924178/100‑10 (MQ=38)
                   aGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgctttttttt                                                                                 <  1:2228309/99‑9 (MQ=38)
                   aGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAAc                                                                                >  2:2031757/1‑100 (MQ=255)
                    gAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgcttttttttga                                                                               <  2:547199/100‑11 (MQ=38)
                    gAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACt                                                                               <  2:1300511/100‑1 (MQ=255)
                     aaCACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTa                                                                              <  1:190991/100‑1 (MQ=255)
                       cacGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATg                                                                            <  2:2059903/100‑1 (MQ=255)
                         cgcgCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGc                                                                          <  2:1459399/100‑1 (MQ=255)
                          gcgcTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCa                                                                         <  1:1096846/100‑1 (MQ=255)
                               gCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTTTATAAAAATCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGtt                                                                    >  1:1761283/1‑100 (MQ=255)
                               gCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGtt                                                                    >  1:692839/1‑100 (MQ=255)
                                cGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCGCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGttt                                                                   >  2:1815297/1‑100 (MQ=255)
                                cGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGttt                                                                   >  2:451182/1‑100 (MQ=255)
                                 ggAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTc                                                                  >  2:602163/1‑100 (MQ=255)
                                 ggAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTc                                                                  >  2:440840/1‑100 (MQ=255)
                                    aaTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGCTTCGtt                                                               >  1:1669419/1‑100 (MQ=255)
                                    aaTTTCCTTCCTCGGGAGATAACTGGCCTAAGTCGGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGtt                                                               <  1:1299640/100‑1 (MQ=255)
                                     aTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTg                                                              <  2:742943/100‑1 (MQ=255)
                                      tttCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTgg                                                             <  2:947231/100‑1 (MQ=255)
                                      tttCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTgg                                                             <  2:2285303/100‑1 (MQ=255)
                                      tttCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTgg                                                             <  2:1664780/100‑1 (MQ=255)
                                       ttccttccTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt                                                            <  2:1346402/100‑1 (MQ=255)
                                         ccttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTta                                                          >  2:653301/1‑100 (MQ=255)
                                          cttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGtaga                                                         >  2:1326596/1‑97 (MQ=255)
                                           ttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtata                                                        <  2:1570177/100‑1 (MQ=255)
                                             cctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATAcc                                                      <  1:2362007/100‑1 (MQ=255)
                                              ctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat                                                         <  1:521961/96‑1 (MQ=255)
                                              ctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCt                                                     >  2:1359220/1‑100 (MQ=255)
                                               tCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTg                                                    >  1:1447907/1‑100 (MQ=255)
                                                 gggAGTTAATTGGTTTAATTCCGAAGCGTTTTTTAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTgg                                                   <  2:2429424/99‑1 (MQ=255)
                                                 gggAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt                                                  <  1:442186/100‑1 (MQ=255)
                                                   gagaTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGc                                                <  2:2356884/100‑1 (MQ=255)
                                                    agaTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCt                                               >  2:1457913/1‑100 (MQ=255)
                                                    agaTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCt                                               <  2:2559896/100‑1 (MQ=255)
                                                       tAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGc                                            <  1:333899/100‑1 (MQ=255)
                                                       tAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGc                                            >  1:958513/1‑100 (MQ=255)
                                                         aCTGTTCGAATGCCGCAGCCGTATTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAgt                                            <  1:183060/98‑2 (MQ=255)
                                                         aCTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGa                                          <  1:734359/100‑1 (MQ=255)
                                                          cTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGAt                                         >  2:2026318/1‑100 (MQ=255)
                                                          cTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGAt                                         >  2:194038/1‑100 (MQ=255)
                                                           tGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATa                                        <  2:1496061/100‑1 (MQ=255)
                                                            ggTCTAATTCCGCAGCCGTTTTTCAAAACAAAGCCCCCTGCGAAGGGGACAAAGCAAACTATGGCAATGTTCCGTTGGTTATACCTGGTGCTAGCGATa                                        >  1:621100/1‑99 (MQ=255)
                                                              tCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAat                                     >  2:749949/1‑100 (MQ=255)
                                                              tCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAat                                     >  1:219041/1‑100 (MQ=255)
                                                              tCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAat                                     >  2:2412825/1‑100 (MQ=255)
                                                               cTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAata                                    >  2:2440661/1‑100 (MQ=255)
                                                               cTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAAATATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAata                                    >  1:1655792/1‑100 (MQ=255)
                                                                   ttCCGCAGCCGTTTTTTTTTAAAAATCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTCCTAGCGATAAATATCcg                                >  1:175764/1‑100 (MQ=255)
                                                                       gaaGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCt                            >  2:1758924/3‑100 (MQ=255)
                                                                             gTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTg                      <  1:1966841/100‑1 (MQ=255)
                                                                              tttttCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGc                     <  1:2590824/100‑1 (MQ=255)
                                                                               ttttCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCa                    >  2:1515468/1‑100 (MQ=255)
                                                                                tttCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCat                   <  2:283049/100‑1 (MQ=255)
                                                                                tttCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCat                   >  2:1657441/1‑100 (MQ=255)
                                                                                  aaaaaaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGtn                      <  1:2138644/93‑2 (MQ=255)
                                                                                   cAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATAtt                >  1:1464307/1‑100 (MQ=255)
                                                                                   cAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATAtt                >  2:1328337/1‑100 (MQ=255)
                                                                                    aaaaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTg               >  2:943298/1‑100 (MQ=255)
                                                                                     aaaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGc              >  1:1770383/1‑100 (MQ=255)
                                                                                      tgctttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGcc             >  2:1629906/7‑100 (MQ=255)
                                                                                      tgctttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAATATCCGCGCTGGCGTGCATATTGCCg            >  1:2083279/7‑100 (MQ=255)
                                                                                      aaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGcc             <  1:1020454/100‑1 (MQ=255)
                                                                                       gctttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAACGATAAATATCCGCGCTGGCGTGCATATTGCCg            >  2:267237/6‑100 (MQ=255)
                                                                                        ctttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCTCTGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGt           >  1:759922/5‑100 (MQ=255)
                                                                                        ctttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGt           >  2:651572/5‑100 (MQ=255)
                                                                                        ctttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGt           >  2:968747/5‑100 (MQ=255)
                                                                                           aGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTAc        <  2:1899463/100‑1 (MQ=255)
                                                                                           aGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTAc        >  2:1862854/1‑100 (MQ=255)
                                                                                            gCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACt       >  1:936993/1‑100 (MQ=255)
                                                                                               cccTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTccc    >  2:1046353/1‑100 (MQ=255)
                                                                                               cccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATGCCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTc      <  1:1568601/98‑1 (MQ=255)
                                                                                               cccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTctc    >  1:1223497/1‑98 (MQ=255)
                                                                                                ccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTcccc   >  2:308486/1‑100 (MQ=255)
                                                                                                ccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTcccc   <  1:1189149/100‑1 (MQ=255)
                                                                                                 cTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTaagatcgg  >  2:1189367/1‑93 (MQ=255)
                                                                                                   |                                                                                                 
GCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCG  >  NC_000913/4424366‑4424562

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: