Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,483,906 T→A 15.5% I51N (ATT→AAT)  evgA → response regulator in two‑component regulatory system with EvgS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,9060TA15.5% 305.9 / 44.3 116I51N (ATT→AAT) evgAresponse regulator in two‑component regulatory system with EvgS
Reads supporting (aligned to +/- strand):  ref base T (46/52);  new base A (10/8);  total (56/60)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

Only 100 of 123 total aligned reads displayed.

TTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATG  >  NC_000913/2483809‑2484004
                                                                                                 |                                                                                                  
ttATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTg                                                                                                  <  1:1991248/100‑1 (MQ=255)
ttATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTg                                                                                                  <  2:493741/100‑1 (MQ=255)
 tATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGa                                                                                                 <  1:926531/100‑1 (MQ=255)
 tATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGa                                                                                                 <  1:2230492/100‑1 (MQ=255)
     gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                             <  2:2682531/100‑1 (MQ=255)
     gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                             <  2:1918629/100‑1 (MQ=255)
     gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                             >  2:696967/1‑100 (MQ=255)
      aTCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCg                                                                                            >  1:482322/1‑100 (MQ=255)
        cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGACACACTTAAGCCTGATATCGTCATCATTGCTGTCGat                                                                                          >  2:2727318/1‑100 (MQ=255)
        cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                          <  2:960032/100‑1 (MQ=255)
        cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGat                                                                                          <  2:960842/100‑1 (MQ=255)
         aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGata                                                                                         >  1:2373081/1‑100 (MQ=255)
         aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGata                                                                                         <  1:1809516/100‑1 (MQ=255)
          ttccgatctttGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGatat                                                                                        <  2:1025690/91‑1 (MQ=255)
          aaaaCGATATTGACATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGatat                                                                                        >  2:1955342/1‑100 (MQ=255)
          aaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGatat                                                                                        >  1:1518137/1‑100 (MQ=255)
             aCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcag                                                                                     <  1:2309980/100‑3 (MQ=38)
              cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcccc                                                                                    <  1:2200048/100‑1 (MQ=255)
              cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                    <  1:57996/100‑4 (MQ=38)
              cGATATTGAAATCTTAGCAGAGTTGACGGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                    <  1:2870916/100‑4 (MQ=38)
               gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                   <  1:2150604/100‑5 (MQ=38)
               gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                   <  1:1233899/100‑5 (MQ=38)
               gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                   <  2:1395391/100‑5 (MQ=38)
                atatTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCgg                                                                                  <  2:2532822/100‑1 (MQ=255)
                  atTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAg                                                                                <  2:2394615/100‑1 (MQ=255)
                   ttGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGTTGTCGATATccc                                                                                     >  1:349187/1‑94 (MQ=255)
                   ttGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGatatagatcggaa                                                                               >  1:1025688/1‑91 (MQ=255)
                   ttGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGt                                                                               <  1:825025/100‑1 (MQ=255)
                     gAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTa                                                                             >  1:1976684/1‑100 (MQ=255)
                       aaTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAc                                                                           <  1:2412386/100‑1 (MQ=255)
                        aTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                          <  1:2667838/100‑1 (MQ=255)
                         tCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                         <  2:2445911/100‑1 (MQ=255)
                         tCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                         <  2:1036988/100‑1 (MQ=255)
                          cTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGt                                                                        >  2:1433448/1‑100 (MQ=255)
                              gCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATcc                                                                    >  2:2644866/1‑100 (MQ=255)
                                   gTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtg                                                               <  1:614720/100‑1 (MQ=255)
                                     tGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGtt                                                             <  1:2590582/100‑1 (MQ=255)
                                      gACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTa                                                            >  1:1724411/1‑100 (MQ=255)
                                      gACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTa                                                            >  2:2685148/1‑100 (MQ=255)
                                      atcTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTa                                                            <  1:1643161/98‑1 (MQ=255)
                                        cTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                          <  2:2495499/100‑1 (MQ=255)
                                          gAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTaagat                                                        >  2:1643165/1‑96 (MQ=255)
                                          gAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGaaa                                                        >  2:2463588/1‑100 (MQ=255)
                                          gAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGaaa                                                        >  1:982759/1‑100 (MQ=255)
                                            aGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACg                                                      >  2:1742531/1‑100 (MQ=255)
                                             ggCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGc                                                     >  2:1270197/1‑100 (MQ=255)
                                             ggCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGc                                                     >  1:1117759/1‑100 (MQ=255)
                                             ggCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGc                                                     <  2:313653/100‑1 (MQ=255)
                                               cGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTg                                                   <  2:2748820/100‑1 (MQ=255)
                                               cGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTg                                                   >  2:2774803/1‑100 (MQ=255)
                                                 gAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                 <  1:2401391/100‑1 (MQ=255)
                                                 gAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                 <  1:1516868/100‑1 (MQ=255)
                                                  atgTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                 <  1:806939/97‑1 (MQ=255)
                                                  aaGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAgg                                                <  2:1749119/100‑1 (MQ=255)
                                                         gTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGcc                                         <  1:2206110/100‑1 (MQ=255)
                                                         gTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGcc                                         <  1:301405/100‑1 (MQ=255)
                                                            cAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAat                                      <  1:141246/100‑1 (MQ=255)
                                                            cAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAat                                      <  2:924297/100‑1 (MQ=255)
                                                            cAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAat                                      <  2:1373906/100‑1 (MQ=255)
                                                             aGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAata                                     >  2:1036428/1‑100 (MQ=255)
                                                              gCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAatat                                    >  1:2749000/1‑100 (MQ=255)
                                                               cGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAAtata                                   <  1:1245387/100‑1 (MQ=255)
                                                                 ggTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGc                                 >  2:1063915/1‑100 (MQ=255)
                                                                   tGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCgg                               <  2:809364/100‑1 (MQ=255)
                                                                   tGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCgg                               <  2:217817/100‑1 (MQ=255)
                                                                     gAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGa                              >  2:2464354/1‑99 (MQ=255)
                                                                      aaaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAat                            >  1:967447/1‑100 (MQ=255)
                                                                         cacTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattat                         <  2:1599941/100‑1 (MQ=255)
                                                                         cacTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattat                         <  2:1599039/100‑1 (MQ=255)
                                                                          acTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatt                        >  1:1504897/1‑100 (MQ=255)
                                                                           cTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                       >  2:2020240/1‑100 (MQ=255)
                                                                           cTTAAGCCTGATATCATCATCATTGATATCTTTATCCCCGGAGTTAACTGTATCAAGGTGTTATAAACTCTCAGGAAGCGCCAATATAGCGGAattatta                       <  1:2641636/100‑1 (MQ=255)
                                                                             tAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattattatt                     >  1:952296/1‑100 (MQ=255)
                                                                             tAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattattatt                     >  1:1727472/1‑100 (MQ=255)
                                                                             tAAGACTGATATCGTCATCATTGATGTCGATATCACCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattattatt                     >  1:520137/1‑100 (MQ=255)
                                                                             ctccggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGCAGCGCCAATATAGCGGAattattatt                     >  1:1432513/8‑100 (MQ=25)
                                                                             ctccggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGAATCCAGGTGTTAGAAACGCTGAGGAAACGCCAATATAGCGGAattatta                       >  1:1432113/8‑98 (MQ=25)
                                                                               aGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATtattattat                   <  1:741664/100‑1 (MQ=255)
                                                                                 ggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                 >  2:1825810/4‑100 (MQ=38)
                                                                                 ggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                 >  2:1037841/4‑100 (MQ=38)
                                                                                 ccTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                 <  1:1270195/100‑1 (MQ=255)
                                                                                 ccTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                 >  2:912312/1‑100 (MQ=255)
                                                                                  gggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGt                >  2:262580/3‑100 (MQ=38)
                                                                                   tGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtc               >  2:2137521/1‑100 (MQ=255)
                                                                                   tGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtc               <  2:735050/100‑1 (MQ=255)
                                                                                   ggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtc               >  2:2508496/2‑100 (MQ=38)
                                                                                    gATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtct              >  2:679982/1‑100 (MQ=39)
                                                                                      tatCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTcc            <  1:869643/100‑1 (MQ=255)
                                                                                      tatCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTcc            >  2:2231547/1‑100 (MQ=255)
                                                                                      tatCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTcc            >  1:2545002/1‑100 (MQ=255)
                                                                                        tCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGc          >  1:2552004/1‑100 (MQ=255)
                                                                                        tCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGc          <  2:2489937/100‑1 (MQ=255)
                                                                                         cGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCt         <  2:1434000/100‑1 (MQ=255)
                                                                                          gTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTa        >  2:1048890/1‑100 (MQ=255)
                                                                                            catcatTGATGTCGATATCCCCGGTGTTAACGGTATCCAGGTGTTAGAAACGCGGAGTAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTaaa      >  1:2040807/1‑100 (MQ=255)
                                                                                             atcatTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTaaaa     >  2:1181083/1‑100 (MQ=255)
                                                                                              tcatTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTaaaaa    <  1:996050/100‑1 (MQ=255)
                                                                                               catTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAAt   <  2:1105244/100‑1 (MQ=255)
                                                                                                atTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATg  <  2:1413860/100‑1 (MQ=255)
                                                                                                atTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATg  <  1:949180/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
TTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATG  >  NC_000913/2483809‑2484004

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