Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,061,993 T→C 5.4% Y277H (TAT→CAT)  yihM → putative sugar phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,061,9930TC 5.4% 261.5 / 3.9 74Y277H (TAT→CAT) yihMputative sugar phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base C (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAAGAGGCAAACCATCTTC  >  NC_000913/4061851‑4062134
                                                                                                                                              |                                                                                                                                             
gCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGc                                                                                                                                          <  2:364871/148‑1 (MQ=255)
gCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGc                                                                                                                                          >  1:402086/1‑148 (MQ=255)
 ccGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGcc                                                                                                                                         >  1:318122/1‑148 (MQ=255)
   gTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAg                                                                                                                                       <  3:939286/148‑1 (MQ=255)
    ttCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGa                                                                                                                                      <  1:916491/148‑1 (MQ=255)
        aaGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAg                                                                                                                                  <  3:988272/148‑1 (MQ=255)
         aGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGt                                                                                                                                 <  2:318122/148‑1 (MQ=255)
            gTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGaa                                                                                                                              <  4:346967/148‑1 (MQ=255)
                   aCGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGcc                                                                                                                                         >  4:539802/1‑130 (MQ=255)
                   aCGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGcc                                                                                                                                         <  3:539802/130‑1 (MQ=255)
                        ccTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTCATctggcg                                                                                                                                       >  3:598976/1‑122 (MQ=255)
                        ccTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTCATctggcg                                                                                                                                       >  1:601355/1‑122 (MQ=255)
                        ccTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTCATctggcg                                                                                                                                       <  2:601355/127‑6 (MQ=255)
                        ccTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTCATctggcg                                                                                                                                       <  4:598976/127‑6 (MQ=255)
                            aTCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGaaa                                                                                                              <  2:402086/148‑1 (MQ=255)
                                              gaGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGa                                                                                                                               <  2:579516/113‑1 (MQ=255)
                                              gaGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGa                                                                                                                               >  1:579516/1‑113 (MQ=255)
                                                      ggTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCa                                                                                                                        >  3:954827/1‑112 (MQ=255)
                                                      ggTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCa                                                                                                                        <  4:954827/112‑1 (MQ=255)
                                                               atatGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGa                                                                                                                                   <  3:632979/92‑1 (MQ=255)
                                                               atatGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGa                                                                                                                                   >  4:632979/1‑92 (MQ=255)
                                                                   ggCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGt                                                                                         <  1:1058529/130‑1 (MQ=255)
                                                                   ggCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGt                                                                                         >  2:1058529/1‑130 (MQ=255)
                                                                      cTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGaa                                                                                                               >  2:381437/1‑105 (MQ=255)
                                                                      cTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGaa                                                                                                               <  1:381437/105‑1 (MQ=255)
                                                                        ggAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGc                                                                                                                                          >  1:476166/1‑76 (MQ=255)
                                                                        ggAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGc                                                                                                                                          <  2:476166/76‑1 (MQ=255)
                                                                          aagaagAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACAcc                                                                                                                         >  1:283520/1‑91 (MQ=255)
                                                                          aagaagAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACAcc                                                                                                                         <  2:283520/91‑1 (MQ=255)
                                                                          aagaagAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaaga                                                                               >  3:20862/1‑133 (MQ=255)
                                                                          aagaagAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaaga                                                                               <  4:20862/133‑1 (MQ=255)
                                                                             aagaGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAaca                                                                                                                           <  2:1017759/86‑1 (MQ=255)
                                                                             aagaGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAaca                                                                                                                           >  1:1017759/1‑86 (MQ=255)
                                                                             aagaGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTaa                                                             >  2:1021577/1‑148 (MQ=255)
                                                                              agagGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAAt                                                            <  1:1021577/148‑1 (MQ=255)
                                                                                    tGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTa                                                                                                      >  4:56459/1‑100 (MQ=255)
                                                                                    tGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTa                                                                                                      <  3:56459/100‑1 (MQ=255)
                                                                                        tACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTg                                                                                                                     >  4:1076863/1‑81 (MQ=255)
                                                                                        tACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTg                                                                                                                     <  3:1076863/81‑1 (MQ=255)
                                                                                            aTGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTAc                                                                                                     <  4:386674/93‑1 (MQ=255)
                                                                                            aTGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTAc                                                                                                     >  3:386674/1‑93 (MQ=255)
                                                                                                       ttCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTaa                                                             <  1:744303/122‑1 (MQ=255)
                                                                                                       ttCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTaa                                                             >  2:744303/1‑122 (MQ=255)
                                                                                                         ccGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGtacta                                 >  1:204211/1‑148 (MQ=255)
                                                                                                          cGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGtactac                                >  4:318926/1‑148 (MQ=255)
                                                                                                             tttGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGAt                                                                     <  1:1020011/108‑1 (MQ=255)
                                                                                                             tttGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGAt                                                                     >  2:1020011/1‑108 (MQ=255)
                                                                                                               tGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTAcaacaa                           <  2:204211/148‑1 (MQ=255)
                                                                                                                   gACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAAt                                                            <  3:882464/111‑1 (MQ=255)
                                                                                                                   gACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAAt                                                            >  4:882464/1‑111 (MQ=255)
                                                                                                                   gACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCa                       >  1:1082890/1‑148 (MQ=255)
                                                                                                                    aCGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCAtg                                               <  1:1055916/123‑1 (MQ=255)
                                                                                                                    aCGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCAtg                                               >  2:1055916/1‑123 (MQ=255)
                                                                                                                    aCGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTAcaacaa                           <  1:501793/143‑1 (MQ=255)
                                                                                                                    aCGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTAcaacaa                           >  2:501793/1‑143 (MQ=255)
                                                                                                                    aCGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCaa                      <  3:318926/148‑1 (MQ=255)
                                                                                                                         gatgatAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAag                 <  2:1082890/148‑1 (MQ=255)
                                                                                                                           tgatAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAagag               >  2:338845/1‑148 (MQ=255)
                                                                                                                             atAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATc                                                           <  3:1140631/102‑1 (MQ=255)
                                                                                                                             atAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATc                                                           >  4:1140631/1‑102 (MQ=255)
                                                                                                                               aaTCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCa                                                 >  1:892522/1‑110 (MQ=255)
                                                                                                                               aaTCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCa                                                 <  2:892522/110‑1 (MQ=255)
                                                                                                                                aTCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAAGAGGCaaa          >  1:767449/1‑148 (MQ=255)
                                                                                                                                 tCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAagag               <  2:322839/142‑1 (MQ=255)
                                                                                                                                 tCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAAGTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAagag               >  1:322839/1‑142 (MQ=255)
                                                                                                                                  ccGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaagaaaaagaa                                                                        >  4:562758/1‑84 (MQ=255)
                                                                                                                                  ccGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaagaaaaagaa                                                                        <  3:562758/84‑1 (MQ=255)
                                                                                                                                  ccGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCACTTAATCAGATAAATCAtg                                               <  2:863786/109‑1 (MQ=255)
                                                                                                                                  ccGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCAtg                                               >  1:863786/1‑109 (MQ=255)
                                                                                                                                     tGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaagaaaa                                                                            >  2:569300/1‑77 (MQ=255)
                                                                                                                                     tGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTaaagaaaa                                                                            <  1:569300/77‑1 (MQ=255)
                                                                                                                                        aTCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAAGAGGCAAACCAtcttc  <  1:338845/148‑1 (MQ=255)
                                                                                                                                         tCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGta                                    >  2:277984/1‑113 (MQ=255)
                                                                                                                                         tCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGta                                    <  1:277984/113‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GCCGTTCAAAGCGTTATTAACGCACCTTATCTGCCTGGGTGATGATGAGCCGCAGGTGACGGCATATGGCCTGGAAGAAGAGGTTGATTACTATGCGCCGGCGTTCCGCTTTGAAGACGAAGATGATAATCCGTGGATCCCTTATCGCCAGATGAGTGAAACACCATTGCCAGAAAATCATTTACTGGATGCGCGGTTACGTAAAGAAAAAGAAGATGCAATTAATCAGATAAATCATGTGCGTAACGTACTACAACAAATCAAACAAGAGGCAAACCATCTTC  >  NC_000913/4061851‑4062134

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