Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,338,285 A→G 15.1% S59S (AGT→AGC mlaD ← OM lipid asymmetry maintenance protein; membrane‑anchored ABC family periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,338,2850AG15.1% 382.1 / 48.1 158S59S (AGT→AGCmlaDOM lipid asymmetry maintenance protein; membrane‑anchored ABC family periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (73/59);  new base G (11/13);  total (86/72)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

TGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACGTTCGCCGCCTTC  >  NC_000913/3338187‑3338384
                                                                                                  |                                                                                                   
tgGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATAc                                                                                                    <  1:721632/100‑1 (MQ=255)
tgGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATAc                                                                                                    <  2:1533088/100‑1 (MQ=255)
 ggTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATGCt                                                                                                   <  1:2111056/100‑1 (MQ=255)
 ggTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACt                                                                                                   >  2:2156905/1‑100 (MQ=255)
  gTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTg                                                                                                  <  2:107538/100‑1 (MQ=255)
  gTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCAACCAATACTg                                                                                                  >  1:2775179/1‑100 (MQ=255)
   ttATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGa                                                                                                 >  2:2390760/1‑100 (MQ=255)
      taACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATGCTGACCg                                                                                              >  2:3181528/1‑100 (MQ=255)
      taACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCg                                                                                              >  1:1304209/1‑100 (MQ=255)
      taACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCg                                                                                              >  2:2156604/1‑100 (MQ=255)
      taACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCg                                                                                              >  1:1761343/1‑100 (MQ=255)
        aCGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCACAACACGACACACAACAACGCCACCAATCCTGACCGga                                                                                            >  1:2193091/1‑100 (MQ=255)
             gttCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAgcg                                                                                       <  2:2534418/100‑1 (MQ=255)
              ttCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACCACGCCACCAATACTGACCGGAGAgcgc                                                                                      >  1:559160/1‑100 (MQ=255)
              ttCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCACAAACAACGCCACCAATACTGACCGGAGAgcgc                                                                                      <  1:314978/100‑1 (MQ=255)
               tCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATGCTGACCGGAGAgcgcg                                                                                     <  1:28623/100‑1 (MQ=255)
               tCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATGCTGACCGGAGAgcgc                                                                                      >  2:803352/1‑99 (MQ=255)
               tCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAgcgcg                                                                                     >  1:1909012/1‑100 (MQ=255)
                cAATTTGCAGCGTTACGCGCGGCAGTGAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAgcgcgc                                                                                    <  2:1432975/100‑1 (MQ=255)
                cAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAgcgcgc                                                                                    <  2:3534397/100‑1 (MQ=255)
                 aaggtCCAGCGTTACGCGCGGCAGGTAGGTTTGCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAgcgcg                                                                                     <  1:3274822/94‑1 (MQ=255)
                 aaTTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACCACAACGCCACCAATGCTGACCGGAGAgcgc                                                                                      >  1:2178791/1‑97 (MQ=255)
                 aaTTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCt                                                                                   >  2:1026979/1‑100 (MQ=255)
                 aaTTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCt                                                                                   >  2:579450/1‑100 (MQ=255)
                    ttCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTc                                                                                <  2:201919/100‑1 (MQ=255)
                     tCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTc                                                                                >  1:2829701/1‑99 (MQ=255)
                       cAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATGCTGACCGGAGAGCGCGCTTTCAgg                                                                             <  1:3127394/100‑1 (MQ=255)
                       cAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAgg                                                                             <  2:551984/100‑1 (MQ=255)
                        aGCTTTACGTGTGGCAGATAGGTTTTCTGGTCTATAGTAATATCCGCCACCAGACCAACAACAACGCCAACAATGCTGACAGGAGAGCGCGCTTTCAggt                                                                            <  1:3260818/100‑2 (MQ=255)
                        aGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGgc                                                                            <  2:1954391/100‑1 (MQ=255)
                           gTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAAAAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGgccgc                                                                         >  1:1870916/1‑100 (MQ=255)
                            ttACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGgccgcc                                                                        <  2:2728564/100‑1 (MQ=255)
                            ttACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGgccgcc                                                                        <  1:3638530/100‑1 (MQ=255)
                             tACGCGCGTCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCTATACTGACCGGAGAGCGCGCTTTCAGGCCGCCa                                                                       <  1:2508193/100‑1 (MQ=255)
                               cgcgcgGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAAt                                                                     >  1:1342053/1‑100 (MQ=255)
                                 cgcgGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGt                                                                   <  2:620799/100‑1 (MQ=255)
                                 cgcgGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGt                                                                   >  2:2884063/1‑100 (MQ=255)
                                  gcgGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATGCTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGt                                                                   <  1:2446964/99‑1 (MQ=255)
                                     gCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATGCTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATc                                                               >  1:1795206/1‑100 (MQ=255)
                                     gCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATc                                                               >  2:1225980/1‑100 (MQ=255)
                                     gCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATc                                                               <  2:3655962/100‑1 (MQ=255)
                                      cAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCg                                                              <  2:990148/100‑1 (MQ=255)
                                          tAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACg                                                          >  2:1507851/1‑100 (MQ=255)
                                            ggTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTc                                                        >  1:2378838/1‑100 (MQ=255)
                                            ggTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTc                                                        >  2:2626025/1‑100 (MQ=255)
                                            ggTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTc                                                        >  2:3013104/1‑100 (MQ=255)
                                            ggTTTTCGGGCCCAGCGTAATATCCGCCCCCCGCCCCACAACAACGCCACCAATACTGACCGGAGAGCACGCTTTCAGGCCGCCAATGTTATCGAACGTc                                                        <  2:1946823/100‑1 (MQ=255)
                                             gTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATGCTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCg                                                       >  1:156842/1‑100 (MQ=255)
                                             gTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCg                                                       <  2:1117338/100‑1 (MQ=255)
                                              ttttCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGc                                                      <  2:3758059/100‑1 (MQ=255)
                                                ttCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCAt                                                    <  1:180911/100‑1 (MQ=255)
                                                 tCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATa                                                   >  1:2509062/1‑100 (MQ=255)
                                                   gggTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATaaa                                                 >  1:726782/1‑100 (MQ=255)
                                                    ggTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAg                                                <  2:3298265/100‑1 (MQ=255)
                                                    ggTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAg                                                <  1:2124398/100‑1 (MQ=255)
                                                     gTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGc                                               >  1:3022786/1‑100 (MQ=255)
                                                      tCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCg                                              <  2:2798680/100‑1 (MQ=255)
                                                        cAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATGCTGACCGGAGCGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCg                                              <  2:937695/98‑1 (MQ=255)
                                                        cAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCgtg                                            >  2:3036130/1‑100 (MQ=255)
                                                          gCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTa                                          <  1:2152632/100‑1 (MQ=255)
                                                               atatCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCg                                     >  2:432597/1‑100 (MQ=255)
                                                               atatCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCg                                     >  1:1394137/1‑100 (MQ=255)
                                                               atatCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCg                                     >  2:931427/1‑100 (MQ=255)
                                                                tatCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCgg                                    >  2:1184288/1‑100 (MQ=255)
                                                                  tCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGtt                                  <  2:869801/100‑1 (MQ=255)
                                                                   ccGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTc                                 >  2:1185402/1‑100 (MQ=255)
                                                                    cGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCa                                >  2:732003/1‑100 (MQ=255)
                                                                       cACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTa                             <  1:1040093/100‑1 (MQ=255)
                                                                        aCCCGACCCACAACAACGCCACCAATACGGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTa                             <  1:1037808/99‑1 (MQ=255)
                                                                         cccGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACg                           >  2:1818643/1‑100 (MQ=255)
                                                                          ncGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGt                          >  2:142010/2‑100 (MQ=255)
                                                                          ccGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGt                          >  2:3656163/1‑100 (MQ=255)
                                                                           cGACCCACAACAACGCCACCAATACTGACCGGCGAGCGCACTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTa                         <  1:2670600/100‑1 (MQ=255)
                                                                            gACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTAt                        <  1:2044349/100‑1 (MQ=255)
                                                                              cccACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGCTTCAGTACGTAt                        >  1:3751823/1‑98 (MQ=255)
                                                                               ccACAACAACGCCACCAATGCTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGa                     >  1:1629195/1‑100 (MQ=255)
                                                                                cacaACAACGCCACCAATGCTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGAc                    >  1:406693/1‑100 (MQ=255)
                                                                                cacaACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGAc                    >  2:2886820/1‑100 (MQ=255)
                                                                                cacaACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGAc                    >  1:1680110/1‑100 (MQ=255)
                                                                                cacaACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGAc                    >  1:1309534/1‑100 (MQ=255)
                                                                                  caacaaCGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGt                  >  1:826992/1‑100 (MQ=255)
                                                                                  caacaaCGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGt                  >  2:1728919/1‑100 (MQ=255)
                                                                                   aacaacGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTc                 <  1:3533508/100‑1 (MQ=255)
                                                                                   aacaacGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTc                 <  2:2118476/100‑1 (MQ=255)
                                                                                   aacaacGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTc                 <  1:2591304/100‑1 (MQ=255)
                                                                                    acaacGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCa                >  1:2369608/1‑100 (MQ=255)
                                                                                    acaacGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCa                <  2:3455401/100‑1 (MQ=255)
                                                                                    acaacGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCa                >  2:3634298/1‑100 (MQ=255)
                                                                                      aacGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACg              >  1:3128575/1‑100 (MQ=255)
                                                                                        cGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACGtt            >  1:1374342/1‑100 (MQ=255)
                                                                                          ccaccaATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCTGTTCGGTACGTATGGCCGTCACGTTc           >  1:71435/1‑99 (MQ=255)
                                                                                           caccaATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACGTTcgc         >  1:2632250/1‑100 (MQ=255)
                                                                                           caccaATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACGTTcgc         <  1:894026/100‑1 (MQ=255)
                                                                                             ccaATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACGTTcgccg       >  2:1636289/1‑100 (MQ=255)
                                                                                             ccaATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACGTTcgccg       >  2:2035868/1‑100 (MQ=255)
                                                                                             ccaATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACGTTcgccg       >  1:3648585/1‑100 (MQ=255)
                                                                                               aaTACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTATGTCGGTTCAGTACGTATGGACGTCACGTTCgccgcc     <  2:2680149/100‑1 (MQ=255)
                                                                                               aaTACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACGTTCgccgcc     >  2:2429247/1‑100 (MQ=255)
                                                                                                aTACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACGTTCGCCGCCt    <  2:2504898/100‑1 (MQ=255)
                                                                                                  aCTGACCGGAGAGCGCGCTTTAAGGCCGCCAATGTTATCGAACGTCGAATGACGCGGGTAGGTCGGTTCAGTACGCATGGACGTAACGTTCAcctccctc  >  1:3760973/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGGTTATAACGTTGTTCAATTTCCAGCGTTACGCGCGGCAGATAGGTTTTCGGGTCCAGCGTAATATCCGCCACCCGACCCACAACAACGCCACCAATACTGACCGGAGAGCGCGCTTTCAGGCCGCCAATGTTATCGAACGTCGCATAAAGCGTGTAGGTCGGTTCAGTACGTATGGACGTCACGTTCGCCGCCTTC  >  NC_000913/3338187‑3338384

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: