Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,509,243 C→T 15.5% intergenic (‑179/‑43) 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,2430CT15.5% 223.3 / 29.2 97intergenic (‑179/‑43)yncJ/hicAuncharacterized protein/mRNA interferase toxin of the HicAB toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (36/46);  new base T (10/5);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

Only 100 of 103 total aligned reads displayed.

CAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAAT  >  NC_000913/1509144‑1509342
                                                                                                   |                                                                                                   
cAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAacac                                                                                                     >  1:2101391/1‑100 (MQ=255)
 aaaGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACt                                                                                                    <  1:1073915/100‑1 (MQ=255)
  aaGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTa                                                                                                   <  2:798069/100‑1 (MQ=255)
   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:2686893/100‑4 (MQ=255)
             gAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACATTAGTGTTTAcctc                                                                                        <  1:810184/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:1824622/100‑1 (MQ=255)
                  ttccttccTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGa                                                                                   <  2:912120/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)
                                                                                                   |                                                                                                   
CAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAAT  >  NC_000913/1509144‑1509342

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