Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,556,769 (T)9→8 intergenic (‑449/‑224) yjiC ← / → iraD uncharacterized protein/RpoS stabilzer after DNA damage, anti‑RssB factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,556,7610T.93.9% 589.8 / 32.4 149intergenic (‑441/‑232)yjiC/iraDuncharacterized protein/RpoS stabilzer after DNA damage, anti‑RssB factor
Reads supporting (aligned to +/- strand):  ref base T (7/2);  new base . (48/91);  total (55/94)
Fisher's exact test for biased strand distribution p-value = 1.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 179 total aligned reads displayed.

TATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGA  >  NC_000913/4556663‑4556856
                                                                                                  |                                                                                               
  tGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttt                                                                                              >  1:680509/1‑100 (MQ=255)
  tGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttt                                                                                              >  2:320121/1‑100 (MQ=255)
  tGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttt                                                                                              >  2:1010127/1‑100 (MQ=255)
     ggAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttttn                                                                                            <  1:1530933/99‑2 (MQ=255)
      gAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttttttt                                                                                          <  2:647678/100‑1 (MQ=255)
      gAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttttttt                                                                                          >  2:932354/1‑100 (MQ=255)
       aaattatttatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttttg                                                                                         <  1:31880/91‑2 (MQ=255)
       aaataaatAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttttg                                                                                         <  2:256534/100‑2 (MQ=255)
        aataaataATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTATCGGttttttttga                                                                                        <  1:2072514/100‑3 (MQ=255)
         ataaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttttgag                                                                                       <  2:76880/100‑4 (MQ=255)
         ataaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttttgag                                                                                       <  1:937166/100‑4 (MQ=255)
           aaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGtt                                                                                    >  2:716435/1‑100 (MQ=255)
             ataatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttttttttttttt                                                                                    >  2:1360657/1‑94 (MQ=255)
             ataatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTtg                                                                                  <  2:589734/100‑1 (MQ=255)
              taatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTtga                                                                                 <  1:1099080/100‑1 (MQ=255)
               aatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTtgat                                                                                <  1:979266/100‑1 (MQ=255)
                   tCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTgatgat                                                                             <  2:2247774/99‑1 (MQ=255)
                     ccACCTGAATTAAGTAGCAGCGGTTACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGAtttt                                                                          <  2:2056893/100‑1 (MQ=255)
                     ccACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGAtttt                                                                          <  2:2333681/100‑1 (MQ=255)
                     ccACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGAtttn                                                                          <  1:658489/100‑2 (MQ=255)
                     ccACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACGATTCGGAATATTCAGAATTTACTAGCGGTTTTTTTTGAGTTTGATGAtttt                                                                          <  1:2929665/100‑1 (MQ=255)
                      cACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTAGTTTGATGAttttt                                                                         >  2:2111741/1‑100 (MQ=255)
                      cACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGAttttt                                                                         >  1:2642055/1‑100 (MQ=255)
                        ccGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTTTTTTATTATTTTTAt                                                                       >  2:2548037/1‑100 (MQ=255)
                        ccGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTAt                                                                       >  1:2356175/1‑100 (MQ=255)
                          ggAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTAAGTTATTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCa                                                                     <  2:2848888/100‑1 (MQ=255)
                           gAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCa                                                                     <  2:805162/99‑1 (MQ=255)
                           gAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCa                                                                     >  1:1775572/1‑99 (MQ=255)
                            aaTTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCaa                                                                    <  2:1547077/99‑1 (MQ=255)
                             aTTAAGTAGCAGCGGTGACTTGTATTTTCATCGACGTAATTGTCCGGAATATTCAGATTATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATc                                                                      <  2:2440230/96‑1 (MQ=255)
                             aTTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAAt                                                                  <  1:897354/100‑1 (MQ=255)
                               tAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATGTTCAGAATGTACTAGCGGGTTGTTTTGAGTTTGATGAGTTTTATCaaa                                                                   >  1:1685261/1‑97 (MQ=255)
                               tAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGGTTTTTATCAAATCa                                                                <  2:1477438/100‑1 (MQ=255)
                               tAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCa                                                                <  2:1243678/100‑1 (MQ=255)
                               tAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCa                                                                >  2:1222782/1‑100 (MQ=255)
                                aaGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCAt                                                               <  2:2440313/100‑1 (MQ=255)
                                aaGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCAt                                                               >  1:1926634/1‑100 (MQ=255)
                                 aGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATg                                                              <  1:1870911/100‑1 (MQ=255)
                                   tAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGAc                                                            >  2:1395235/1‑100 (MQ=255)
                                     gcagcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACtt                                                          >  1:1888572/1‑100 (MQ=255)
                                         cGGTGACAAGTATTTTCATCTACGTATCTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACtttttt                                                      <  2:2583562/100‑1 (MQ=255)
                                         cGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACtttttt                                                      <  2:1546152/100‑1 (MQ=255)
                                          ggTGATAAGTATTTTTATCTACGTAATTATTCCGAATTTTCAGAATATACTCGCCATGTTTTTTGCGTTTCACGATTTTTATCAAATCATGACTTTTTTc                                                     <  2:385743/100‑1 (MQ=255)
                                           gTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTcc                                                    <  1:278892/100‑1 (MQ=255)
                                           gTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTcc                                                    <  1:2109861/100‑1 (MQ=255)
                                            tGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCt                                                   <  2:1646909/100‑1 (MQ=255)
                                            tGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCt                                                   <  1:1611464/100‑1 (MQ=255)
                                              aCAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTg                                                  <  2:679298/99‑1 (MQ=255)
                                               cAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTTTTTTTTTATTTTTATCAAATCATGACTTTTTTCCt                                                   >  2:2718923/1‑97 (MQ=255)
                                               cAAGTATTTTCATCTACGTAACGATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaa                                                >  2:1795323/1‑100 (MQ=255)
                                                aaGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaaa                                               <  1:2879792/100‑1 (MQ=255)
                                                aaGTATTTTCATATACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaa                                                <  1:1795342/99‑1 (MQ=255)
                                                 aGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaaaa                                              <  2:2664695/100‑1 (MQ=255)
                                                 aGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaaaa                                              >  1:2980495/1‑100 (MQ=255)
                                                  gTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAg                                             <  1:748846/100‑1 (MQ=255)
                                                  gTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAg                                             >  2:297037/1‑100 (MQ=255)
                                                    aTTTTAATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTc                                           <  1:1003769/100‑1 (MQ=255)
                                                     ttttTATCTATGTAATTATTTGTAATATTCAGAATATACTAGTGGTTTTTTTTGAGTTTGATGATTTTTATTAAATCATGATTTTTTTTCTTAAAAGtt                                           <  2:206827/99‑2 (MQ=255)
                                                     ttttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCa                                          <  2:1598564/100‑1 (MQ=255)
                                                        tCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTa                                       <  1:1265314/100‑1 (MQ=255)
                                                         cATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAg                                      >  1:1737145/1‑100 (MQ=255)
                                                         cATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAg                                      <  2:17859/100‑1 (MQ=255)
                                                         cATCTACGTAAATATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAg                                      <  1:2487088/100‑1 (MQ=255)
                                                          aTCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGa                                     <  1:2543473/100‑1 (MQ=255)
                                                             tACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGatat                                  >  1:544402/1‑100 (MQ=255)
                                                                gTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCa                               <  1:716438/100‑1 (MQ=255)
                                                                 tAATTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCa                               <  1:1105453/99‑1 (MQ=255)
                                                                  aaCTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATa                             >  1:81878/1‑100 (MQ=255)
                                                                    cTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAgg                           <  2:1926612/100‑1 (MQ=255)
                                                                    cTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAgg                           >  1:1369264/1‑100 (MQ=255)
                                                                     tATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGc                          <  1:1114984/100‑1 (MQ=255)
                                                                     tATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGc                          <  1:2532820/100‑1 (MQ=255)
                                                                       ttCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCaa                        >  1:2187648/1‑100 (MQ=255)
                                                                        tCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAg                       <  1:758896/100‑1 (MQ=255)
                                                                        tCGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAg                       <  2:558358/100‑1 (MQ=255)
                                                                         cGGAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTACTGAAAAGTAAGTAGATATTCATAGGCAAGt                      >  2:1074703/1‑100 (MQ=255)
                                                                          ggAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGt                      <  1:284126/99‑1 (MQ=255)
                                                                           gAATATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTaa                    <  2:1862815/100‑1 (MQ=255)
                                                                            aaTATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATTAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAg                   >  1:2772590/1‑100 (MQ=255)
                                                                            aaTATTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTaa                    <  2:1966026/99‑1 (MQ=255)
                                                                             atatTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAgg                  <  1:1642075/100‑1 (MQ=255)
                                                                              tatTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGt                 <  1:431076/100‑1 (MQ=255)
                                                                               atTCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGtt                <  1:2983623/100‑1 (MQ=255)
                                                                                ttCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGttt               <  1:536890/100‑1 (MQ=255)
                                                                                ttCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGttt               <  1:1078283/100‑1 (MQ=255)
                                                                                ttCAGAATATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGttn               <  1:1542921/100‑2 (MQ=255)
                                                                                    gAATATACTAGCGGTTTTTTTTTTTTTTTATGATTTTTATCAAATCATGACTTTTTTCCTGTAAAGTCAGTAGATATTCATAGGCATGTAAGGTTTtata           >  2:2483483/1‑100 (MQ=255)
                                                                                     aaTATACTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATAc          <  1:2779807/100‑1 (MQ=255)
                                                                                       tataCTAGCGGTTTTTTTTTTTTTTTTTTATTTTTATCAAAGCATGACTTTTTTCTTTGTGAGTTAGTTGATATTCATAGGCTTTTCCGGTTTTATACtt        >  2:349508/1‑100 (MQ=255)
                                                                                       tataCTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACtt        <  2:1852523/100‑1 (MQ=255)
                                                                                        ataCTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACttt       <  1:2144513/100‑1 (MQ=255)
                                                                                         taCTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTg      <  2:2565824/100‑1 (MQ=255)
                                                                                         taCTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTg      <  2:1521835/100‑1 (MQ=255)
                                                                                          aCTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGc     <  1:2420770/100‑1 (MQ=255)
                                                                                           cTAGCGGTTTTTTTTTTGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCt    >  2:95640/1‑100 (MQ=255)
                                                                                           cTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCt    >  2:39739/1‑100 (MQ=255)
                                                                                           cTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCt    >  1:2692710/1‑100 (MQ=255)
                                                                                           cTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCt    <  1:1501275/100‑1 (MQ=255)
                                                                                           cTAGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTACTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCt    <  1:2015324/100‑1 (MQ=255)
                                                                                             aGCGGTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGa  >  1:1055436/1‑100 (MQ=255)
                                                                                                  |                                                                                               
TATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGA  >  NC_000913/4556663‑4556856

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: