Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,142 T→A 16.1% 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,1420TA16.1% 346.1 / 55.8 137intergenic (+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/61);  new base A (7/15);  total (61/76)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

Only 100 of 160 total aligned reads displayed.

CAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCG  >  NC_000913/2286043‑2286240
                                                                                                   |                                                                                                  
cAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                                                    >  1:1494969/1‑100 (MQ=255)
cAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                                                    >  2:4108199/1‑100 (MQ=255)
 aaCTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatta                                                                                                   <  2:338123/100‑1 (MQ=255)
 aaCTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatta                                                                                                   >  1:642126/1‑100 (MQ=255)
   cTTACGTGGTGAAAGAGTGAAAGATGAAAAGCCACTGTTAGGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtatt                                                                                                 <  2:2613951/100‑2 (MQ=255)
   cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 <  2:419059/100‑1 (MQ=255)
   cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 <  1:2038290/100‑1 (MQ=255)
   cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 >  2:2849720/1‑100 (MQ=255)
    ttACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaa                                                                                                <  1:3476405/100‑1 (MQ=255)
     tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                               >  1:1949763/1‑100 (MQ=255)
     tACGTCGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGATAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                               <  2:10171/100‑1 (MQ=255)
       cGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                             >  1:709466/1‑100 (MQ=255)
       cGTGGTGAAAGAGTGAAAGAAAAAAAACCACAGTTAAGTTTGTTATTGCAAGTACTAACAGACGAGAAACGTTTTATCGCTAACTGATTAATTATAAATc                                                                                             >  2:588069/1‑100 (MQ=255)
        gtggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                             <  1:3765244/99‑1 (MQ=255)
        gtggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCa                                                                                            >  2:2119938/1‑100 (MQ=255)
         tggtgAAAGAGTGAAAGATGGAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                           >  1:3508549/1‑100 (MQ=255)
          ggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGGAAGTGTTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                           <  2:3911130/99‑1 (MQ=255)
          ggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                                                          <  2:3746433/100‑1 (MQ=255)
            tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTa                                                                                        <  1:212516/100‑1 (MQ=38)
            tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATGAATTATAAATAAGTTa                                                                                        >  2:2419586/1‑100 (MQ=255)
             gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                       >  1:3942401/1‑100 (MQ=255)
              aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                      <  1:3985085/100‑1 (MQ=38)
              aaagagtgaaagaTGAAAAAACACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                      <  1:106508/100‑1 (MQ=38)
              aaaaagtgaaaaaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGCAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                             <  2:2729086/93‑1 (MQ=255)
               aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                     <  2:823996/100‑1 (MQ=38)
               aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCg                                                                                     <  1:655287/100‑1 (MQ=38)
               aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCg                                                                                     <  2:2806236/100‑1 (MQ=255)
               aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCGAACTGATTTATAATTAATCAGTTAGc                                                                                      <  2:682403/99‑1 (MQ=38)
                agagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCGa                                                                                    <  1:388366/100‑1 (MQ=38)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                                                   >  2:658384/1‑100 (MQ=255)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGaa                                                                                   >  2:662364/1‑100 (MQ=255)
                    tgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAata                                                                                <  1:2357372/100‑1 (MQ=255)
                     gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGGTAGCGaaaga                                                                                >  2:3242389/1‑97 (MQ=255)
                      aaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATc                                                                              <  2:3809997/100‑1 (MQ=255)
                       tctaTGAAAACCCACAGTTAAGTTTGTTATTGCAAGTGTTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCTAAATATCt                                                                             <  2:2828922/97‑1 (MQ=255)
                       aagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCt                                                                             >  1:984465/1‑100 (MQ=255)
                        agaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTTACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCtt                                                                            <  1:4027234/100‑1 (MQ=255)
                         gaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTa                                                                           <  1:3221538/100‑1 (MQ=255)
                          aTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGTTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATACATCAGTTAGCGAAatcactaga                                                                          >  1:2828938/1‑93 (MQ=255)
                           tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACt                                                                         >  2:3119598/1‑100 (MQ=255)
                           tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACt                                                                         >  1:1335607/1‑100 (MQ=255)
                             aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGATTAAAGACTAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                       <  2:3602238/100‑1 (MQ=255)
                              aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                      <  1:2936693/100‑1 (MQ=255)
                              aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                      <  2:3104027/100‑1 (MQ=255)
                              aaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGc                                                                      <  1:2012104/100‑1 (MQ=255)
                                    acaGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCgg                                                                <  2:2099694/100‑1 (MQ=255)
                                     caGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgt                                                               <  1:3967289/100‑1 (MQ=255)
                                      aGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtg                                                              <  2:1880344/100‑1 (MQ=255)
                                      aGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtg                                                              <  2:1879256/100‑1 (MQ=255)
                                      aGTTAAGTTTGTGATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtg                                                              <  2:1076718/100‑1 (MQ=255)
                                       gTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtgt                                                             <  1:2683575/100‑1 (MQ=255)
                                        ttAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGtgtg                                                            >  1:2753819/1‑100 (MQ=255)
                                        ttAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGtgtg                                                            <  1:1103571/100‑1 (MQ=255)
                                         tAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTgg                                                           <  1:1759044/100‑1 (MQ=255)
                                             tttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCTCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaaa                                                       >  2:3392828/1‑100 (MQ=255)
                                              ttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAAc                                                      <  2:1879108/100‑1 (MQ=255)
                                                gTTATTGCAAGTGCTGACAGACGAGCAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                    >  2:2661452/1‑100 (MQ=255)
                                                gTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                    >  2:3705462/1‑100 (MQ=255)
                                                gTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                    >  1:1866176/1‑100 (MQ=255)
                                                 ttATTGGAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                    <  2:2263181/99‑1 (MQ=255)
                                                   aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                                 <  1:106285/100‑1 (MQ=255)
                                                   aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                                 >  2:427041/1‑100 (MQ=255)
                                                    ttGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagt                                                <  1:1970727/100‑1 (MQ=255)
                                                     tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGTGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                               <  2:4001520/100‑1 (MQ=255)
                                                     tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                               >  1:1070093/1‑100 (MQ=255)
                                                       cAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAtt                                             <  1:2659503/100‑1 (MQ=255)
                                                           tGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagca                                         <  1:1351869/100‑1 (MQ=255)
                                                           tGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagca                                         >  1:2113935/1‑100 (MQ=255)
                                                              taaCAGACGAGAAACGTTTTATCGCTAACTGATTAATTATAACTCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGcc                                      >  1:2786629/3‑100 (MQ=255)
                                                              tGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGcc                                      >  2:1568002/1‑100 (MQ=255)
                                                               gACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCa                                     >  2:3375516/1‑100 (MQ=255)
                                                                 cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                   <  2:1578389/100‑1 (MQ=255)
                                                                 cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                   <  1:2576874/100‑1 (MQ=255)
                                                                   gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                 >  1:3298527/1‑100 (MQ=255)
                                                                     cGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTc                               >  2:1154496/1‑100 (MQ=255)
                                                                      gagaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCg                              >  2:1222289/1‑100 (MQ=255)
                                                                          aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCAc                          >  1:3144143/1‑100 (MQ=255)
                                                                                atttCGCTAACTGATTTATAATTAATCAGTTAGCTAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                    >  2:667355/4‑100 (MQ=37)
                                                                                 tATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcgc                   >  1:2129810/1‑100 (MQ=255)
                                                                                    cGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                 >  2:4040668/1‑99 (MQ=255)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  2:1507076/1‑100 (MQ=38)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  1:443271/1‑100 (MQ=38)
                                                                                     gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  1:1020812/1‑100 (MQ=255)
                                                                                     gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  1:72804/1‑100 (MQ=255)
                                                                                      cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCACTCGGTGCGGAACACGGTAGCACTAGCCGCCACGAGCCCGCACGTCGCgccgccc              <  2:2408781/100‑2 (MQ=255)
                                                                                      cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt              >  2:3206767/1‑100 (MQ=255)
                                                                                        aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCGACGAGTCGGCACGTAGCGCAGCCTgg            <  1:2521005/100‑1 (MQ=255)
                                                                                         aCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGt           <  2:259213/100‑1 (MQ=255)
                                                                                          cTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCAACGAGTCGGCACGTAGCGCAGCCTGGt           >  1:2231898/1‑99 (MQ=255)
                                                                                           tGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAg         >  2:3528620/1‑100 (MQ=255)
                                                                                            gATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgc        <  2:2915265/100‑1 (MQ=255)
                                                                                             attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcg       >  1:1263223/1‑100 (MQ=255)
                                                                                             attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcg       >  1:2435915/1‑100 (MQ=255)
                                                                                              ttaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcgc      >  1:1928422/1‑100 (MQ=255)
                                                                                              ttaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcgc      <  1:662370/100‑1 (MQ=255)
                                                                                                aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc    <  1:3044801/100‑1 (MQ=255)
                                                                                                aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc    >  1:1442870/1‑100 (MQ=255)
                                                                                                 agtataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGGTGAAAACGGTAGTATTAGCAGCCACGAGGCGGCACGCAGCGCAGCCTGGTATCGCAc    >  2:3503062/3‑99 (MQ=255)
                                                                                                  ttaTAAATCATTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCTGCACTTTGCGCAGCCTGGTAGCGCACCg  <  2:1389797/100‑1 (MQ=255)
                                                                                                  ttaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCg  >  1:3840966/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCG  >  NC_000913/2286043‑2286240

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: