Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 805,312 G→T 6.7% E604* (GAG→TAG)  ybhJ → aconitase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913805,3120GT 6.7% 267.1 / 2.5 93E604* (GAG→TAG) ybhJaconitase family protein
Reads supporting (aligned to +/- strand):  ref base G (35/49);  new base T (2/6);  total (38/55)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGA  >  NC_000913/805163‑805459
                                                                                                                                                     |                                                                                                                                                   
aaCTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCg                                                                                                                                                     <  1:1467204/150‑1 (MQ=255)
aaCTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCg                                                                                                                                                     >  2:1798566/1‑150 (MQ=255)
   tGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGc                                                                                                                                                  >  1:1042356/1‑150 (MQ=255)
     cttccttccGCTGACACCTCTTCTTATCGTTCAAGTCCGATTGGGCGGGCGGAGGTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTg                                                                                                                                                <  1:352018/149‑1 (MQ=255)
     aTTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTg                                                                                                                                                >  1:1228991/1‑150 (MQ=255)
     aTTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTg                                                                                                                                                <  2:245103/150‑1 (MQ=255)
     aTTCCATCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTAATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTg                                                                                                                                                <  1:357355/150‑1 (MQ=255)
      ttccttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGa                                                                                                                                               <  2:1115802/150‑1 (MQ=255)
       tccttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGAc                                                                                                                                              >  1:1831919/1‑150 (MQ=255)
       tcattccGGTGAAACCTCTTCTTAACGTTCAAATCCGATTGGTCTGGCGGAGATAACCATGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGCCTGCTGAGCTGGAAAAACAGCGACTGGCGGGGAATGTCAGCGAGCTCAc                                                                                                                                              >  1:1856576/4‑150 (MQ=255)
         cttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACag                                                                                                                                            >  1:1391433/1‑150 (MQ=255)
          ttccGGTGAAACCTCTTCTTCTCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACatt                                                                                                                                           <  1:680674/150‑3 (MQ=255)
           tccGGTTAAACCTCTTCTTATCTTTCAAAACAGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCTATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACagag                                                                                                                                          <  2:1668787/150‑1 (MQ=255)
           tccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACagag                                                                                                                                          >  2:497107/1‑150 (MQ=255)
            ccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAgg                                                                                                                                         <  1:199411/150‑1 (MQ=255)
            ccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAgg                                                                                                                                         <  1:855111/150‑1 (MQ=255)
             cGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattccc                                                                                                                                        <  2:1642601/150‑6 (MQ=255)
              ggTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGACGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgtg                                                                                                                                       <  1:1421187/150‑1 (MQ=255)
              ggTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgtg                                                                                                                                       >  2:1803039/1‑150 (MQ=255)
                        cttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCa                                                                                                                             <  1:1105685/150‑1 (MQ=255)
                         ttcttATCGTTCAAATCCGATTGGGCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGAAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAt                                                                                                                            <  2:1043459/150‑1 (MQ=255)
                          tcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                           <  2:557185/150‑1 (MQ=255)
                           cttTTCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGATGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTa                                                                                                                          <  1:1432359/150‑1 (MQ=255)
                             tATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGGGATCCCGGTTTTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAg                                                                                                                        <  1:773677/150‑1 (MQ=255)
                             tATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAg                                                                                                                        >  2:447252/1‑150 (MQ=255)
                               tCGTTCAAATCCGCTTGGTCAGGCGGAGTTTACCCGGTCTCGCCGCGACCCCGGTTATGTTAGCACACGAAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCa                                                                                                                      <  2:1382385/150‑1 (MQ=255)
                                 gTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCCGCGAGCTGACAGAGGTGTTTGAGCGCATTAAGCAGa                                                                                                                    >  2:647815/1‑150 (MQ=255)
                                 gTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGa                                                                                                                    >  2:211991/1‑150 (MQ=255)
                                   tCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGAtt                                                                                                                  <  2:1769510/150‑1 (MQ=255)
                                     aaaTCCGATTGTGCTGGCGTAGTTTACCCTGTCTCGCCGCGATCGCGGGTATGTTAGCAGAAGTAAATCGACTGCTGAGCTGTAAAAGCAGCGTCTGGCGGGGAATGTCAGCTAGCTGACAGAGTTGTTTGCGCGCATTAAGCAGATTGc                                                                                                                <  1:771731/150‑1 (MQ=255)
                                     aaaTCCGATTGGTCTGGCGGAGGTTACCCTGGCTCGCCGCGATCCCGGTAATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGGTTGCGCGCATTAAGCAGATTGc                                                                                                                <  2:391266/150‑1 (MQ=255)
                                         ccGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGt                                                                                                            >  2:1310475/1‑150 (MQ=255)
                                           gaagggccTGTCGTAGTTAACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCCGAACTAAAGCGACTGCTGACCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCa                                                                                                          <  2:1170226/143‑1 (MQ=255)
                                           gATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCa                                                                                                          >  1:658279/1‑150 (MQ=255)
                                            aTTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATTTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGTCGGGGAATGTCAGCCAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAg                                                                                                         >  1:111953/1‑150 (MQ=255)
                                            aTTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAg                                                                                                         <  1:1723406/150‑1 (MQ=255)
                                               ggTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                      <  2:1765491/150‑1 (MQ=255)
                                                 tCTTGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGACGTAAAGCGACTGCTGAGCGGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa                                                                                                    <  2:1399614/150‑1 (MQ=255)
                                                 tCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGCCAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa                                                                                                    >  2:1446598/1‑150 (MQ=255)
                                                       ggAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGa                                                                                              >  1:1145704/1‑150 (MQ=255)
                                                       ggAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGa                                                                                              <  2:1019224/150‑1 (MQ=255)
                                                        gAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCTTTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                             <  1:497107/150‑1 (MQ=255)
                                                        gAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                             <  2:1088642/150‑1 (MQ=255)
                                                         aGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATc                                                                                            <  2:1797834/150‑1 (MQ=255)
                                                         aGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATc                                                                                            <  2:339035/150‑1 (MQ=255)
                                                            ttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGGGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGc                                                                                         <  2:1748507/150‑1 (MQ=255)
                                                               cccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCCGGAGCATATTGATCCGCTGc                                                                                      >  2:1025855/1‑150 (MQ=255)
                                                               cccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGc                                                                                      >  2:241106/1‑150 (MQ=255)
                                                                   gTCTCGCCGCGATCCCGTTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                                                                  <  2:727159/150‑1 (MQ=255)
                                                                   gTCTCGCCGCGATCACGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGATTGCGCGCATTAAGCAGAATGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                                                                  <  1:1710278/150‑1 (MQ=255)
                                                                    tctcCCCGCGATCCCGGTTATGTTAGCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACt                                                                                 <  2:415373/150‑1 (MQ=255)
                                                                             gATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGCCGGGGAATGTCAGCTAGCTGACAGAGGTGTTTCCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgc                                                                        >  2:496002/1‑150 (MQ=255)
                                                                              aTCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgca                                                                       >  2:1502567/1‑150 (MQ=255)
                                                                               tCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcag                                                                      <  1:1310475/150‑1 (MQ=255)
                                                                               tCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcag                                                                      >  2:1034111/1‑150 (MQ=255)
                                                                               tCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcag                                                                      >  1:1130046/1‑150 (MQ=255)
                                                                                  cGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                                                                   <  1:1798566/150‑1 (MQ=255)
                                                                                     ttATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGCGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGt                                                                <  1:22807/150‑1 (MQ=255)
                                                                                        tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa                                                             >  1:666392/1‑150 (MQ=255)
                                                                                         gTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGGAGCATGGTCTAt                                                            >  1:341215/1‑150 (MQ=255)
                                                                                         gTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTAt                                                            >  2:750897/1‑150 (MQ=255)
                                                                                         gTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTAt                                                            >  1:1455512/1‑150 (MQ=255)
                                                                                          ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTCGCAGCATGGTCTAtt                                                           >  2:123447/1‑149 (MQ=255)
                                                                                           tAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGAGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGc                                                          >  1:698212/1‑150 (MQ=255)
                                                                                              cAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGt                                                       <  2:1177129/150‑1 (MQ=255)
                                                                                                gAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGa                                                     <  2:1674592/150‑1 (MQ=255)
                                                                                                 aaGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaa                                                    <  2:1070810/150‑1 (MQ=255)
                                                                                                 aaGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaa                                                    <  2:529410/150‑1 (MQ=255)
                                                                                                   ctAACGCGACTGCGGAGCTGGAAATTCAGCGTCTGGAGGTGAATGTCAGCGACCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGGTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAc                                                  <  2:81543/149‑1 (MQ=255)
                                                                                                     aaaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCa                                                >  1:317063/1‑150 (MQ=255)
                                                                                                       agctaatgctgAGCTGGAAAATGAGCGTCTGTCGGGGAATTTCAGCTAGCTGACAGTTTTGTTTGCGCTAATTAAGCAGTTTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTTTCAGCATGGTCTATGCGGTGAAACCAgg                                              <  1:1103622/144‑1 (MQ=255)
                                                                                                              gctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGt                                       <  1:1755862/150‑1 (MQ=255)
                                                                                                                    cTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTgcg                                 >  2:1676757/1‑150 (MQ=255)
                                                                                                                      ggAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGcgcg                               >  2:718354/1‑150 (MQ=255)
                                                                                                                      ggAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCACTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGcgcg                               <  2:1304914/150‑1 (MQ=255)
                                                                                                                        aaaaTCAGCGTCTGGCGGGGAATGGCAGCTAGCTGACAGAGGTGTTTTCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATTGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTg                             <  1:890789/150‑1 (MQ=255)
                                                                                                                           aTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAAc                          <  2:1042356/150‑1 (MQ=255)
                                                                                                                             cAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAg                        >  2:1733103/1‑150 (MQ=255)
                                                                                                                              aGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAgg                       >  2:504861/1‑150 (MQ=255)
                                                                                                                                cGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAATCTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAggcg                     >  2:439436/1‑150 (MQ=255)
                                                                                                                                cGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAggcg                     <  2:715157/150‑1 (MQ=255)
                                                                                                                                  tCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAggcggc                   >  2:615873/1‑150 (MQ=255)
                                                                                                                                  gcTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAggcggc                   <  1:1608203/149‑1 (MQ=255)
                                                                                                                                   cTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGgcggcg                  <  2:1145704/150‑1 (MQ=255)
                                                                                                                                   cTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGgcggcg                  >  1:1578224/1‑150 (MQ=255)
                                                                                                                                     ggCGGGGAATGTCAGCGCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAg                >  1:1490049/1‑150 (MQ=255)
                                                                                                                                     ggCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAg                <  1:447252/150‑1 (MQ=255)
                                                                                                                                         gggAATGTCCGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGGTCTGCGCGTGAACAGGCGGCGAGCTGc            <  1:1564331/150‑1 (MQ=255)
                                                                                                                                         gggAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGc            <  1:1553657/150‑1 (MQ=255)
                                                                                                                                           gAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCa          >  2:1365811/1‑150 (MQ=255)
                                                                                                                                            tctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAg         >  2:794263/3‑150 (MQ=255)
                                                                                                                                            aaTGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAg         <  2:1784641/150‑1 (MQ=255)
                                                                                                                                                tCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtg     <  2:1114497/150‑1 (MQ=255)
                                                                                                                                                  cgCGAGCTGACAGTGGTGTTTGCTCGCATTAAGCAGATTGCTGGTCATGAGCATATTGATCCGCTGCATACTGAAATTGGCAGGATGGTCTTTGGGGTGAAACGAGGCGATGTTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGtgtg   <  2:1856576/149‑1 (MQ=255)
                                                                                                                                                   gCGAGCTGACAGAGGTGTTTGCGCGCAGTCAGCAGATTGCTGGTCAGGAGCATACTGATCCGCTGCAAACTGAAATTGGCAGCAGGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGa  <  1:47744/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGA  >  NC_000913/805163‑805459

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