Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,533,536 A→G A83A (GCT→GCC)  yjhQ ← GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,533,5360AG97.9% 495.5 / ‑3.3 142A83A (GCT→GCC) yjhQGNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base G (74/64);  total (74/68)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

GTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCT  >  NC_000913/4533441‑4533629
                                                                                               |                                                                                             
gTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGa                                                                                           >  2:615869/1‑100 (MQ=255)
gTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGa                                                                                           <  2:2739337/100‑1 (MQ=255)
 tCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAg                                                                                          <  2:666795/100‑1 (MQ=255)
   tGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGa                                                                                        >  1:1201171/1‑100 (MQ=255)
   tGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGa                                                                                        >  2:2810926/1‑100 (MQ=255)
    ggCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAt                                                                                       >  2:374218/1‑100 (MQ=255)
    ggCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAt                                                                                       >  2:2827718/1‑100 (MQ=255)
      cAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAtgt                                                                                     <  2:3033639/100‑1 (MQ=255)
       aGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGAtgtg                                                                                    <  1:1518982/100‑1 (MQ=255)
        gCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGCTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGc                                                                                   <  1:2173548/100‑1 (MQ=255)
         cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCa                                                                                  <  1:381069/100‑1 (MQ=255)
         cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCa                                                                                  <  2:2805631/100‑1 (MQ=255)
         cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCa                                                                                  <  1:1091049/100‑1 (MQ=255)
          ccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGGCCGCCAGGGGGGCGAGGATGTGCAt                                                                                 >  2:1007319/1‑100 (MQ=255)
           naTCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATc                                                                                >  2:3061886/2‑100 (MQ=255)
           cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATc                                                                                <  1:753198/100‑1 (MQ=255)
           cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATc                                                                                >  1:2622085/1‑100 (MQ=255)
           cATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATc                                                                                >  1:705685/1‑100 (MQ=255)
              cAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAAc                                                                             <  2:1101664/100‑1 (MQ=255)
                 ccTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGt                                                                          >  2:1994347/1‑100 (MQ=255)
                 ccTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGt                                                                          <  1:2810953/100‑1 (MQ=255)
                  cTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTg                                                                         <  2:577321/100‑1 (MQ=255)
                  cTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTg                                                                         <  1:1163718/100‑1 (MQ=255)
                  cTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTg                                                                         <  2:1271265/100‑1 (MQ=255)
                  cTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTg                                                                         <  1:786167/100‑1 (MQ=255)
                     aGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCAGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAAt                                                                      >  2:2164211/1‑100 (MQ=255)
                      ggTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATc                                                                     >  1:2978756/1‑100 (MQ=255)
                       gTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATcc                                                                    <  1:878878/100‑1 (MQ=255)
                        tGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCa                                                                   >  1:533441/1‑100 (MQ=255)
                        tGCTCGATGCCAGTGCGTATCAACCGTCCCCCCACGCCCATCCCTCGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATcc                                                                    <  1:2115657/99‑1 (MQ=255)
                          cTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCAtt                                                                 >  1:2687434/1‑100 (MQ=255)
                          cTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCAtt                                                                 <  1:1375646/100‑1 (MQ=255)
                              aTGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTtctc                                                             >  1:938131/1‑100 (MQ=255)
                               tGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTcc                                                            <  1:1409805/100‑1 (MQ=255)
                               tGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGCATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTcc                                                            <  2:1695859/100‑1 (MQ=255)
                                 ccAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCtt                                                          >  2:633570/1‑100 (MQ=255)
                                 ccAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCtt                                                          >  2:940019/1‑100 (MQ=255)
                                  cAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCttt                                                         <  2:752989/100‑1 (MQ=255)
                                    gTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGa                                                       <  1:1692055/100‑1 (MQ=255)
                                       cGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAgg                                                    <  2:2151431/100‑1 (MQ=255)
                                        gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                                   >  2:2213058/1‑100 (MQ=255)
                                        gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                                   >  2:438692/1‑100 (MQ=255)
                                        gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                                   <  2:1508946/100‑1 (MQ=255)
                                        gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                                   >  1:2352244/1‑100 (MQ=255)
                                        gTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                                                   >  1:8982/1‑100 (MQ=255)
                                         tATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGtt                                                  >  2:14021/1‑100 (MQ=255)
                                            cAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGcc                                               >  2:369036/1‑100 (MQ=255)
                                             aaCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGccc                                              >  2:481138/1‑100 (MQ=255)
                                             aaCCGTCATCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGccc                                              >  1:1622019/1‑100 (MQ=255)
                                              aCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCg                                             <  2:938122/100‑1 (MQ=255)
                                              aCCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCg                                             >  2:2862355/1‑100 (MQ=255)
                                               ccGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCgg                                            <  2:157029/100‑1 (MQ=255)
                                                  tcctccCACGTCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta                                         <  2:2020152/100‑1 (MQ=255)
                                                  tcctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta                                         >  2:2669586/1‑100 (MQ=255)
                                                  tcctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta                                         <  1:1934734/100‑1 (MQ=255)
                                                  tcctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta                                         <  1:229106/100‑1 (MQ=255)
                                                    ctccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaa                                       >  1:1110035/1‑100 (MQ=255)
                                                     tccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaat                                      >  2:2391424/1‑100 (MQ=255)
                                                     tccCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaa                                       <  1:1998999/99‑1 (MQ=255)
                                                     tccCACGCCCATCCATTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaat                                      >  1:765732/1‑100 (MQ=255)
                                                      cccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaata                                     >  2:470931/1‑100 (MQ=255)
                                                       ccACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaataa                                    >  2:17751/1‑100 (MQ=255)
                                                         aCGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaa                                  <  1:1804229/100‑1 (MQ=255)
                                                          cGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTaaataaaat                                 >  2:2548273/1‑100 (MQ=255)
                                                           gCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTACCCGGGTaaataaaata                                <  2:3083918/100‑1 (MQ=255)
                                                             ccATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATg                              <  1:2827745/100‑1 (MQ=255)
                                                              cATCCCTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGc                             <  2:1520405/100‑1 (MQ=255)
                                                                tCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGGTGTACATCAACGGTGAATCCATTTCCCCTTTGAAGGTTGCCCGGGTAAATAAAATATGcc                            <  2:715590/99‑1 (MQ=255)
                                                                 cccTTGGTACTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCg                          >  2:2248611/1‑100 (MQ=255)
                                                                 cccTTGATATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCg                          >  1:1927899/1‑100 (MQ=255)
                                                                  ccTTGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGa                         >  1:2130532/1‑100 (MQ=255)
                                                                   cTTGGTATTCCGGTATGACCGCCAGGCGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGAc                        >  1:1885912/1‑100 (MQ=255)
                                                                     tGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCATCGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACgg                      <  2:1970282/100‑1 (MQ=255)
                                                                     tGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACgg                      <  1:1421644/100‑1 (MQ=255)
                                                                     tGGTATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACg                       <  2:2360557/99‑1 (MQ=255)
                                                                        tATTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCtt                   <  2:1512962/100‑1 (MQ=255)
                                                                         aTTCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCttt                  >  1:5676/1‑100 (MQ=255)
                                                                          ttCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTg                 <  2:2284970/100‑1 (MQ=255)
                                                                           tCCGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGc                >  1:2781665/1‑100 (MQ=255)
                                                                            ccGGTATGACCGCTAGGGGGGCGAGGATGTGCATCAACAGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGcc               >  2:1900510/1‑100 (MQ=255)
                                                                            ccGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGcc               >  2:266478/1‑100 (MQ=255)
                                                                             cGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCt              >  1:412604/1‑100 (MQ=255)
                                                                             cGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCt              <  2:1261828/100‑1 (MQ=255)
                                                                             cGGTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGATTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCt              >  2:1688586/1‑100 (MQ=255)
                                                                              ggTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCtt             >  1:1557657/1‑100 (MQ=255)
                                                                               gTATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTc            <  2:1470235/100‑1 (MQ=255)
                                                                                tATGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCa           <  1:590245/100‑1 (MQ=255)
                                                                                  tGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATa         <  1:2248634/100‑1 (MQ=255)
                                                                                  tGACCGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATa         >  2:859434/1‑100 (MQ=255)
                                                                                   gACCGCCAGGGGGGAGAGTATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCAt          <  2:2504944/98‑1 (MQ=255)
                                                                                     ccGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCg      >  2:1640920/1‑100 (MQ=255)
                                                                                     ccGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCg      >  1:2644958/1‑100 (MQ=255)
                                                                                     ccGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCg      >  2:1092207/1‑100 (MQ=255)
                                                                                      cGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATGCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCgg     <  2:876253/100‑1 (MQ=255)
                                                                                      cGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATAccag     >  2:1878328/1‑98 (MQ=255)
                                                                                      cGCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCgg     <  1:1893438/100‑1 (MQ=255)
                                                                                       gCCAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCggg    >  2:512108/1‑100 (MQ=255)
                                                                                        ccAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGc   >  2:1760375/1‑100 (MQ=255)
                                                                                        ccAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGc   <  1:1007330/100‑1 (MQ=255)
                                                                                         cAGGGGGGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCt  <  2:1694702/100‑1 (MQ=255)
                                                                                               |                                                                                             
GTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAAATATGCCCGACGGCTTTGCCTTCATACCGGGCT  >  NC_000913/4533441‑4533629

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