Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,373,288 C→T 5.2% intergenic (+126/+77) yjeI → / ← yjeJ DUF4156 family lipoprotein/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,373,2880CT 5.2% 458.2 / 8.7 137intergenic (+126/+77)yjeI/yjeJDUF4156 family lipoprotein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (68/61);  new base T (4/3);  total (73/64)
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.43e-01

Only 100 of 145 total aligned reads displayed.

AGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGA  >  NC_000913/4373140‑4373434
                                                                                                                                                    |                                                                                                                                                  
aGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                                    >  4:808747/1‑149 (MQ=255)
aGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                                    >  1:838244/1‑149 (MQ=255)
   tttATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAACTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCGTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                 >  2:326593/1‑149 (MQ=255)
      ataAGTGCCCGAACTGATGTGAGCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                              <  4:761570/149‑1 (MQ=255)
        aaGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                            <  1:395288/149‑2 (MQ=38)
        aaGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                            <  4:367265/149‑1 (MQ=255)
         aGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                            >  1:275973/1‑147 (MQ=255)
          gTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATg                                                                                                                                          >  4:73102/1‑149 (MQ=255)
            gCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAc                                                                                                                                        >  1:655275/1‑149 (MQ=255)
              ccGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCa                                                                                                                                      >  2:779030/1‑149 (MQ=255)
              ccGAACTGATGTGATCCAGAGGCTGGGTGAAAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCa                                                                                                                                      >  3:773731/1‑149 (MQ=255)
                   cTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGAt                                                                                                                                 <  1:851604/149‑1 (MQ=255)
                    tGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATc                                                                                                                                >  3:140632/1‑149 (MQ=255)
                      aTGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAt                                                                                                                                           >  4:623097/1‑136 (MQ=255)
                      aTGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCAc                                                                                                                              >  2:880359/1‑149 (MQ=255)
                        gtgATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGt                                                                                                                            >  2:211516/1‑149 (MQ=255)
                         tgATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGtt                                                                                                                           >  2:768877/1‑149 (MQ=255)
                           aTCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTAt                                                                                                                         >  1:918784/1‑149 (MQ=255)
                           aTCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGAGCTTCGCGTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTAAt                                                                                                                         >  3:384788/1‑149 (MQ=255)
                            tCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTATGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATc                                                                                                                        >  4:89166/1‑149 (MQ=255)
                             ccAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCa                                                                                                                       >  1:514946/1‑149 (MQ=255)
                             ccAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCa                                                                                                                       <  1:477898/149‑1 (MQ=255)
                             ccAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCa                                                                                                                       <  4:915197/149‑1 (MQ=255)
                              cAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAg                                                                                                                      >  2:670604/1‑149 (MQ=255)
                                     tGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTc                                                                                                               >  3:143843/1‑149 (MQ=255)
                                       ggTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAg                                                                                                             <  4:143843/149‑1 (MQ=255)
                                        gTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAgg                                                                                                            <  4:239294/149‑1 (MQ=255)
                                          gACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCt                                                                                                          >  2:219434/1‑149 (MQ=255)
                                             aaaCGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCta                                                                                                         <  4:694666/147‑2 (MQ=255)
                                             aaaCGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCta                                                                                                         >  3:694666/1‑146 (MQ=255)
                                             aaaCGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGa                                                                                                       >  1:888309/1‑149 (MQ=255)
                                             aaaCGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGa                                                                                                       >  2:632876/1‑149 (MQ=255)
                                             aaaCGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGa                                                                                                       >  2:504210/1‑149 (MQ=255)
                                                 gTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGa                                                                                                       <  1:594024/145‑1 (MQ=255)
                                                 gTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGa                                                                                                       >  2:594024/1‑145 (MQ=255)
                                                    taATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTcc                                                                                                <  1:211516/149‑1 (MQ=255)
                                                    taATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTcc                                                                                                >  2:401540/1‑149 (MQ=255)
                                                          ccTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTc                                                                                          >  2:12506/1‑149 (MQ=255)
                                                           cTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCg                                                                                         >  3:248618/1‑149 (MQ=255)
                                                            tGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTTTCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGc                                                                                        >  1:108551/1‑149 (MQ=255)
                                                            tGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGc                                                                                        >  2:395383/1‑149 (MQ=255)
                                                            tGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGAGCCATCCAGCCtgaa                                                                                        <  1:662812/149‑4 (MQ=38)
                                                                     tcgcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCaaa                                                                               <  4:773731/149‑1 (MQ=255)
                                                                        cttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGAt                                                                            <  2:498747/149‑1 (MQ=255)
                                                                        cttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCAGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCCAAGAt                                                                            >  2:837536/1‑149 (MQ=255)
                                                                          tATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATga                                                                          <  3:691384/149‑1 (MQ=255)
                                                                          tATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATga                                                                          <  2:918784/149‑1 (MQ=255)
                                                                            tCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGTATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATgaga                                                                        <  4:233213/149‑1 (MQ=255)
                                                                            tCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGTATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATgaga                                                                        >  3:233213/1‑149 (MQ=255)
                                                                            tCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATgaga                                                                        <  1:204135/149‑1 (MQ=255)
                                                                            tCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATgaga                                                                        <  3:587408/149‑1 (MQ=255)
                                                                             cAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGat                                                                       <  3:73102/149‑1 (MQ=255)
                                                                                gCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGattat                                                                    >  3:856076/1‑149 (MQ=255)
                                                                                   tACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTc                                                                 >  3:710074/1‑149 (MQ=255)
                                                                                     cATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGc                                                               <  1:401540/149‑1 (MQ=255)
                                                                                      aTAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATACGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGcc                                                              >  2:666004/1‑149 (MQ=255)
                                                                                         aGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgt                                                           <  1:219434/149‑1 (MQ=255)
                                                                                         aGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgt                                                           <  4:856076/149‑1 (MQ=255)
                                                                                           gCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTc                                                                                          <  2:233531/116‑1 (MQ=255)
                                                                                               tGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTTTCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCa                                                     >  1:359254/1‑149 (MQ=255)
                                                                                                 cAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGt                                                   <  1:779030/149‑1 (MQ=255)
                                                                                                  aaCTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGtt                                                  >  1:555239/1‑149 (MQ=255)
                                                                                                  aaCTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGtt                                                  >  4:282610/1‑149 (MQ=255)
                                                                                                    cTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGtttt                                                <  2:316578/149‑1 (MQ=255)
                                                                                                      gaatGAATGCTAGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATAAGTCTTTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCCAAGACGAGATTATTTCGCCTGTGGTGCAGTTTTGc                                              <  4:571832/146‑1 (MQ=38)
                                                                                                      gATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGc                                              <  3:668173/149‑1 (MQ=255)
                                                                                                       aTTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCt                                             >  1:881252/1‑149 (MQ=255)
                                                                                                       aTTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCt                                             >  3:202411/1‑149 (MQ=255)
                                                                                                           aaTGCTTGAGAGTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtggtg                                         >  3:525760/1‑149 (MQ=255)
                                                                                                           aaTGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtggtg                                         >  3:534281/1‑149 (MQ=255)
                                                                                                           aaTGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtggtg                                         <  3:307496/149‑1 (MQ=255)
                                                                                                           aaTGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtggtg                                         >  3:90323/1‑149 (MQ=255)
                                                                                                             tGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTAATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGa                                       <  2:655275/149‑1 (MQ=255)
                                                                                                              gCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgtg                                                          <  3:741801/129‑1 (MQ=255)
                                                                                                              gCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgtg                                                          >  4:741801/1‑129 (MQ=255)
                                                                                                              gCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGAt                                      <  2:355620/149‑1 (MQ=255)
                                                                                                               cTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGAtt                                     <  2:205271/149‑1 (MQ=255)
                                                                                                               cTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGAtt                                     <  4:258746/149‑1 (MQ=255)
                                                                                                                ttGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTg                                    >  4:196475/1‑149 (MQ=255)
                                                                                                                ttGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTg                                    >  4:440070/1‑149 (MQ=255)
                                                                                                                ttGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTg                                    >  4:603388/1‑149 (MQ=255)
                                                                                                                ttGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTg                                    >  1:667510/1‑149 (MQ=255)
                                                                                                                  gagaATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGaa                                  <  3:178300/149‑1 (MQ=255)
                                                                                                                         tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTaa                           <  3:603388/149‑1 (MQ=255)
                                                                                                                         tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTaa                           <  4:212698/149‑1 (MQ=255)
                                                                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtggt                                          >  3:581736/1‑129 (MQ=255)
                                                                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCATAACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCAGTGGTGCAGTTTTGCtggt                                          <  4:581736/129‑1 (MQ=255)
                                                                                                                               ccGGATGCGGCGTAAACGCCGTATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATg                     >  2:581335/1‑149 (MQ=255)
                                                                                                                                 ggATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGAtt                                     <  3:802511/131‑1 (MQ=255)
                                                                                                                                   aTAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAgc                 >  1:99123/1‑149 (MQ=255)
                                                                                                                                       ggCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAgcgcgc             >  1:307735/1‑149 (MQ=255)
                                                                                                                                          gTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGAtt                                     >  3:205828/1‑122 (MQ=255)
                                                                                                                                          gTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTa          >  4:153562/1‑149 (MQ=255)
                                                                                                                                          gTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTa          >  3:891210/1‑149 (MQ=255)
                                                                                                                                           tttACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTaa         <  2:888309/149‑1 (MQ=255)
                                                                                                                                              aCGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCa      >  4:125814/1‑149 (MQ=255)
                                                                                                                                              aCGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCa      <  2:677786/149‑1 (MQ=255)
                                                                                                                                               cgccgcATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCAt     <  3:808747/149‑1 (MQ=255)
                                                                                                                                                gccgcATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCAtt    >  1:146057/1‑149 (MQ=255)
                                                                                                                                                  cgcATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGa  <  4:202411/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGA  >  NC_000913/4373140‑4373434

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