Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,902,543 G→A G216G (GGG→GGA ygcE → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,902,5430GA99.0% 330.7 / NA 102G216G (GGG→GGAygcEputative kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (51/50);  minor base C (0/1);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

GTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTGC  >  NC_000913/2902407‑2902676
                                                                                                                                        |                                                                                                                                     
gTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGc                                                                                                                                     >  1:1072584/1‑139 (MQ=255)
  cATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCg                                                                                                                                   >  1:897659/1‑139 (MQ=255)
  cATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCg                                                                                                                                   <  1:280251/139‑1 (MQ=255)
                  tgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGAtgatga                                                                                                                   >  1:1802126/1‑139 (MQ=255)
                    cAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGaaa                                                                                                                 <  1:244107/139‑1 (MQ=255)
                        tGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAAcggc                                                                                                             <  2:897659/139‑1 (MQ=255)
                                   gTCCTCGGACATATCACACCACAAGCCGGACTGGCTACACATTTCCCGGCAGGACTGACATTCGTTAGTACCACTAGTGACAAACCGGTAGAAGCTCTGGGCGCCGGATTACTGGCTGATGAGACGGCGGTAATTtctc                                                                                                  <  1:374860/139‑1 (MQ=255)
                                    tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCg                                                                                                                                   <  2:982265/105‑1 (MQ=255)
                                    tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCg                                                                                                                                   >  1:982265/1‑105 (MQ=255)
                                       tCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGGTGTTTGTACCACCAGTGATAAACCGGTAGTAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTgg                                                                                              <  1:189722/139‑1 (MQ=255)
                                        cGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc                                                                                             <  1:1687294/139‑1 (MQ=255)
                                              tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                                >  2:1658559/1‑130 (MQ=255)
                                              tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                                <  1:1658559/130‑1 (MQ=255)
                                               atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGTCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACa                                                                                      >  2:1450084/1‑139 (MQ=255)
                                               atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc                                                                                                  <  1:1894612/127‑1 (MQ=255)
                                               atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc                                                                                                  >  2:1894612/1‑127 (MQ=255)
                                               atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                                <  2:1453856/129‑1 (MQ=255)
                                               atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACa                                                                                      >  2:1196634/1‑139 (MQ=255)
                                               atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACa                                                                                      >  2:2036777/1‑139 (MQ=255)
                                               atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTACCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                                >  1:1453856/1‑129 (MQ=255)
                                                     ccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAt                                                                                <  2:1184125/139‑1 (MQ=255)
                                                        caAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACtt                                                                                         <  1:1625149/127‑1 (MQ=255)
                                                        caAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACtt                                                                                         >  2:1625149/1‑127 (MQ=255)
                                                        caAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtga                                                                             >  2:1483717/1‑139 (MQ=255)
                                                        caAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtga                                                                             >  1:1129707/1‑139 (MQ=255)
                                                         aaGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAAcggc                                                                                                             >  2:954968/1‑106 (MQ=255)
                                                         aaGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAAcggc                                                                                                             <  1:954968/106‑1 (MQ=255)
                                                         aaGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                                                            >  2:569227/1‑139 (MQ=255)
                                                         aaGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                                                            >  2:1396380/1‑139 (MQ=255)
                                                           gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                                >  1:1749186/1‑117 (MQ=255)
                                                           gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                                <  2:1749186/117‑1 (MQ=255)
                                                           gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatga                                                                          >  1:1749606/1‑139 (MQ=255)
                                                           gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatga                                                                          >  2:399628/1‑139 (MQ=255)
                                                           gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatga                                                                          >  1:587115/1‑139 (MQ=255)
                                                             cGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgatg                                                                        <  1:569227/139‑1 (MQ=255)
                                                                         acaTTTCCCGGCAGGACTGCGCGTTGGTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGCTTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCa                                                            <  1:1196634/139‑1 (MQ=255)
                                                                         acaTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCa                                                            <  1:399628/139‑1 (MQ=255)
                                                                             ttCCCGGCAGGACTGTCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCagtgctttgc                                                        <  1:1450084/139‑10 (MQ=255)
                                                                               cccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCg                                                      <  1:1396380/139‑1 (MQ=255)
                                                                                 cGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                                                    <  1:2036777/139‑1 (MQ=255)
                                                                                     aGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTgg                                                                                              <  2:896847/93‑1 (MQ=255)
                                                                                     aGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTgg                                                                                              >  1:896847/1‑93 (MQ=255)
                                                                                     aGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGAt                                                >  1:1458858/1‑139 (MQ=255)
                                                                                         cTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAg                                            >  1:500632/1‑139 (MQ=255)
                                                                                           gCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTa                                          <  2:500632/139‑1 (MQ=255)
                                                                                                    ttgtACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTgg                                 <  2:879223/139‑1 (MQ=255)
                                                                                                     tgtACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGc                                <  2:1129707/139‑1 (MQ=255)
                                                                                                      gtACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                                <  2:929113/74‑1 (MQ=255)
                                                                                                      gtACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                                >  1:929113/1‑74 (MQ=255)
                                                                                                             caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAAtt                                                                                                      >  1:1568583/1‑61 (MQ=255)
                                                                                                             caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAAtt                                                                                                      <  2:1568583/61‑1 (MQ=255)
                                                                                                             caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATg                        >  1:280104/1‑139 (MQ=255)
                                                                                                              aGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAAt                                                                                                       >  2:1871219/1‑59 (MQ=255)
                                                                                                              aGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAAt                                                                                                       <  1:1871219/59‑1 (MQ=255)
                                                                                                               gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTa                                          <  2:877552/119‑1 (MQ=255)
                                                                                                               gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTa                                          >  1:877552/1‑119 (MQ=255)
                                                                                                               gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGtc                      <  2:1875663/139‑1 (MQ=255)
                                                                                                                tGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCa                                                                 >  2:132498/1‑95 (MQ=255)
                                                                                                                tGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCa                                                                 <  1:132498/95‑1 (MQ=255)
                                                                                                                 gATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTgg                                                                                              >  1:1364689/1‑65 (MQ=255)
                                                                                                                 gATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTgg                                                                                              <  2:1364689/65‑1 (MQ=255)
                                                                                                                 gATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGtctt                    >  2:799533/1‑139 (MQ=255)
                                                                                                                    aaaCCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTa                 <  2:1458858/139‑1 (MQ=255)
                                                                                                                     aaCCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTAt                <  2:587115/139‑1 (MQ=255)
                                                                                                                        cGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCa                                                            >  1:1616383/1‑92 (MQ=255)
                                                                                                                        cGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCa                                                            <  2:1616383/92‑1 (MQ=255)
                                                                                                                        cGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGa                                                 >  2:1556463/1‑103 (MQ=255)
                                                                                                                        cGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGa                                                 <  1:1556463/103‑1 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGTATTACTGGATGATGaa                                                                                                                  <  2:426175/36‑1 (MQ=37)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGaa                                                                                                                  >  1:426175/1‑36 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatg                                                                           >  1:627527/1‑75 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatg                                                                           <  2:627527/75‑1 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                                                    <  1:23167/98‑1 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                                                    >  2:23167/1‑98 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTg                                  <  1:1186651/116‑1 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTg                                  >  2:1186651/1‑116 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGAtt                           <  1:1557042/123‑1 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGAtt                           >  2:1557042/1‑123 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGc           >  2:507902/1‑139 (MQ=255)
                                                                                                                          gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGc           >  2:1606129/1‑139 (MQ=255)
                                                                                                                           tAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc                                                                                             <  1:2000111/56‑1 (MQ=255)
                                                                                                                           tAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc                                                                                             >  2:2000111/1‑56 (MQ=255)
                                                                                                                            aGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACt                                                                                          <  1:1924978/58‑1 (MQ=255)
                                                                                                                            aGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACt                                                                                          >  2:1924978/1‑58 (MQ=255)
                                                                                                                            aGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCaa         >  2:1161479/1‑139 (MQ=255)
                                                                                                                            aGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCaa         >  1:1100986/1‑139 (MQ=255)
                                                                                                                               aGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACa      <  1:799533/139‑1 (MQ=255)
                                                                                                                                  tctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc                                                                                             >  1:1334925/1‑49 (MQ=255)
                                                                                                                                  tctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTg   <  1:1606129/139‑1 (MQ=255)
                                                                                                                                  tctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTg   <  2:1749606/139‑1 (MQ=255)
                                                                                                                                  tctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTGCTCTTGGc                                                                                             <  2:1334925/49‑1 (MQ=38)
                                                                                                                                   ctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTGc  >  2:633198/1‑139 (MQ=255)
                                                                                                                                    tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc                                                                                                  >  1:1367094/1‑42 (MQ=255)
                                                                                                                                    tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc                                                                                                  <  2:1367094/42‑1 (MQ=255)
                                                                                                                                    tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                                                    >  2:237625/1‑88 (MQ=255)
                                                                                                                                    tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                                                    <  1:237625/88‑1 (MQ=255)
                                                                                                                                    tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTa                                     <  2:354442/103‑1 (MQ=255)
                                                                                                                                    tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTa                                     >  1:354442/1‑103 (MQ=255)
                                                                                                                                     gggagCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGa                                                 >  2:1502719/1‑90 (MQ=255)
                                                                                                                                     gggagCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGtctt                    <  2:1503088/119‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
GTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTGC  >  NC_000913/2902407‑2902676

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