Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,509,246 A→G 15.2% intergenic (‑182/‑40) yncJ ← / → hicA uncharacterized protein/mRNA interferase toxin of the HicAB toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,509,2460AG15.2% 214.1 / 32.6 99intergenic (‑182/‑40)yncJ/hicAuncharacterized protein/mRNA interferase toxin of the HicAB toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (38/46);  new base G (10/5);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

Only 100 of 103 total aligned reads displayed.

AGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGC  >  NC_000913/1509147‑1509345
                                                                                                   |                                                                                                   
aGAAATTCTTGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCGATTTTCGAGCAAATTACACACGGAGGTAAACACTaa                                                                                                     <  2:2802261/100‑1 (MQ=255)
 gAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAAt                                                                                                    <  2:950704/100‑1 (MQ=255)
 gAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAAt                                                                                                    >  2:710861/1‑100 (MQ=255)
 gAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGAAGGTAAAAACaaa                                                                                                     >  1:609707/1‑99 (MQ=255)
  aaaTTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATg                                                                                                   <  1:1072961/100‑1 (MQ=255)
    aTTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGTAAATTACACACGGAGGTAAACACTAATGtt                                                                                                 <  2:60717/100‑1 (MQ=255)
    aTTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGtt                                                                                                 >  1:1570045/1‑100 (MQ=255)
      tCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTa                                                                                               <  1:2006543/100‑1 (MQ=255)
       cTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGAAGGTAAACACTAATGTTTAc                                                                                              >  2:338568/1‑100 (MQ=255)
         ggAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACtt                                                                                            >  2:2874405/1‑100 (MQ=255)
          gAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACATTAGTGTTTAcctc                                                                                           <  1:810184/100‑4 (MQ=255)
          gAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACATTAGTGTTTAcctc                                                                                           <  1:2686893/100‑4 (MQ=255)
           aaTCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACATTAGTGTTTAcctcc                                                                                          <  2:371287/100‑5 (MQ=255)
           aaTCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttg                                                                                          >  1:912126/1‑100 (MQ=255)
            aTCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTCTTTTCGAGCAAATTACACACGGAGGTAAACATTAGTGTTTAcctcc                                                                                          <  2:1674959/99‑5 (MQ=255)
            aTCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttgt                                                                                         >  2:45763/1‑100 (MQ=255)
            aTCTTCCTTCATGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGGAAATTACACACGGAAGTAAACATTAGTGTGTAcctcc                                                                                          <  2:451140/99‑5 (MQ=255)
             tCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttgtt                                                                                        <  1:2727879/100‑1 (MQ=255)
               ttccttccTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGa                                                                                      <  2:912120/100‑1 (MQ=255)
               ttccttccTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGa                                                                                      <  2:1824622/100‑1 (MQ=255)
                tccttcctGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGat                                                                                     <  2:1253358/100‑1 (MQ=255)
                 ccttcctGATTTTGCGTTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGata                                                                                    <  1:2997591/100‑1 (MQ=255)
                  cttcctGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGatat                                                                                   <  1:1595675/100‑1 (MQ=255)
                  cttcctGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGatat                                                                                   >  1:2360318/1‑100 (MQ=255)
                  cttcctGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGatat                                                                                   <  2:732089/100‑1 (MQ=255)
                  cttcctGATTTTACATTGCAATATGACGTAGCGGCGATTTAGTGATATGATCGAGACAATTACACACTGAGGTAAACACTAATGGTTAGTTTGGTGata                                                                                    >  2:350383/1‑99 (MQ=255)
                   ttcctGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGAtata                                                                                  <  1:91240/100‑1 (MQ=255)
                    tcctGATTTTGCATTGCATTCTGCCGTTGCGGCTATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATAc                                                                                 <  1:1538147/100‑1 (MQ=255)
                       tGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCa                                                                              <  1:1532214/100‑1 (MQ=255)
                              gcattgcattCTGCCGTTGCGGTGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAg                                                                       <  1:712767/100‑1 (MQ=255)
                                  tgcattCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAg                                                                   <  1:549061/100‑1 (MQ=255)
                                  tgcattCTGCCGTTGCGGCGATTTAGTGCGATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAg                                                                   <  1:2255713/100‑1 (MQ=255)
                                  tgcattCTGCCGTTGCGGCGATTTAGAGCTATTTTCGAGCAAATTACACACGGAGGTAAAAACTAATGTTTACTTTGTTGATATACTAAGCGGCAGGGAg                                                                   >  1:1847505/1‑100 (MQ=255)
                                   gcattCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGCGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGa                                                                    >  1:1752487/1‑98 (MQ=255)
                                   gcattCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGa                                                                    >  2:2604339/1‑98 (MQ=255)
                                    cattCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGc                                                                 >  2:1785949/1‑100 (MQ=255)
                                        cTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAAATACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGAtg                                                             <  2:2067608/100‑1 (MQ=255)
                                         tGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACATACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGAtgt                                                            <  1:327198/100‑1 (MQ=255)
                                         tGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGTGGCAGGGAGGCGAtgt                                                            >  1:546939/1‑100 (MQ=255)
                                          gCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGAtgtg                                                           <  2:1549523/100‑1 (MQ=255)
                                          gCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGAtgtg                                                           <  2:478393/100‑1 (MQ=255)
                                            cGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACATGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGaa                                                         <  1:246717/100‑1 (MQ=255)
                                            cGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGaa                                                         >  2:1745442/1‑100 (MQ=255)
                                              ttGCGGCGATTTAGTGCTATTTTCGAGTAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAAc                                                       <  1:1070970/100‑1 (MQ=255)
                                              ttGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGTAGGCGATGTGAAAc                                                       <  2:1009329/100‑1 (MQ=255)
                                              ttGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAAc                                                       <  2:522506/100‑1 (MQ=255)
                                               tgcgtggagttAGTGCTATTTTCGAGCAAATTACATACGGAGGTAAACACTAATGTTTACTTTGTTGATATAATCAGCTGTAGGGAGGCGATGTGAAACa                                                      <  2:2810235/91‑1 (MQ=255)
                                               tGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACa                                                      <  1:1927568/100‑1 (MQ=255)
                                               tGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACa                                                      <  1:3012942/100‑1 (MQ=255)
                                                gcggcgATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACaa                                                     >  2:1181672/1‑100 (MQ=255)
                                                gcggcgATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACaa                                                     >  1:1927935/1‑100 (MQ=255)
                                                  ggttaTTTATTGCTATTTTCGATTAAATTACACACTTATGTAAACACTAATGTGTACTTTGTGGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAg                                                   <  1:555388/96‑1 (MQ=255)
                                                   gcgATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGc                                                  <  1:1213442/100‑1 (MQ=255)
                                                    cgATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCg                                                 <  2:2293900/100‑1 (MQ=255)
                                                    cgATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCg                                                 <  2:3000149/100‑1 (MQ=255)
                                                      aTTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAg                                               <  2:1013407/100‑1 (MQ=255)
                                                       tttAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGt                                              >  2:2874485/1‑100 (MQ=255)
                                                           gTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAg                                          >  2:2279804/1‑100 (MQ=255)
                                                              cTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACg                                       <  2:2101562/100‑1 (MQ=255)
                                                                 ttttCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTg                                    <  2:2401606/100‑1 (MQ=255)
                                                                   ttCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGc                                  <  1:2206396/100‑1 (MQ=255)
                                                                   ttCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAAACGt                                      >  1:2203216/1‑96 (MQ=255)
                                                                    tCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCt                                 <  2:3020664/100‑1 (MQ=255)
                                                                       aGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGa                              >  1:2004922/1‑100 (MQ=255)
                                                                        gCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGaa                             >  2:410675/1‑100 (MQ=255)
                                                                        gCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCAAGTTCAAACGTTGGCTCGaa                             >  1:1284428/1‑100 (MQ=255)
                                                                        gCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAACGAGAACAGACGTTGGCTCGaa                             >  1:1774713/1‑100 (MQ=255)
                                                                          aaaTTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAAtc                           >  2:2958965/1‑100 (MQ=255)
                                                                          aaaTTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAAtc                           <  1:2279830/100‑1 (MQ=255)
                                                                            aTTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAAtctc                         <  2:2365082/100‑1 (MQ=255)
                                                                            aTTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAAtctc                         >  1:2242047/1‑100 (MQ=255)
                                                                               agacacGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAgg                      <  2:2852477/98‑1 (MQ=255)
                                                                                tccgatctggTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAggg                     <  1:98610/92‑1 (MQ=255)
                                                                                 tcaacaaagTAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAggg                     >  2:2652578/9‑99 (MQ=38)
                                                                                  cacGGATGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGc                    >  2:1206995/1‑99 (MQ=255)
                                                                                  cacGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCg                   >  1:649613/1‑100 (MQ=255)
                                                                                  cacGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGAGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCATACGTTGGCTCGAATCTCAggg                     >  2:1724483/1‑98 (MQ=255)
                                                                                  caacaaagTAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCg                   >  2:917471/8‑100 (MQ=255)
                                                                                   aacaaagTAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCg                   >  2:1334476/7‑99 (MQ=255)
                                                                                   aacaaagTAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGt                  >  1:2350714/7‑100 (MQ=255)
                                                                                   aacaaagTAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGt                  >  1:212881/7‑100 (MQ=255)
                                                                                    cGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTAAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTc                 >  1:833352/1‑100 (MQ=255)
                                                                                    acaaagTAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTc                 >  1:2451581/6‑100 (MQ=255)
                                                                                     ggAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGtgg                <  2:475059/100‑3 (MQ=255)
                                                                                     caaagTAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCg                >  2:581313/5‑100 (MQ=255)
                                                                                      aaagTAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGAGTCGa               >  1:1367449/4‑100 (MQ=255)
                                                                                      aaagTAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGAGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGa               >  2:1870593/4‑100 (MQ=255)
                                                                                        ggTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATg             >  1:2762494/1‑100 (MQ=255)
                                                                                        ggTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGATCAGACGTTGGCTCGAATCTCAGGGCGTCGATg             >  1:2813387/1‑100 (MQ=255)
                                                                                         gTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGt            <  1:1778793/100‑1 (MQ=255)
                                                                                          tAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAACGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTa           >  1:1902518/1‑100 (MQ=255)
                                                                                            aaCACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGc         >  2:1724212/1‑100 (MQ=255)
                                                                                            aaCACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGc         >  2:2032289/1‑100 (MQ=255)
                                                                                             acacTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCg        <  2:984997/100‑1 (MQ=255)
                                                                                                cTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAAt     >  1:363393/1‑100 (MQ=255)
                                                                                                cTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTAAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAAt     >  1:865172/1‑100 (MQ=255)
                                                                                                cTAATGTTTACTTTGTTGATATACACAGCGACAGGGAGGCGATGTGAAACAAAGCGAGTTCAAACGTTGGCTAGAATCTCAGGGCGTCGATGTAGCGAAt     >  2:454493/1‑100 (MQ=255)
                                                                                                 tAATGTTTACTTTGTTGATATACTCAGCGGGAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAAt     <  1:25096/99‑1 (MQ=255)
                                                                                                  aaTGTTTACTTTGTTGATATAATCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTGAGACGTTGGCTCGAATCTCAGGGCGTCGATGAAGCGAATgg   <  2:378028/100‑1 (MQ=255)
                                                                                                   aTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGgc  <  1:1483543/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGC  >  NC_000913/1509147‑1509345

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