Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,142 T→A 17.8% intergenic (+6/‑69) yejM → / → proL essential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,1420TA17.8% 359.7 / 63.5 145intergenic (+6/‑69)yejM/proLessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (54/65);  new base A (11/15);  total (65/80)
Fisher's exact test for biased strand distribution p-value = 8.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 160 total aligned reads displayed.

CAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACC  >  NC_000913/2286043‑2286239
                                                                                                   |                                                                                                 
   cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatta                                                                                                  >  1:2133411/1‑98 (MQ=255)
    ttACGTGTTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaa                                                                                               <  2:3085284/100‑1 (MQ=255)
      aCGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAAt                                                                                             <  1:3405003/100‑1 (MQ=255)
       cGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                            <  1:3110088/100‑1 (MQ=255)
        gtggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCa                                                                                           >  1:1505334/1‑100 (MQ=255)
        gtggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCa                                                                                           >  1:1504850/1‑100 (MQ=255)
         tggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                          <  1:737711/100‑1 (MQ=255)
          ggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                                                         >  2:191358/1‑100 (MQ=255)
            tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGATGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTa                                                                                       <  2:1953458/100‑1 (MQ=255)
            tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTa                                                                                       <  2:3354086/100‑1 (MQ=38)
             gaaagagtgacagaTGGAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCCCTAACTGATTAATTATAAATCAGTTAg                                                                                      >  2:724297/1‑100 (MQ=255)
             gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAg                                                                                      <  2:777062/100‑1 (MQ=38)
             gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaaagatcgga                                                                                      >  2:1636474/1‑92 (MQ=255)
             gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                      >  2:1291759/1‑100 (MQ=255)
              aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                     <  1:1543270/100‑1 (MQ=38)
              aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                     <  2:1484958/100‑1 (MQ=38)
               aagagtgaaagaTGAAAAACCACAGTTACGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                    <  1:3027119/100‑1 (MQ=38)
               aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                    <  2:43131/100‑1 (MQ=38)
                agagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGa                                                                                   <  2:3263166/100‑1 (MQ=255)
                 gagtgaatgaTAAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACCGATTTATAATTAATCAGTTAGc                                                                                     <  2:1097025/97‑1 (MQ=37)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAgatcg                                                                                  >  1:1700925/1‑96 (MQ=38)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                                                  <  1:2894545/100‑1 (MQ=255)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                                                  >  2:2053361/1‑100 (MQ=255)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                                                  >  2:2951043/1‑100 (MQ=255)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATCGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                                                  <  1:2555661/100‑1 (MQ=255)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTAGTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCGat                                                                                  <  2:1042496/100‑2 (MQ=37)
                     gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                                                              >  1:2286288/1‑100 (MQ=255)
                     gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                                                              >  2:3575390/1‑100 (MQ=255)
                     gaaagaTGAAAAACAACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATGAGTTATAAATCAGTTAGCGaaa                                                                                 >  2:1699656/1‑97 (MQ=255)
                        acaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                    <  2:2809690/89‑1 (MQ=38)
                            gAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACtt                                                                       <  2:1164779/100‑1 (MQ=255)
                             aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                      <  1:2380202/100‑1 (MQ=255)
                             aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                      <  1:309686/100‑1 (MQ=255)
                             aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                      <  1:139467/100‑1 (MQ=255)
                             aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                      >  1:2579147/1‑100 (MQ=255)
                              aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                     <  1:715885/100‑1 (MQ=255)
                              aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                     >  1:2532144/1‑100 (MQ=255)
                               aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCa                                                                    <  2:409010/100‑1 (MQ=255)
                                  ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATc                                                                 >  1:822013/1‑100 (MQ=255)
                                   cacaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                                >  1:1793011/1‑100 (MQ=255)
                                   cacaGTGAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                                <  1:2120429/100‑1 (MQ=255)
                                    acaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCgn                                                               <  1:166389/100‑2 (MQ=255)
                                     caGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgt                                                              >  1:93675/1‑100 (MQ=255)
                                        ttAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGtgtg                                                           <  2:2155393/100‑1 (MQ=255)
                                         tAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTgg                                                          <  1:2642149/100‑1 (MQ=255)
                                            gTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaag                                                       >  1:3610136/1‑99 (MQ=255)
                                              ttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAACTCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAAc                                                     <  2:1990246/100‑1 (MQ=255)
                                              ttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAAc                                                     <  2:3207583/100‑1 (MQ=255)
                                              ttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAACCGGTGTGGAAAAc                                                     >  2:1744349/1‑100 (MQ=255)
                                               tGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGCTTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACg                                                    >  2:2527762/1‑100 (MQ=255)
                                                gTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                   <  2:2356447/100‑1 (MQ=255)
                                                gTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                   <  2:1241820/100‑1 (MQ=255)
                                                cgatctGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                   <  1:271932/95‑1 (MQ=255)
                                                 ttATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                  <  1:279432/100‑1 (MQ=255)
                                                 ttATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                  <  1:528656/100‑1 (MQ=255)
                                                  tATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgta                                                 <  1:1310319/100‑1 (MQ=255)
                                                   aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                                <  2:2080677/100‑1 (MQ=255)
                                                     tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTACTTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                              >  1:2728849/1‑100 (MQ=255)
                                                      gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACggagatcg                                             >  2:271930/1‑94 (MQ=255)
                                                      gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATATTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                             >  1:1079534/1‑100 (MQ=255)
                                                         aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTag                                          <  1:1940624/100‑1 (MQ=255)
                                                          tctcTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                         <  2:3624051/97‑1 (MQ=255)
                                                          gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                         <  1:1076496/100‑1 (MQ=255)
                                                          gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagc                                         <  2:2192243/100‑1 (MQ=255)
                                                             cTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagcagc                                      >  1:2541069/1‑100 (MQ=255)
                                                             cTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagcagc                                      >  1:2350791/1‑100 (MQ=255)
                                                             cTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagcag                                       >  2:1609005/1‑99 (MQ=255)
                                                               gACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCa                                    >  1:2023634/1‑100 (MQ=255)
                                                                 cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                  <  2:2837960/100‑1 (MQ=255)
                                                                 cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                  <  2:3584529/100‑1 (MQ=255)
                                                                 cAGACGAGAAGCGTTTTATCGCTAACGGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                  <  2:11019/100‑1 (MQ=255)
                                                                         aaGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCa                          <  1:916965/100‑1 (MQ=255)
                                                                         aaGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCa                          <  1:274954/100‑1 (MQ=255)
                                                                          aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCAc                         >  1:404607/1‑100 (MQ=255)
                                                                            cGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGt                       <  2:1766085/100‑1 (MQ=255)
                                                                            agatatttCGCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGt                       >  2:916702/8‑100 (MQ=25)
                                                                             gTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTa                      <  1:2191128/100‑1 (MQ=255)
                                                                             gTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATAGCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTa                      <  2:2338774/100‑1 (MQ=255)
                                                                                ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                   <  1:2843768/100‑1 (MQ=255)
                                                                                  ttCGCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCa                 >  1:683535/2‑100 (MQ=38)
                                                                                  aTCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCa                 <  1:2174379/100‑1 (MQ=255)
                                                                                   tCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                >  1:560754/1‑100 (MQ=255)
                                                                                   tCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                <  1:425251/100‑1 (MQ=255)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCTATCGGAGTGGAAAACGGAAGTATTAGCAGCCACGAGACGGCACGTAGCGCAGc               >  2:1543682/1‑99 (MQ=37)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc              >  1:1370825/1‑100 (MQ=38)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATAGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc              >  1:2851866/1‑100 (MQ=38)
                                                                                     gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc              >  1:1669327/1‑100 (MQ=255)
                                                                                     gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGc               >  2:3499526/1‑99 (MQ=255)
                                                                                      cTAACTGATTTATAATTAATAAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt             >  2:2223204/1‑100 (MQ=38)
                                                                                        aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                >  2:3045265/1‑95 (MQ=255)
                                                                                             attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcg      <  1:2750469/100‑1 (MQ=255)
                                                                                             attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcg      <  1:1331669/100‑1 (MQ=255)
                                                                                             attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAGTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcg      <  2:3166284/100‑1 (MQ=255)
                                                                                              ttaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcgc     >  1:3348722/1‑100 (MQ=255)
                                                                                               taattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCa    <  1:1706644/100‑1 (MQ=255)
                                                                                                aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc   >  1:2945760/1‑100 (MQ=255)
                                                                                                aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc   <  1:191358/100‑1 (MQ=255)
                                                                                                aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc   >  1:1434361/1‑100 (MQ=255)
                                                                                                 attaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAcc  >  1:377338/1‑100 (MQ=255)
                                                                                                 attaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAcc  >  2:272667/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
CAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACC  >  NC_000913/2286043‑2286239

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