Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 231,007 C→A 7.4% intergenic (‑126/‑115) yafC ← / → yafD LysR family putative transcriptional regulator/endo/exonuclease/phosphatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,0070CA 7.4% 334.5 / 9.8 111intergenic (‑126/‑115)yafC/yafDLysR family putative transcriptional regulator/endo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (50/52);  new base A (3/6);  total (53/58)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

TTCTTCCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTATGTTGC  >  NC_000913/230861‑231156
                                                                                                                                                  |                                                                                                                                                     
ttcttcCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTg                                                                                                                                                    <  1:1725050/150‑1 (MQ=255)
  cttcCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCt                                                                                                                                                  >  2:1025555/1‑150 (MQ=255)
  cttcCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCt                                                                                                                                                  >  1:1706834/1‑150 (MQ=255)
    tcCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTa                                                                                                                                                >  1:24532/1‑150 (MQ=255)
    tcCGACGTGGCTTTCATTTTTGCTCATGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTa                                                                                                                                                <  2:1722761/150‑1 (MQ=255)
    tcCGACCTGCCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTACAAAATGAGCAAAATTCTGACCTGTGTAAGCACTTGCTTa                                                                                                                                                >  1:210486/1‑150 (MQ=255)
     ccGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTAc                                                                                                                                               >  2:1541901/1‑150 (MQ=255)
     ccGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTAc                                                                                                                                               <  2:386385/150‑1 (MQ=255)
      cGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGTATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACa                                                                                                                                              <  2:1647888/150‑1 (MQ=255)
        aCGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGCTGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACacc                                                                                                                                            <  2:89738/150‑3 (MQ=255)
        aCGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACacc                                                                                                                                            <  2:753319/150‑3 (MQ=255)
          gTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAATAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACaccgg                                                                                                                                          <  1:1347165/150‑5 (MQ=255)
          gTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAAc                                                                                                                                          <  1:1481245/150‑1 (MQ=255)
          gTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACaccgg                                                                                                                                          <  1:1215345/150‑5 (MQ=255)
          gTGGCTTTCATTTTTGCTCATGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAAc                                                                                                                                          <  2:135075/150‑1 (MQ=255)
            ggCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACaa                                                                                                                                        <  1:1170822/150‑1 (MQ=255)
             gCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACaccggtca                                                                                                                                       <  2:150475/150‑8 (MQ=255)
             gCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGAACTTGCTTACATAACAat                                                                                                                                       >  2:613271/1‑150 (MQ=255)
              cTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAata                                                                                                                                      >  2:1680500/1‑150 (MQ=255)
              cTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAATTGATTACATAACAata                                                                                                                                      >  2:543290/1‑150 (MQ=255)
              cTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACAGGTGTAAGCACTTGCTTACATAACAata                                                                                                                                      >  1:1255129/1‑150 (MQ=255)
              atttCATTTTTGCACCTGTTGAAAAATAGAATAGATATTTTGTATTTATTTGTAATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAata                                                                                                                                      <  1:197810/149‑1 (MQ=255)
                 tCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATAc                                                                                                                                   <  2:1597671/150‑1 (MQ=255)
                     ttttGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAAtt                                                                                                                               <  2:893273/150‑1 (MQ=255)
                      tttGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTg                                                                                                                              <  2:657702/150‑1 (MQ=255)
                      tttGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATACTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTg                                                                                                                              <  1:736867/150‑1 (MQ=255)
                         gCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTc                                                                                                                           <  2:864971/150‑1 (MQ=255)
                          cTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCg                                                                                                                          >  2:391024/1‑150 (MQ=255)
                           tCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGt                                                                                                                         >  1:444598/1‑150 (MQ=255)
                             cTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTg                                                                                                                       >  1:684749/1‑150 (MQ=255)
                             cTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTg                                                                                                                       <  1:214460/150‑1 (MQ=255)
                             cTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGACATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGttt                                                                                                                       >  2:1535561/1‑149 (MQ=255)
                              tGTTGCAAAATAGAATATATATTTTGAATTTATTTTTCATTAAACCATCAGGATGCTTGATATGTCATGCTATTTAATGTTCTCCATAATGATCATAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGa                                                                                                                      <  1:1169752/150‑1 (MQ=255)
                              tGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGa                                                                                                                      <  1:401945/150‑1 (MQ=255)
                              tGTTGCAAAATAGAAGAGATATTTTGAATTTAATTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTAACATAACAATATACAATTGCTCGTTGa                                                                                                                      <  1:267356/150‑1 (MQ=255)
                               gTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGC‑‑GATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAag                                                                                                                   >  2:1170652/1‑150 (MQ=255)
                                 tGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAag                                                                                                                   >  2:572311/1‑150 (MQ=255)
                                  gCAAAATAGAAGAGATATTTTGAATGTATTTGTGATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACTATTGATAGATGaaata                                                                                                                  >  1:248990/1‑148 (MQ=255)
                                   cAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAagag                                                                                                                 >  2:1767040/1‑150 (MQ=255)
                                       aTAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAg                                                                                                             >  1:460910/1‑150 (MQ=255)
                                        tAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCTAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGc                                                                                                            <  1:522173/150‑1 (MQ=255)
                                        tAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGc                                                                                                            <  1:724023/150‑1 (MQ=255)
                                         agaagaGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTATTGTTATCAATAATGAGAAAAATTCTGACAGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCt                                                                                                           >  1:66096/1‑150 (MQ=255)
                                            agagAAATATTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGAGATTTAATGTTCTCCATAATGAGAAAAATTCTGACCGGTGTAAGAACTTTCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTaaa                                                                                                        >  1:383089/1‑150 (MQ=255)
                                               gaTATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATc                                                                                                     >  2:41694/1‑150 (MQ=255)
                                               gaTATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATc                                                                                                     <  1:868799/150‑1 (MQ=255)
                                                atatTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATcc                                                                                                    <  1:1086878/150‑1 (MQ=255)
                                                  atTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCt                                                                                                  <  2:184483/150‑1 (MQ=255)
                                                   ttttGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCta                                                                                                 <  1:136778/150‑1 (MQ=255)
                                                    tttGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCtat                                                                                                >  2:1020283/1‑150 (MQ=255)
                                                     ttGAATTTATTTGTCAGTAAACCATCAGGATCCTTGATATGTCATGAGATTTAATGTTGTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTTAAATAGTGAGCTAAAATCCCtata                                                                                               <  2:83408/150‑1 (MQ=255)
                                                      tGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATaa                                                                                              <  1:170175/150‑1 (MQ=255)
                                                         atttatttTTCATTAAACCATCAGGACGCGTGATATGTCATGCGATTAAATGTTATCCATAATGAGCAAAATTCTGACCGGTGATAGCTCTTGCTTACATAACAATATTAAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAAAAg                                                                                           <  1:279389/150‑1 (MQ=255)
                                                          tttatttGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGt                                                                                          <  2:1020331/150‑1 (MQ=255)
                                                          tttatttGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAATTATGACCGGTTTAATCACTTGCTTACATAACAATATAAAATTTCTAGTTGAAAGAGTTAGCAAAAATACCTATAACAGTa                                                                                         >  1:1004323/1‑150 (MQ=255)
                                                           ttatttGTCATTAAACCATCAGGATGCGTGATATGTCATGAGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTa                                                                                         <  1:849595/150‑1 (MQ=255)
                                                             ttttGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGa                                                                                       <  1:1207561/149‑1 (MQ=255)
                                                                    aTTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATAAAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTcc                                                                                >  1:1147989/1‑150 (MQ=255)
                                                                      tAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCg                                                                              <  1:180575/150‑1 (MQ=255)
                                                                      tAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCg                                                                              <  1:391031/150‑1 (MQ=255)
                                                                           cATCAGGAAGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTTAGCTAAAATCCCTATAACAGTAGAACACTCCCGAGTGc                                                                         >  2:153392/1‑150 (MQ=255)
                                                                             tCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCgg                                                                       <  2:67033/150‑1 (MQ=255)
                                                                                ggATGAGTGATATGTCATGCGATTTAATGTTCTCCATAAGGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACTTAACAATATACAATTGCTCTTTGAAAGAGTGAGCTAAAATCCATATAACAGTAGATCCCTGCAGAGTGCTGAAg                                                                    >  1:271771/1‑150 (MQ=255)
                                                                                   tGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGt                                                                 <  2:791797/150‑1 (MQ=255)
                                                                                    gCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGtt                                                                >  1:1067143/1‑150 (MQ=255)
                                                                                       tGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGAc                                                             >  2:434485/1‑150 (MQ=255)
                                                                                       tGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGAc                                                             <  2:24532/150‑1 (MQ=255)
                                                                                       tGATATGTCATGCGATTTAATGTTCTCCATAATGAGAAAAATTCTGACCGGTGTAAGCACTTGCTTACATTACAATATACATTTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGAc                                                             <  2:210486/150‑1 (MQ=255)
                                                                                        gATATGTCATGCGATTTAATGTTCTCCCTAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACg                                                            <  2:1327408/150‑1 (MQ=255)
                                                                                        gATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACg                                                            <  2:1501893/150‑1 (MQ=255)
                                                                                               cATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATaga                                                     >  1:1584936/1‑150 (MQ=255)
                                                                                                 tGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAgg                                                   <  1:1816017/150‑1 (MQ=255)
                                                                                                  gCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGAAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGt                                                  >  2:17744/1‑150 (MQ=255)
                                                                                                      tttAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCa                                              <  2:466264/150‑1 (MQ=255)
                                                                                                               ctcCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCaaaa                                     >  2:1714212/1‑150 (MQ=255)
                                                                                                                 ccATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGt                                   >  2:561207/1‑150 (MQ=255)
                                                                                                                  cATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTg                                  <  1:1541901/150‑1 (MQ=255)
                                                                                                                      aTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaa                              >  1:1523201/1‑150 (MQ=255)
                                                                                                                      aTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaa                              <  1:1680500/150‑1 (MQ=255)
                                                                                                                       tGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaaa                             >  1:619424/1‑150 (MQ=255)
                                                                                                                       tGAGCAAAATTCTGACCGGTGTAAGCAATTTCTTACATAACAATATACAATTTCTAGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaaa                             <  1:598371/150‑1 (MQ=255)
                                                                                                                        gAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaaaa                            <  2:639557/150‑1 (MQ=255)
                                                                                                                         aGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaaaaa                           >  1:1825888/1‑150 (MQ=255)
                                                                                                                          gCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaaaaaa                          >  2:590756/1‑150 (MQ=255)
                                                                                                                          gCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaaaaaa                          >  1:324468/1‑150 (MQ=255)
                                                                                                                           cAAAATTCTGACCGGTGTAAGCACTTGCTTACATAATAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAac                         >  2:113821/1‑150 (MQ=255)
                                                                                                                                ttCTGACCGGTGTAAGAACTTGCTTACATAACAATATACAATTGCTCGTTTAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCATCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACAACTa                    >  1:1685589/1‑150 (MQ=255)
                                                                                                                                  cTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTAta                  >  2:1530828/1‑149 (MQ=255)
                                                                                                                                  cTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATg                  >  2:1533608/1‑150 (MQ=255)
                                                                                                                                   tGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGc                 >  2:1754238/1‑150 (MQ=255)
                                                                                                                                      gttatgtAAGCAAGTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCAt              >  1:991859/5‑150 (MQ=255)
                                                                                                                                      gttatgtAAGCAAGTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCAt              >  2:40694/5‑150 (MQ=255)
                                                                                                                                         gtgtAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATgcg           <  1:1095790/150‑1 (MQ=255)
                                                                                                                                           gtAAGCACTTGATTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAACTCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCt         <  2:248990/150‑1 (MQ=255)
                                                                                                                                            tAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTa        >  2:1031006/1‑150 (MQ=255)
                                                                                                                                            tAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGCTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTa        <  2:1486631/150‑1 (MQ=255)
                                                                                                                                             aaGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGAAGAAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTAt       >  1:1307186/1‑150 (MQ=255)
                                                                                                                                               gCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTATGt     <  2:684749/150‑1 (MQ=255)
                                                                                                                                                 aCTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTATGTTg   <  2:1257759/150‑1 (MQ=255)
                                                                                                                                                  cTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCACTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTATGTTGc  >  1:910856/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTCTTCCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTATGTTGC  >  NC_000913/230861‑231156

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