Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 407,203 Δ1 bp 19.7% intergenic (+33/‑225) yaiA → / → aroM OxyR‑regulated conserved protein/AroM family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913407,2030G.19.7% 330.1 / 88.6 122intergenic (+33/‑225)yaiA/aroMOxyR‑regulated conserved protein/AroM family protein
Reads supporting (aligned to +/- strand):  ref base G (46/52);  new base . (10/14);  total (56/66)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

Only 100 of 126 total aligned reads displayed.

TGAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTGGAATC  >  NC_000913/407104‑407302
                                                                                                   |                                                                                                   
tGAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTg                                                                                                     <  2:891942/100‑1 (MQ=255)
 gAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGCGGACCCTGACAACGAGTGACCTCATCAGACTGTTCGGAAGGGATTCTGAGTGc                                                                                                    >  1:544612/1‑100 (MQ=255)
 gAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGc                                                                                                    >  2:332513/1‑100 (MQ=255)
  aaCATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGcc                                                                                                   <  2:905925/100‑1 (MQ=255)
  aaCATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGcc                                                                                                   <  1:1268237/100‑1 (MQ=255)
    cATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCCCATCAGACTGCTCGGAAGGGATTCTGAGTGCCAc                                                                                                 >  1:85627/1‑100 (MQ=255)
     aTCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTCCACTa                                                                                               >  2:753137/1‑100 (MQ=255)
      tCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTa                                                                                               >  1:2244655/1‑100 (MQ=255)
        cGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACa                                                                                             >  1:1230596/1‑100 (MQ=255)
           ccGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAgg                                                                                          >  1:783249/1‑100 (MQ=255)
           ccGCTCAGGAAGCGATGAATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTCCACTACAAggg                                                                                         >  2:458612/1‑100 (MQ=255)
             gCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGa                                                                                        <  1:1286775/100‑1 (MQ=255)
             gCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGa                                                                                        >  1:1758564/1‑100 (MQ=255)
             gCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGa                                                                                        >  1:1761092/1‑100 (MQ=255)
             gCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGa                                                                                        >  2:601524/1‑100 (MQ=255)
              cTCAGGAAGCGATGGAGGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAg                                                                                           <  2:2180268/96‑1 (MQ=255)
               tCAGGACGCGATGGATGCCAACAAACGCTAAGAGGACCCTGACAACGCCTGACCGCATCCCTCCGCTCGGAAGGGATTCTGCGTGACACTACa                                                                                             <  2:1072885/93‑1 (MQ=255)
               tCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATc                                                                                      >  2:387420/1‑100 (MQ=255)
                cAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTCCACTACAAGGGATCTg                                                                                    >  1:1912465/1‑100 (MQ=255)
                 aGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTg                                                                                    >  1:373230/1‑100 (MQ=255)
                   gAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCg                                                                                  <  2:67714/100‑1 (MQ=255)
                   gAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCg                                                                                  >  1:1064421/1‑100 (MQ=255)
                    aaGCGATGGATGCGAAAGAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTCCACTACAAGGGATCTGCGTc                                                                                <  1:106245/100‑1 (MQ=255)
                    aaGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGt                                                                                 >  2:544974/1‑100 (MQ=255)
                    aaGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTCCACTACAAGGGATCTGCGTc                                                                                <  2:1821147/100‑1 (MQ=255)
                      gCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTca                                                                               >  1:1897119/1‑100 (MQ=255)
                        atctgatgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTcaca                                                                             <  1:840818/96‑1 (MQ=255)
                           ggatgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACAttt                                                                          <  2:722784/100‑1 (MQ=255)
                            gatgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTcacaagat                                                                         >  2:840809/1‑96 (MQ=255)
                            gatgCGAAAAAACGCTATGAGGACACTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCCTCACAtttt                                                                         >  2:2294075/1‑100 (MQ=255)
                             atgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACAttttt                                                                        <  1:1397517/100‑1 (MQ=255)
                              tgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTc                                                                       <  1:159598/100‑1 (MQ=255)
                              tgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTc                                                                       <  1:544977/100‑1 (MQ=255)
                              tgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTc                                                                       >  2:843078/1‑100 (MQ=255)
                              tgCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTCCACTACAAGGGATCTGCGTCACATTTTTCa                                                                      <  2:1962732/100‑1 (MQ=255)
                                cGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCAt                                                                     <  1:523595/100‑1 (MQ=255)
                                 gAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATa                                                                    >  2:899824/1‑100 (MQ=255)
                                     aaaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTc                                                                <  2:1775421/100‑1 (MQ=255)
                                      aaCGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCa                                                               <  2:1857091/100‑1 (MQ=255)
                                         gCTATGAGGACCCTGACAAAGAGTGACCGCATCAGAATGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGt                                                            >  1:1689047/1‑100 (MQ=255)
                                           tATGAGTACCCTGACATAGAGTGACCGCATCAGACTGCTCTGAAGGTATTCTGAGTCCACTACAAAGGATCTGCGTCACATTTTTCATAATTCATGtttt                                                         <  2:2001078/100‑1 (MQ=255)
                                                ggACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCtaa                                                    >  2:1241180/1‑100 (MQ=255)
                                                 gACACTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTTTTATTGCCACTACAATGTATCTGCGTCTCATTTTTCATAGTTCATGTTTTTc                                                       <  2:396231/97‑1 (MQ=255)
                                                   cccTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCtaata                                                  >  1:2016780/1‑100 (MQ=255)
                                                    ccTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCtaataa                                                 <  1:2234473/100‑1 (MQ=255)
                                                    ccTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCtaataa                                                 >  1:1371469/1‑100 (MQ=255)
                                                       gACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTa                                              >  2:1879353/1‑100 (MQ=255)
                                                        aCAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAAtt                                               >  2:2082737/1‑98 (MQ=255)
                                                        aCAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAg                                             <  2:1352740/100‑1 (MQ=255)
                                                         cAAAGAGTGACCTCATCAGACTACCCGGAAGGGCTTCTGAGTGCCACTACAAGAGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGa                                            >  2:1437083/1‑100 (MQ=255)
                                                           aaGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAat                                          <  2:2199505/100‑1 (MQ=255)
                                                            agagTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAata                                         <  2:1469127/100‑1 (MQ=255)
                                                             gagTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATATGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAatat                                        >  1:1178651/1‑100 (MQ=255)
                                                                tGACCGCATCAGACTGCTCGGAAGGGATTTTGAGTGCCACTACTATGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATAtt                                       <  2:1484763/98‑1 (MQ=255)
                                                                  aCCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaac                                   <  1:875600/100‑1 (MQ=255)
                                                                  aCCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaac                                   <  1:1520800/100‑1 (MQ=255)
                                                                  aCCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaac                                   <  2:1623439/100‑1 (MQ=255)
                                                                   ccGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaaca                                  >  2:1713332/1‑100 (MQ=255)
                                                                   ccGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACCACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaaca                                  >  2:1276985/1‑100 (MQ=255)
                                                                   ccGCATCAGACTGCTCGGAAGGGATTCTGAGTCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaacaa                                 >  2:436102/1‑100 (MQ=255)
                                                                     gCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaacaat                                >  1:449232/1‑100 (MQ=255)
                                                                     gCATCAGACTGCTCGGAAGGATTTCTAAGTCCACTACTAGGTATATGCGTCACATTTTTTATAATTCATGTTTTTCTAATAATTAGAATATTAaacaata                               <  2:811261/100‑1 (MQ=255)
                                                                      cATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaacaata                               >  2:1969093/1‑100 (MQ=255)
                                                                       aTCAGACTGCTCGTAATGGATTCTGTGTGACACTACAAGGTATCTTTGTCACATTTTTCATAATTCATGTTTTTCTAATATTTAGAATATTAaacaataa                              <  2:1124457/100‑1 (MQ=255)
                                                                        tCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaacaataac                             <  1:905116/100‑1 (MQ=255)
                                                                        tCAGACTGCTCGGAAGGGATTCTGAGTCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaacaataaca                            <  1:1870384/100‑1 (MQ=255)
                                                                        tCAGACTGCTCGGAAGGGATTCTGAGTCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAaacaataaca                            <  1:753145/100‑1 (MQ=255)
                                                                           gACTTCTCGGAAGGGATTATGAGTCCACTACAAGTGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATc                         <  2:1014710/100‑1 (MQ=255)
                                                                            aCTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATc                         <  1:442804/100‑1 (MQ=255)
                                                                             cTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATcc                        <  2:1985977/100‑1 (MQ=255)
                                                                              tGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATc                         <  2:1453387/98‑1 (MQ=255)
                                                                              tGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATc                         <  1:671222/98‑1 (MQ=255)
                                                                              tGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCa                       <  1:1111703/100‑1 (MQ=255)
                                                                               gCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCAt                      >  1:37274/1‑100 (MQ=255)
                                                                               gCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATACAt                      >  2:823977/1‑100 (MQ=255)
                                                                                cTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCAtt                     >  2:453495/1‑100 (MQ=255)
                                                                                  cGGAAGGGATTCTGAGTGCCACTATAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTAc                   <  2:1032377/100‑1 (MQ=255)
                                                                                  cGGAAGGGATTCTGAGTCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACt                  >  2:1003035/1‑100 (MQ=255)
                                                                                   ggAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACt                  >  2:559346/1‑100 (MQ=255)
                                                                                   ggAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACt                  >  1:2192297/1‑100 (MQ=255)
                                                                                      aGGGATTCTGAGTGCCACTACAAGTGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGa               <  2:2207497/100‑1 (MQ=255)
                                                                                       gggATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGaa              >  2:1231576/1‑100 (MQ=255)
                                                                                        ggATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAAt             >  2:1649217/1‑100 (MQ=255)
                                                                                        ggATTCTGAGTCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATc            <  1:1844475/100‑1 (MQ=255)
                                                                                         gATTCTGAGTGCCAATACAAGGGAACTGAGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATACATTAATTGAAt             >  2:455239/1‑99 (MQ=255)
                                                                                          aTTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGaa              <  1:1057856/97‑1 (MQ=255)
                                                                                          aTTCTGAGTGCCACTACAAGGAATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCa           >  2:1921103/1‑100 (MQ=255)
                                                                                             cTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCAttt        >  1:536397/1‑100 (MQ=255)
                                                                                             cTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTAATGGAATCAtt         >  2:1651264/1‑99 (MQ=255)
                                                                                              tGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTg       <  2:1485646/100‑1 (MQ=255)
                                                                                              tGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTg       <  1:1467626/100‑1 (MQ=255)
                                                                                               gAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTgg      <  2:1412947/100‑1 (MQ=255)
                                                                                               gAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTgg      <  2:725727/100‑1 (MQ=255)
                                                                                               gAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTgg      <  1:1596256/100‑1 (MQ=255)
                                                                                               gAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTgg      <  1:1565381/100‑1 (MQ=255)
                                                                                                 gTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTGGaa    >  2:1212157/1‑100 (MQ=255)
                                                                                                  tGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTGGAAt   <  1:103504/100‑1 (MQ=255)
                                                                                                   gCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTGGAATc  >  2:156947/1‑100 (MQ=255)
                                                                                                   gCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTGGAATc  <  2:544609/100‑1 (MQ=255)
                                                                                                   gCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTAATGGAATCATTTGGa     >  2:2348696/1‑97 (MQ=255)
                                                                                                   |                                                                                                   
TGAACATCCGACCGCTCAGGAAGCGATGGATGCGAAAAAACGCTATGAGGACCCTGACAAAGAGTGACCGCATCAGACTGCTCGGAAGGGATTCTGAGTGCCACTACAAGGGATCTGCGTCACATTTTTCATAATTCATGTTTTTCTAATAATTAGAATATTAAACAATAACAATCCATTACTGGAATCATTTGGAATC  >  NC_000913/407104‑407302

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: