Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,504:1 +T 21.4% intergenic (‑180/+167) yjhU ← / ← yjhF putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,520,5041.T21.4% 331.0 / 116.2 126intergenic (‑180/+167)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base . (60/39);  new base T (10/17);  total (70/56)
Fisher's exact test for biased strand distribution p-value = 4.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

Only 100 of 128 total aligned reads displayed.

TGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGA  >  NC_000913/4520406‑4520602
                                                                                                    |                                                                                                  
tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACTATTTGATTGATCGCCCACTATTCCCCCATACTGTAAAGCCGGAGACCT‑G‑c                                                                                                   >  1:1663526/1‑100 (MQ=255)
   ttttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTcc                                                                                                <  1:2032966/100‑1 (MQ=255)
   ttttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTcc                                                                                                <  1:974474/100‑1 (MQ=255)
     ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTc                                                                                                 >  1:1624095/1‑97 (MQ=255)
      tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTGCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCgg                                                                                              >  2:1253934/1‑99 (MQ=255)
       ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCt                                                                                            >  2:2060066/1‑100 (MQ=255)
        tattgATATTTTTTTCTTTTACGTGTTATAAACCAGCATTCCATTTGGGCGTTCGCCCGCTTTTCCTCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGttt                                                                                           <  1:2582992/96‑1 (MQ=16)
        taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCtt                                                                                           >  1:2139667/1‑100 (MQ=255)
        taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCtt                                                                                           <  1:1247448/100‑1 (MQ=38)
        taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCtt                                                                                           <  2:1514396/100‑1 (MQ=38)
        taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCtt                                                                                           <  1:834831/100‑1 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCttt                                                                                          <  1:1547201/100‑1 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCttt                                                                                          <  2:1315128/100‑1 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCttt                                                                                          <  2:1315169/100‑1 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCttt                                                                                          <  2:1696345/100‑1 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCtta                                                                                          >  2:834834/1‑99 (MQ=38)
         aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCt                                                                                            <  2:320447/98‑1 (MQ=38)
         aCAGATAGTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCCGTTCACCCATACTGTATAGCCGGAG‑CATGG‑CTCCGGCttt                                                                                          <  1:1107263/100‑1 (MQ=255)
         aCAGAGATTTGATGTTTTTACGTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCt                                                                                            <  1:1515731/98‑1 (MQ=25)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTa                                                                                         >  1:473178/1‑100 (MQ=255)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTa                                                                                         >  1:798626/1‑100 (MQ=255)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCtt                                                                                           <  1:1423236/98‑1 (MQ=38)
          cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCTTTa                                                                                         <  2:1729297/100‑1 (MQ=38)
           aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCTTTAc                                                                                        <  2:721165/100‑1 (MQ=38)
            gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACa                                                                                       >  2:1549185/1‑100 (MQ=255)
            gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACa                                                                                       <  2:1860585/100‑1 (MQ=255)
            gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCTTTACa                                                                                       <  2:1291536/100‑1 (MQ=38)
             atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAg                                                                                      <  1:1673858/100‑1 (MQ=255)
             atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCTTTACAg                                                                                      <  2:2479533/100‑1 (MQ=38)
             atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CAT‑GTCTCCGGCTTTACAg                                                                                      <  2:2411086/100‑1 (MQ=38)
               atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTa                                                                                    >  2:1044551/1‑100 (MQ=255)
                ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTAt                                                                                   >  2:1374657/1‑100 (MQ=255)
                   tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTAc                                                                                <  2:2139681/100‑1 (MQ=255)
                   tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTAc                                                                                <  2:2350065/100‑1 (MQ=255)
                      tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTga                                                                             >  2:2106305/1‑100 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTgag                                                                            >  2:1846771/1‑100 (MQ=255)
                       cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTgag                                                                            <  2:686799/100‑1 (MQ=255)
                        ttttACATGAAATAAACCAGCATACCATGTGAGCGATCGCTCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTgaga                                                                           <  2:1509322/100‑1 (MQ=255)
                        ttttACATGAAATAAACCAGAATACCATTTTAGGGCTCACCCGCTGATCCCACATACAGTAAGGCCGGAGACAG‑A‑CTCGGGCTTTACAGTCTTACt                                                                               >  1:2459652/1‑96 (MQ=255)
                           tACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATa                                                                        >  1:886273/1‑100 (MQ=255)
                            acattgaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑A‑CTCCGGCTTTACAGTATTACTGAGAAt                                                                         <  1:2288105/92‑1 (MQ=255)
                             cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTGCGGCTTTACAGTATTACTGAGAATaca                                                                      >  1:2591706/1‑100 (MQ=255)
                             cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                      >  2:1504189/1‑100 (MQ=255)
                             cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                      >  2:2355535/1‑100 (MQ=255)
                             cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                      >  2:846937/1‑100 (MQ=255)
                             cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                      >  2:88748/1‑100 (MQ=255)
                              aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCTCATACTGTGAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTAATGAGAATacac                                                                     >  2:1126488/1‑100 (MQ=255)
                                 aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                                  >  1:2482064/1‑100 (MQ=255)
                                  aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                                 <  2:1249070/100‑1 (MQ=255)
                                      aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                             >  2:1856157/1‑100 (MQ=255)
                                        ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCTCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                             <  1:328811/98‑1 (MQ=255)
                                           gaaTACCATTTGAACGATCGCCCGCTGTTCCACCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATaa                                                        >  2:1746457/3‑100 (MQ=255)
                                             aTACCTTTTGAGCGATCGCCCGCTGTTCCTCCATACTGTAGAGCCGGAGACAT‑T‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                      <  1:2293035/100‑1 (MQ=255)
                                                  aTTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                                 <  2:1699295/100‑1 (MQ=255)
                                                        gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATACTCGAAGGTTCa                                             >  1:336463/1‑98 (MQ=255)
                                                        gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                           <  2:473180/100‑1 (MQ=255)
                                                         cGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAata                                          <  1:1874593/100‑1 (MQ=255)
                                                         cGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAata                                          <  1:1948161/100‑1 (MQ=255)
                                                          gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                         >  1:700642/1‑100 (MQ=255)
                                                                 cGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCt                                  >  2:909387/1‑100 (MQ=255)
                                                                 cGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCt                                  >  1:1693221/1‑100 (MQ=255)
                                                                     gTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                              <  2:532405/100‑1 (MQ=255)
                                                                     gTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                              >  2:1058464/1‑100 (MQ=255)
                                                                      ttCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATTAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                             <  1:1121709/100‑1 (MQ=255)
                                                                         ccccATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                          >  1:1325299/1‑100 (MQ=255)
                                                                         ccccATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                          <  2:2471330/100‑1 (MQ=255)
                                                                         ccccATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                          <  2:1901622/100‑1 (MQ=255)
                                                                          tctaTACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                         <  2:552448/97‑1 (MQ=255)
                                                                          cccATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                         >  1:103455/1‑100 (MQ=255)
                                                                           ccATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGc                        >  1:1308249/1‑100 (MQ=255)
                                                                             aTACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCgaga                      >  1:552446/1‑97 (MQ=255)
                                                                              tACTGTAAAGCCGGAG‑CATAT‑CACAGGCTTAACAGTATTACTGAGAATAAACCCATCACGCATACTCGAAGGTTCAATATCGACCCTTATTCATCGcta                      >  1:315106/1‑97 (MQ=19)
                                                                               aCTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                    <  1:2366242/100‑1 (MQ=255)
                                                                               aCTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                    <  2:2279819/100‑1 (MQ=255)
                                                                               aCTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                    >  2:1805823/1‑100 (MQ=255)
                                                                                cTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                   >  2:1465351/1‑100 (MQ=255)
                                                                                cTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                   <  1:88747/100‑1 (MQ=255)
                                                                                cTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                   >  2:2210814/1‑100 (MQ=255)
                                                                                cTGTAAAGCCGGAG‑CAT‑GTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                   >  2:324495/1‑100 (MQ=38)
                                                                                 tGTAAAGCCGGAG‑CAT‑GTCTCCGGCTTTACAGTATTACTGAGAATACACTCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa                  >  1:1918287/1‑100 (MQ=37)
                                                                                   tAAAGCCGGAG‑CAT‑GTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt                >  1:1889492/1‑100 (MQ=38)
                                                                                    aaaGCCGGAGACCT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc               >  1:1135316/1‑100 (MQ=255)
                                                                                    aaaGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTAGTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                 >  1:647248/1‑98 (MQ=255)
                                                                                    aaaGCCGGAG‑CAT‑GTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc               >  2:2222094/1‑100 (MQ=38)
                                                                                    aaaGCCGGAG‑CAT‑GTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc               >  2:498261/1‑100 (MQ=38)
                                                                                    aaaGCCGGAG‑CAT‑GTCTCCGGCTTTACAGTATTACTAAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc               >  2:1829254/1‑100 (MQ=38)
                                                                                     aaGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg              <  2:2476931/100‑1 (MQ=255)
                                                                                     aaGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg              >  2:2193085/1‑100 (MQ=255)
                                                                                      aGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                           >  1:624992/1‑86 (MQ=255)
                                                                                      aGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg              <  2:2271497/99‑1 (MQ=255)
                                                                                      aGCCGGAG‑CAT‑GTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATTGCCAGCGAtt                >  2:816081/1‑97 (MQ=38)
                                                                                       gCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc            >  2:347210/1‑100 (MQ=255)
                                                                                       gCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc            >  2:733284/1‑100 (MQ=255)
                                                                                       gCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc            <  1:568722/100‑1 (MQ=255)
                                                                                        ccGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCa           <  1:1561150/100‑1 (MQ=255)
                                                                                         cGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc          >  2:789489/1‑100 (MQ=255)
                                                                                         cGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTACAATATCGCCCCATTTTAATCGCCCACGATTCGtaa           >  2:860071/1‑97 (MQ=255)
                                                                                            agaCAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGt       <  1:2230279/100‑1 (MQ=255)
                                                                                             gaCAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg      >  1:395833/1‑100 (MQ=255)
                                                                                               cAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa    >  2:1801337/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
TGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAT‑G‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGA  >  NC_000913/4520406‑4520602

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: ‑