Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,139 A→T 17.8% intergenic (+3/‑72) 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,1390AT17.8% 360.8 / 66.3 146intergenic (+3/‑72)yejM/proLessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (55/65);  new base T (11/15);  total (66/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 = 5.52e-01

Only 100 of 162 total aligned reads displayed.

TACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAC  >  NC_000913/2286041‑2286238
                                                                                                  |                                                                                                   
tACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                    >  1:942209/1‑100 (MQ=255)
tACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                    >  1:819926/1‑100 (MQ=255)
 aCAACTTACGTGGTGAAAGAGTGAAAGATGAATAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGatta                                                                                                     <  2:3143058/98‑1 (MQ=255)
   aaCTTACGTGGTGAAAGAAAAAAAAAAAAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatta                                                                                                 <  2:1579424/100‑1 (MQ=255)
     cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGGCAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                               >  1:2814299/1‑100 (MQ=255)
     cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatta                                                                                                 >  1:2133411/1‑98 (MQ=255)
        aCGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAAt                                                                                            <  1:3405003/100‑1 (MQ=255)
         cGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                           >  1:409012/1‑100 (MQ=255)
          gtggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCa                                                                                          >  1:1504850/1‑100 (MQ=255)
          gtggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCa                                                                                          >  1:1505334/1‑100 (MQ=255)
           tggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                         <  2:2139196/100‑1 (MQ=255)
           tggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                         >  1:3568161/1‑100 (MQ=255)
             cgatctgagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGtt                                                                                       <  2:1700917/94‑1 (MQ=38)
               gaaagagtgacagaTGGAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCCCTAACTGATTAATTATAAATCAGTTAg                                                                                     >  2:724297/1‑100 (MQ=255)
               gaaagagtgaaagaTGAAAAGCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAg                                                                                     <  1:1765611/100‑1 (MQ=38)
               gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAg                                                                                     <  2:777062/100‑1 (MQ=38)
               gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaaagatcgga                                                                                     >  2:1636474/1‑92 (MQ=255)
               gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                     >  2:1933954/1‑100 (MQ=255)
                aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                    <  1:1543270/100‑1 (MQ=38)
                aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                    <  1:2792461/100‑1 (MQ=38)
                 aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                   <  2:43131/100‑1 (MQ=38)
                 aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCg                                                                                   <  1:2555119/100‑1 (MQ=255)
                  agagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGa                                                                                  <  2:3263166/100‑1 (MQ=255)
                  agagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGa                                                                                  >  1:595503/1‑100 (MQ=255)
                   gagtgaatgaTAAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACCGATTTATAATTAATCAGTTAGc                                                                                    <  2:1097025/97‑1 (MQ=37)
                   gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTACCTGATTTATAATTAATCAGTTAGCGat                                                                                 <  1:730255/100‑2 (MQ=38)
                   gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAgatcg                                                                                 >  1:1700925/1‑96 (MQ=38)
                   gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                                                 <  1:2894545/100‑1 (MQ=255)
                   gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATCGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                                                 <  1:2555661/100‑1 (MQ=255)
                       gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                                                             >  2:3575390/1‑100 (MQ=255)
                       gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                                                             >  1:2358599/1‑100 (MQ=255)
                       gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                                                             >  1:2286288/1‑100 (MQ=255)
                        aaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATc                                                                            >  1:2594421/1‑100 (MQ=255)
                         aagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCt                                                                           >  2:3001176/1‑100 (MQ=255)
                         aagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCt                                                                           <  2:308941/100‑1 (MQ=255)
                          acaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                   <  2:2809690/89‑1 (MQ=38)
                               aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                     <  1:139467/100‑1 (MQ=255)
                               aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                     <  1:309686/100‑1 (MQ=255)
                               aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                     <  1:2380202/100‑1 (MQ=255)
                                aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGCGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                    >  2:1993641/1‑100 (MQ=255)
                                aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                    >  1:2532144/1‑100 (MQ=255)
                                aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                    <  1:715885/100‑1 (MQ=255)
                                   aCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAAt                                                                 <  1:3085687/100‑1 (MQ=255)
                                    ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATc                                                                >  1:822013/1‑100 (MQ=255)
                                     cacaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                               >  1:1793011/1‑100 (MQ=255)
                                      acaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCgn                                                              <  1:166389/100‑2 (MQ=255)
                                      acaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCgg                                                              <  1:2265924/100‑1 (MQ=255)
                                       caGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgt                                                             >  1:93675/1‑100 (MQ=255)
                                         gTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtgt                                                           <  1:2856887/100‑1 (MQ=255)
                                           tAAGTTTGTTATTGCAAGTGCTGACAGATGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTgg                                                         <  1:2124247/100‑1 (MQ=255)
                                           tAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTgg                                                         <  1:2642149/100‑1 (MQ=255)
                                           tAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTgg                                                         <  2:2865765/100‑1 (MQ=255)
                                              gTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaag                                                      >  1:3610136/1‑99 (MQ=255)
                                              gTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaa                                                      >  1:3268474/1‑100 (MQ=255)
                                                ttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAACTCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAAc                                                    <  2:1990246/100‑1 (MQ=255)
                                                ttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAAc                                                    <  2:3207583/100‑1 (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:528656/100‑1 (MQ=255)
                                                   ttATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAAACGGTGTGGAAAACGgt                                                 <  2:2286272/100‑1 (MQ=255)
                                                    tATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgta                                                <  1:1310319/100‑1 (MQ=255)
                                                     aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                               <  2:2080677/100‑1 (MQ=255)
                                                        gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                            >  2:3317565/1‑100 (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)
                                                               cTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagcaga                                     >  1:3624078/1‑99 (MQ=255)
                                                                 gACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGGGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCa                                   <  2:2737145/100‑1 (MQ=255)
                                                                 gACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCa                                   >  1:2023634/1‑100 (MQ=255)
                                                                   cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                 <  2:3584529/100‑1 (MQ=255)
                                                                    aGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                 <  2:1079526/99‑1 (MQ=255)
                                                                     gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                               >  2:1864717/1‑100 (MQ=255)
                                                                          gaAGCGTTTTATCGCTACCTGATTAATTATAAATCAGTTAGGGAAATATCTTACGTGCAATCGGTGGGGAAAACGGTAGTATTAGCAGCCACGAGTCGGc                          <  2:117960/100‑1 (MQ=255)
                                                                          gaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGc                          >  1:2962774/1‑100 (MQ=255)
                                                                           aaGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCa                         <  1:916965/100‑1 (MQ=255)
                                                                           aaGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCa                         <  1:274954/100‑1 (MQ=255)
                                                                              cGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGt                      <  2:1766085/100‑1 (MQ=255)
                                                                              agatatttCGCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGt                      >  2:916702/8‑100 (MQ=25)
                                                                                  ttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                  <  1:2843768/100‑1 (MQ=255)
                                                                                   tATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcgc                 >  2:1107588/1‑100 (MQ=255)
                                                                                    ttCGCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCa                >  1:539/2‑100 (MQ=38)
                                                                                    aTCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCa                <  1:2174379/100‑1 (MQ=255)
                                                                                      cGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg               >  2:2356586/1‑99 (MQ=255)
                                                                                       gCTAACTGGTTTATAATTAATCAGTTAGCGAAGTTTCTTACTGGCATTGGGTTTTGAATACGGTAGTATTAGCAGCGACGAGTCGTCAGGTTgc                   >  2:2172135/1‑94 (MQ=12)
                                                                                       gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCTATCGGAGTGGAAAACGGAAGTATTAGCAGCCACGAGACGGCACGTAGCGCAGc              >  2:1543682/1‑99 (MQ=37)
                                                                                       gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc             >  2:1345276/1‑100 (MQ=38)
                                                                                       gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc             >  1:1370825/1‑100 (MQ=38)
                                                                                       gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc             >  1:1669327/1‑100 (MQ=255)
                                                                                       gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc             <  2:161412/100‑1 (MQ=255)
                                                                                       gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGc              >  2:3499526/1‑99 (MQ=255)
                                                                                        cTAACTGATTTATAATTAATAAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt            >  2:2223204/1‑100 (MQ=38)
                                                                                          aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg               >  2:3045265/1‑95 (MQ=255)
                                                                                          aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTg           <  2:372884/99‑1 (MQ=255)
                                                                                              gATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgc      <  1:273569/100‑1 (MQ=255)
                                                                                               attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcg     <  1:1331669/100‑1 (MQ=255)
                                                                                                ttaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcgc    >  2:2612834/1‑100 (MQ=255)
                                                                                                  aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc  >  1:1434361/1‑100 (MQ=255)
                                                                                                  aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc  <  1:191358/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAC  >  NC_000913/2286041‑2286238

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: