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,5430GA98.8% 270.5 / ‑3.7 84coding (648/1479 nt)ygcEputative kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (41/42);  minor base . (1/0);  total (42/42)
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.11e-01

CATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTG  >  NC_000913/2902409‑2902675
                                                                                                                                      |                                                                                                                                    
cATATGCTGTTTGATGTGCAAATGCCTGGCATCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCg                                                                                                                                  <  1:231703/139‑1 (MQ=255)
         tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTgg                                                                                                                         >  2:237350/1‑139 (MQ=255)
         tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTgg                                                                                                                         >  2:1001452/1‑139 (MQ=255)
         tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTgg                                                                                                                         >  1:1101252/1‑139 (MQ=255)
         tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTgg                                                                                                                         >  2:904192/1‑139 (MQ=255)
              tgtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGa                                                                                                                                <  1:828043/127‑1 (MQ=255)
              tgtgCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGa                                                                                                                                >  2:828043/1‑127 (MQ=255)
                  cAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGaaa                                                                                                                <  1:1001452/139‑1 (MQ=255)
                   aaaTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAAc                                                                                                               >  2:881533/1‑139 (MQ=255)
                       gCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGAGAGCCGGATTACTGGATGATGAAAcggcg                                                                                                           >  2:717309/1‑139 (MQ=255)
                            gCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAAt                                                                                                      >  2:304543/1‑139 (MQ=255)
                                  tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAAc                                                                                                               <  1:581305/124‑1 (MQ=255)
                                  tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAAc                                                                                                               >  2:581305/1‑124 (MQ=255)
                                  tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTctct                                                                                                >  1:1113289/1‑139 (MQ=255)
                                  tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTctct                                                                                                >  1:231114/1‑139 (MQ=255)
                                   ccTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                                                               <  2:478880/139‑1 (MQ=255)
                                    cTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg                                                                                              <  2:1113289/139‑1 (MQ=255)
                                      cGGTCATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc                                                                                            <  1:717309/139‑1 (MQ=255)
                                            tatCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCCCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACa                                                                                     >  1:461711/1‑139 (MQ=255)
                                             atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACa                                                                                     <  1:304543/139‑1 (MQ=255)
                                              tCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCAc                                                                                          >  2:158920/1‑133 (MQ=255)
                                              tCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCAc                                                                                          <  1:158920/133‑1 (MQ=255)
                                              tCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                                                    <  1:863001/139‑1 (MQ=255)
                                              tCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                                                    <  1:665802/139‑1 (MQ=255)
                                                 caccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCCCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCa                                                                                <  2:461711/139‑1 (MQ=255)
                                                      caAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCg                                                                                                                                  >  2:893527/1‑85 (MQ=255)
                                                      caAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCg                                                                                                                                  <  1:893527/85‑1 (MQ=255)
                                                       aaGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCg                                                                                  <  2:382821/132‑1 (MQ=255)
                                                       aaGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCg                                                                                  >  1:382821/1‑132 (MQ=255)
                                                       aaGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAt                                                                               <  2:265088/135‑1 (MQ=255)
                                                       aaGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAt                                                                               >  1:265088/1‑135 (MQ=255)
                                                         gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc                                                                                                 >  2:840401/1‑115 (MQ=255)
                                                         gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc                                                                                                 <  1:840401/115‑1 (MQ=255)
                                                         gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatga                                                                         >  1:75455/1‑139 (MQ=255)
                                                         gCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatga                                                                         >  1:659227/1‑139 (MQ=255)
                                                              aCTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCa                                                                                <  2:61684/127‑1 (MQ=255)
                                                              aCTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCa                                                                                >  1:61684/1‑127 (MQ=255)
                                                                            tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGga                                                                                                                        >  2:617820/1‑73 (MQ=255)
                                                                            tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGga                                                                                                                        <  1:617820/73‑1 (MQ=255)
                                                                            tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAAt                                                                                                      <  1:966166/91‑1 (MQ=255)
                                                                            tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAAt                                                                                                      >  2:966166/1‑91 (MQ=255)
                                                                               cGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                                                   <  2:75455/139‑1 (MQ=255)
                                                                                 gCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGaa                                                                     >  2:52956/1‑119 (MQ=255)
                                                                                 gCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGaa                                                                     <  1:52956/119‑1 (MQ=255)
                                                                                  cAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGG‑‑‑AGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATcc                                             >  1:651966/1‑139 (MQ=255)
                                                                                    ggACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc                                              <  1:904192/139‑1 (MQ=255)
                                                                                        tGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGt                                          >  2:257751/1‑139 (MQ=255)
                                                                                         gCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTa                                         <  2:1101252/139‑1 (MQ=255)
                                                                                                 tttgtACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTg                                 <  2:1013281/139‑1 (MQ=255)
                                                                                                    gtACCACCAGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACAGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGcc                              >  1:862026/1‑139 (MQ=255)
                                                                                                          ccaGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg                                                                                              <  2:996339/69‑1 (MQ=255)
                                                                                                          ccaGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg                                                                                              >  1:996339/1‑69 (MQ=255)
                                                                                                          ccaGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTAt                        >  2:966713/1‑139 (MQ=255)
                                                                                                           caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTa                                                                                       >  2:155914/1‑75 (MQ=255)
                                                                                                           caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTa                                                                                       <  1:155914/75‑1 (MQ=255)
                                                                                                           caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCa                                                           >  1:65137/1‑103 (MQ=255)
                                                                                                           caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCa                                                           <  2:65137/103‑1 (MQ=255)
                                                                                                           caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTg                                 >  2:248352/1‑129 (MQ=255)
                                                                                                           caGTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTg                                 <  1:248352/129‑1 (MQ=255)
                                                                                                             gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTAt                        <  1:827559/136‑1 (MQ=255)
                                                                                                             gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTAt                        >  2:827559/1‑136 (MQ=255)
                                                                                                             gTGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGtc                     >  2:20509/1‑139 (MQ=255)
                                                                                                              tGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                                                    <  1:759337/75‑1 (MQ=255)
                                                                                                              tGATAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGCCACTTACAt                                                                                    >  2:759337/1‑75 (MQ=255)
                                                                                                                aTAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCACAt                                                                               <  1:1091680/78‑1 (MQ=255)
                                                                                                                aTAAACCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCACAt                                                                               >  2:1091680/1‑78 (MQ=255)
                                                                                                                   aaCCGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTAt               <  1:237350/139‑1 (MQ=255)
                                                                                                                     ccGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc                                              <  1:98598/106‑1 (MQ=255)
                                                                                                                     ccGGTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc                                              >  2:98598/1‑106 (MQ=255)
                                                                                                                        gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGaaa                                                                                                                <  2:552005/37‑1 (MQ=38)
                                                                                                                        gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGaaa                                                                                                                >  1:552005/1‑37 (MQ=255)
                                                                                                                        gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAAcg                                                                                                              <  1:1003260/39‑1 (MQ=38)
                                                                                                                        gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAAcg                                                                                                              >  2:1003260/1‑39 (MQ=255)
                                                                                                                        gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                                                    <  1:669169/65‑1 (MQ=255)
                                                                                                                        gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                                                    >  2:669169/1‑65 (MQ=255)
                                                                                                                        gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCa                                                                                <  1:481150/69‑1 (MQ=255)
                                                                                                                        gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCa                                                                                >  2:481150/1‑69 (MQ=255)
                                                                                                                        gTAGAAGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGc          <  2:659227/139‑1 (MQ=255)
                                                                                                                             aGCTCTGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACa     <  2:231114/139‑1 (MQ=255)
                                                                                                                                tctGGGAGCGGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTg  <  1:711432/139‑1 (MQ=255)
                                                                                                                                tctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCAc                                                          >  2:866043/1‑83 (MQ=255)
                                                                                                                                tctGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCAc                                                          <  1:866043/83‑1 (MQ=255)
                                                                                                                                  tGGGAGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc                                              <  2:722836/93‑1 (MQ=255)
                                                                                                                                  tGGGAGCCAGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATc                                              >  1:722836/1‑93 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACTGGCCGATTATGTCTTCTATTCCGCAAACATTG  >  NC_000913/2902409‑2902675

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