Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,500 Δ1 bp 23.0% intergenic (‑176/+171) 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,5000A.23.0% 329.5 / 99.9 126intergenic (‑176/+171)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base A (60/37);  new base . (11/18);  total (71/55)
Fisher's exact test for biased strand distribution p-value = 3.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

Only 100 of 129 total aligned reads displayed.

GCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTG  >  NC_000913/4520403‑4520598
                                                                                                 |                                                                                                    
gCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGgaga                                                                                                      >  1:2337509/1‑98 (MQ=255)
 cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACa                                                                                                    >  2:1093683/1‑99 (MQ=255)
 cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGAGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGCGACAt                                                                                                   >  2:1093684/1‑100 (MQ=255)
 cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATAAAATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                   >  1:387858/1‑100 (MQ=255)
   tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACTATTTGATTGATCGCCCACTATTCCCCCATACTGTAAAGCCGGAGACCTGc                                                                                               >  1:1663526/1‑100 (MQ=255)
   tGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGc                                                                                               >  2:1618334/1‑100 (MQ=255)
      ttttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTcc                                                                                            <  1:2032966/100‑1 (MQ=255)
        ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTc                                                                                             >  1:1624095/1‑97 (MQ=255)
         tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTGCCCCATACTGTAAAGCCGGAGACATGCTCCgg                                                                                          >  2:1253934/1‑99 (MQ=255)
          ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCt                                                                                        >  2:2060066/1‑100 (MQ=255)
           taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCtt                                                                                       >  1:2139667/1‑100 (MQ=255)
           taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                       <  1:1247448/100‑1 (MQ=38)
           taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                       <  2:1514396/100‑1 (MQ=38)
            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                      <  2:1315128/100‑1 (MQ=38)
            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                      <  1:1547201/100‑1 (MQ=38)
            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                      <  2:1315169/100‑1 (MQ=38)
            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                      <  2:1696345/100‑1 (MQ=38)
            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtta                                                                                      >  2:834834/1‑99 (MQ=38)
            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCt                                                                                        <  2:320447/98‑1 (MQ=38)
            aCAGATAGTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCCGTTCACCCATACTGTATAGCCGGAGCATGGCTCCGGCttt                                                                                      <  1:1107263/100‑1 (MQ=255)
            aCAGAGATTTGATGTTTTTACGTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCt                                                                                        <  1:1515731/98‑1 (MQ=25)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                     <  1:802847/100‑1 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                     >  1:798626/1‑100 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTGTa                                                                                     >  2:2230294/1‑100 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTa                                                                                     <  2:1729297/100‑1 (MQ=38)
              aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTAc                                                                                    <  2:721165/100‑1 (MQ=38)
               gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACa                                                                                   <  2:1860585/100‑1 (MQ=255)
               gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACa                                                                                   <  2:1291536/100‑1 (MQ=38)
                atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAg                                                                                  <  1:1673858/100‑1 (MQ=255)
                atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                  <  2:2479533/100‑1 (MQ=38)
                  atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTa                                                                                >  2:1044551/1‑100 (MQ=255)
                   ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTAt                                                                               >  2:1374657/1‑100 (MQ=255)
                      tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTAc                                                                            <  2:2139681/100‑1 (MQ=255)
                      tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTAc                                                                            <  2:2350065/100‑1 (MQ=255)
                      tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTAc                                                                            <  2:734182/100‑1 (MQ=255)
                        tgcTTTTACATGCAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTAc                                                                            <  1:583770/96‑1 (MQ=255)
                         tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTga                                                                         >  2:2106305/1‑100 (MQ=255)
                          cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                        <  2:686799/100‑1 (MQ=255)
                           ttttACATGAAATAAACCAGCATACCATGTGAGCGATCGCTCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgaga                                                                       <  2:1509322/100‑1 (MQ=255)
                           ttttACATGAAATAAACCAGAATACCATTTTAGGGCTCACCCGCTGATCCCACATACAGTAAGGCCGGAGACAGACTCGGGCTTTACAGTCTTACt                                                                           >  1:2459652/1‑96 (MQ=255)
                              tACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATa                                                                    >  1:886273/1‑100 (MQ=255)
                               acattgaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATACTCCGGCTTTACAGTATTACTGAGAAt                                                                     <  1:2288105/92‑1 (MQ=255)
                                cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTGCGGCTTTACAGTATTACTGAGAATaca                                                                  >  1:2591706/1‑100 (MQ=255)
                                cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                  >  2:1504189/1‑100 (MQ=255)
                                cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                  >  2:2355535/1‑100 (MQ=255)
                                cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                  >  2:846937/1‑100 (MQ=255)
                                cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                  >  2:88748/1‑100 (MQ=255)
                                   gAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAccc                                                               >  1:912844/1‑100 (MQ=255)
                                    aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                              >  1:2482064/1‑100 (MQ=255)
                                     aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                             <  2:1249070/100‑1 (MQ=255)
                                         aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                         >  2:1856157/1‑100 (MQ=255)
                                           ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCTCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                         <  1:328811/98‑1 (MQ=255)
                                           ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                       >  1:1466885/1‑100 (MQ=255)
                                              gaaTACCATTTGAACGATCGCCCGCTGTTCCACCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATaa                                                    >  2:1746457/3‑100 (MQ=255)
                                                aTACCTTTTGAGCGATCGCCCGCTGTTCCTCCATACTGTAGAGCCGGAGACATTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                  <  1:2293035/100‑1 (MQ=255)
                                                  aCCAATTGAGCGATCGCCCGCTGTTCCTCCATACTGTAAAGCCGAAGACATGCTCCAGTTTAACAGTATTACTGAGAATACACACATCAAGCATAATc                                                  >  1:30127/1‑98 (MQ=255)
                                                     aTTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                             <  2:1699295/100‑1 (MQ=255)
                                                        tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                          <  1:471472/100‑1 (MQ=255)
                                                            cGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAata                                      <  1:1874593/100‑1 (MQ=255)
                                                            cGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAata                                      <  1:1948161/100‑1 (MQ=255)
                                                             gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                     >  1:700642/1‑100 (MQ=255)
                                                                 gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                 >  2:235906/1‑100 (MQ=255)
                                                                 gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                 >  1:2271482/1‑100 (MQ=255)
                                                                    cGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCt                              >  2:909387/1‑100 (MQ=255)
                                                                      cTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttt                            >  1:2388162/1‑100 (MQ=255)
                                                                        gTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                          <  2:532405/100‑1 (MQ=255)
                                                                        gTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                          >  2:1058464/1‑100 (MQ=255)
                                                                         ttCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATTAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                         <  1:1121709/100‑1 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                      >  1:1325299/1‑100 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                      <  2:2471330/100‑1 (MQ=255)
                                                                             tctaTACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                     <  2:552448/97‑1 (MQ=255)
                                                                             cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                     <  2:1886997/100‑1 (MQ=255)
                                                                             cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                     >  1:103455/1‑100 (MQ=255)
                                                                             cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                     >  2:2165709/1‑100 (MQ=255)
                                                                              ccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGc                    >  1:1308249/1‑100 (MQ=255)
                                                                                aTACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCgaga                  >  1:552446/1‑97 (MQ=255)
                                                                                 tACTGTAAAGCCGGAGCATATCACAGGCTTAACAGTATTACTGAGAATAAACCCATCACGCATACTCGAAGGTTCAATATCGACCCTTATTCATCGcta                  >  1:315106/1‑97 (MQ=19)
                                                                                  aCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                <  2:2279819/100‑1 (MQ=255)
                                                                                  aCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                >  2:1805823/1‑100 (MQ=255)
                                                                                   cTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg               <  1:88747/100‑1 (MQ=255)
                                                                                   cTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg               >  2:2210814/1‑100 (MQ=255)
                                                                                   cTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg               >  2:1465351/1‑100 (MQ=255)
                                                                                   cTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg               >  2:324495/1‑100 (MQ=38)
                                                                                    tGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACTCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa              >  1:1918287/1‑100 (MQ=37)
                                                                                     gTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt             >  2:418106/1‑100 (MQ=38)
                                                                                      tAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt            >  1:1889492/1‑100 (MQ=38)
                                                                                       aaaGCCGGAGACCTGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  1:1135316/1‑100 (MQ=255)
                                                                                       aaaGCCGGAGACATGCTCCGGCTTTACAGTATTAGTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt             >  1:647248/1‑98 (MQ=255)
                                                                                       aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  2:2222094/1‑100 (MQ=38)
                                                                                       aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  2:498261/1‑100 (MQ=38)
                                                                                       aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTAAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  2:1829254/1‑100 (MQ=38)
                                                                                        aaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg          <  2:2476931/100‑1 (MQ=255)
                                                                                         aGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATTGCCAGCGAtt            >  2:816081/1‑97 (MQ=38)
                                                                                          gCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGt         >  1:1800499/1‑99 (MQ=255)
                                                                                          gCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc        >  2:733284/1‑100 (MQ=255)
                                                                                          gCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc        >  2:347210/1‑100 (MQ=255)
                                                                                          gCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc        <  1:568722/100‑1 (MQ=255)
                                                                                           ccGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCa       <  1:1561150/100‑1 (MQ=255)
                                                                                            cGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTACAATATCGCCCCATTTTAATCGCCCACGATTCGtaa       >  2:860071/1‑97 (MQ=255)
                                                                                                gaCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg  >  1:395833/1‑100 (MQ=255)
                                                                                                 |                                                                                                    
GCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTG  >  NC_000913/4520403‑4520598

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: