Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,588,028 A→G 5.2% Y280H (TAC→CAC)  yjiA ← metal‑binding GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,588,0280AG 5.2% 308.1 / 4.0 96Y280H (TAC→CAC) yjiAmetal‑binding GTPase
Reads supporting (aligned to +/- strand):  ref base A (45/46);  new base G (3/2);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGTAATC  >  NC_000913/4587882‑4588172
                                                                                                                                                  |                                                                                                                                                
gCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTag                                                                                                                                                 <  3:56079/148‑1 (MQ=255)
 cTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTaga                                                                                                                                                <  3:312845/148‑1 (MQ=255)
  tGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTagag                                                                                                                                               >  5:162070/1‑148 (MQ=255)
    cgcCTTATCCGTCCCATAACCCCTTACTAGCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAgg                                                                                                                                              <  6:162070/147‑1 (MQ=255)
          aTCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTgg                                                                                                                                        <  7:423713/147‑1 (MQ=255)
              ggCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGcc                                                                                                                                   >  2:219566/1‑148 (MQ=255)
               gCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGccc                                                                                                                                  >  8:437730/1‑148 (MQ=255)
               gCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGccc                                                                                                                                  <  8:156215/148‑1 (MQ=255)
                cccATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAAGCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCt                                                                                                                                 <  5:150717/148‑1 (MQ=255)
                   aTAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGa                                                                                                                              <  7:398037/148‑1 (MQ=255)
                      aCCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAca                                                                                                                           <  7:319098/148‑1 (MQ=255)
                      aCCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAca                                                                                                                           <  5:199530/148‑1 (MQ=255)
                        cccTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAcag                                                                                                                          >  8:70276/1‑147 (MQ=255)
                         ccTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAcag                                                                                                                          >  3:198300/1‑146 (MQ=255)
                          cTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAcagcag                                                                                                                       >  4:170306/1‑148 (MQ=255)
                            tACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGc                                                                                                                     >  1:365706/1‑148 (MQ=255)
                            tACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGc                                                                                                                     >  7:175370/1‑148 (MQ=255)
                             aCTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAc                                                                                                                            >  3:87271/1‑140 (MQ=255)
                             aCTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAc                                                                                                                            <  4:87271/140‑1 (MQ=255)
                             aCTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAc                                                                                                                            >  5:342512/1‑140 (MQ=255)
                             aCTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAAc                                                                                                                            <  6:342512/140‑1 (MQ=255)
                             aCTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCg                                                                                                                    >  7:258488/1‑148 (MQ=255)
                              cTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCg                                                                                                                    <  8:449046/147‑1 (MQ=255)
                                               aCGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATACACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAAt                                                                                                  <  2:34622/148‑1 (MQ=255)
                                                  cAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTGGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCca                                                                                               >  5:60448/1‑148 (MQ=255)
                                                  cAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGCCGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCca                                                                                               >  2:361472/1‑148 (MQ=255)
                                                  cAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGt                                                                                                      >  2:414361/1‑141 (MQ=255)
                                                  cAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGt                                                                                                      <  1:414361/141‑1 (MQ=255)
                                                   aGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCcac                                                                                              <  7:133139/148‑1 (MQ=255)
                                                   aGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCcac                                                                                              <  8:258488/148‑1 (MQ=255)
                                                      ccGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGc                                                                                           <  2:365706/148‑1 (MQ=255)
                                                        gAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCAt                                                                                         >  3:393874/1‑148 (MQ=255)
                                                          aTCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGtt                                                                                                                 >  3:172064/1‑122 (MQ=255)
                                                          aTCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGtt                                                                                                                 <  4:172064/122‑1 (MQ=255)
                                                              ctTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCttt                                                                                   >  5:37964/1‑148 (MQ=255)
                                                              ctTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCttt                                                                                   >  4:5105/1‑148 (MQ=255)
                                                               ttCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCtt                                                                                    <  8:302328/146‑1 (MQ=255)
                                                                tcctccGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGACTGTCACAGTCGGCGATGTAGAGGCGCTGGACGCCATGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTg                                                                                  >  7:32593/1‑147 (MQ=255)
                                                                 cctccGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCttt                                                                                   <  1:100750/145‑1 (MQ=255)
                                                                  ctccGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTaa                                                                               <  7:437730/148‑1 (MQ=255)
                                                                    ccGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACg                                                                             >  4:285497/1‑148 (MQ=255)
                                                                    ccGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACg                                                                             >  2:12794/1‑148 (MQ=255)
                                                                     cGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGt                                                                            <  3:5105/148‑1 (MQ=255)
                                                                     cGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGt                                                                            <  1:219566/148‑1 (MQ=255)
                                                                      ggCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTa                                                                           <  3:285497/148‑1 (MQ=255)
                                                                       gCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAg                                                                          <  4:393874/148‑1 (MQ=255)
                                                                          gCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCttt                                                                                   <  2:220506/136‑1 (MQ=255)
                                                                          gCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCttt                                                                                   >  1:220506/1‑136 (MQ=255)
                                                                              gATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGtt                                                                     <  8:175370/146‑1 (MQ=255)
                                                                              gATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTa                                                                   >  5:431074/1‑148 (MQ=255)
                                                                               aTCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTAt                                                                  >  8:329955/1‑148 (MQ=255)
                                                                               aTCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTAt                                                                  >  1:267567/1‑148 (MQ=255)
                                                                                  ccGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATccg                                                               >  2:336766/1‑148 (MQ=255)
                                                                                     aTAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGt                                                                                                             <  8:49469/99‑1 (MQ=255)
                                                                                     aTAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGt                                                                                                             >  7:49469/1‑99 (MQ=255)
                                                                                     aTAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATccgccg                                                            >  8:333055/1‑148 (MQ=255)
                                                                                          cacCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACCGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGCCTcc                                                       >  1:354022/1‑148 (MQ=255)
                                                                                           acCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTc                                                        <  2:99404/146‑1 (MQ=255)
                                                                                           acCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTc                                                        >  1:99404/1‑146 (MQ=255)
                                                                                              aTCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCcag                                                     >  8:396726/1‑146 (MQ=255)
                                                                                                     tATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTa                                                                   <  5:55398/125‑1 (MQ=255)
                                                                                                     tATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTa                                                                   >  6:55398/1‑125 (MQ=255)
                                                                                                       tGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGTAGGtttt                                          <  6:37964/148‑1 (MQ=255)
                                                                                                        gCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGttt                                                                    <  8:399590/121‑1 (MQ=255)
                                                                                                        gCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGttt                                                                    >  7:399590/1‑121 (MQ=255)
                                                                                                          ggTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTGGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCa                                       >  8:208459/1‑148 (MQ=255)
                                                                                                          ggTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTGGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCa                                       >  8:202327/1‑148 (MQ=255)
                                                                                                             ttttCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGt                                                                                                      <  7:114296/82‑1 (MQ=255)
                                                                                                             ttttCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGt                                                                                                      >  8:114296/1‑82 (MQ=255)
                                                                                                             atttCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCcagcagcag                                               >  3:61341/2‑137 (MQ=255)
                                                                                                             atttCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCcagcagcag                                               <  4:61341/136‑1 (MQ=255)
                                                                                                                 cATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCTATCCACAGCTTCCCTTTGTAACTTTGCAGTTTATCCGCCGACTCCAGCACCAGGTTTTCCATCAcgc                                 <  2:312221/147‑1 (MQ=255)
                                                                                                                 cATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCAcgc                                 <  1:12794/147‑1 (MQ=255)
                                                                                                                    cGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTcc                                        >  6:126629/1‑137 (MQ=255)
                                                                                                                    cGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTcc                                        <  5:126629/137‑1 (MQ=255)
                                                                                                                     gCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGaa                            <  7:396726/148‑1 (MQ=255)
                                                                                                                        ccACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATccc                                                                                      <  2:215035/87‑1 (MQ=255)
                                                                                                                        ccACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATccc                                                                                      >  1:215035/1‑87 (MQ=255)
                                                                                                                           cGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCg                      <  6:431074/148‑1 (MQ=255)
                                                                                                                           cGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCg                      <  1:336766/148‑1 (MQ=255)
                                                                                                                             gCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGCCGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCAc                                   >  5:169327/1‑133 (MQ=255)
                                                                                                                             gCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCAc                                   <  6:169327/133‑1 (MQ=255)
                                                                                                                                 gTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGata                >  3:21068/1‑148 (MQ=255)
                                                                                                                                  tCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGat                 <  4:198300/146‑1 (MQ=255)
                                                                                                                                     cAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTg                   >  3:217760/1‑141 (MQ=255)
                                                                                                                                     cAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTg                   <  4:217760/141‑1 (MQ=255)
                                                                                                                                        tCGGCGCTGTGGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCg         <  7:208459/148‑1 (MQ=255)
                                                                                                                                        tCGGCGCTGTGGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCg         <  7:202327/148‑1 (MQ=255)
                                                                                                                                         cGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGCCTCCAGCCGCAGGTGTTCCCTCACGCGGGAAACTGCGCTGAGATCTACCgg        >  3:333753/1‑148 (MQ=255)
                                                                                                                                         cGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCTAGCAGCAGGTTTTCCATc                                     <  7:279931/119‑1 (MQ=255)
                                                                                                                                         cGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCAAGCAGCAGGTTTTCCATc                                     >  8:279931/1‑119 (MQ=255)
                                                                                                                                            cgcTGTAGCGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTACCGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACtt                        >  7:360898/1‑129 (MQ=255)
                                                                                                                                            cgcTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACtt                        <  8:360898/129‑1 (MQ=255)
                                                                                                                                              cTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGTAAt   <  7:329955/148‑1 (MQ=255)
                                                                                                                                                gTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGTAATc  >  7:436946/1‑147 (MQ=255)
                                                                                                                                                gTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGTAATc  >  2:394447/1‑147 (MQ=255)
                                                                                                                                                 tAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAAc                          >  2:97476/1‑122 (MQ=255)
                                                                                                                                                 tAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAAc                          <  1:97476/122‑1 (MQ=255)
                                                                                                                                                 tAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGTAATc  >  2:87515/1‑146 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGCTGGACGCCCTGGAACAGCAGGCGGTTAGGTTCGCCGTCAATCCACAGCATCCCTTTGTAACGTAGCAGTTTATCCGCCGACTCCAGCAGCAGGTTTTCCATCACGCGGGAAACTTCGCTGATATCTACCGGGTAATC  >  NC_000913/4587882‑4588172

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