Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 231,008 T→G 6.3% intergenic (‑127/‑114) 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,0080TG 6.3% 349.0 / 10.9 112intergenic (‑127/‑114)yafC/yafDLysR family putative transcriptional regulator/endo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base T (52/53);  new base G (2/5);  total (54/58)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 113 total aligned reads displayed.

TTCTTCCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTATGTTGCC  >  NC_000913/230861‑231157
                                                                                                                                                   |                                                                                                                                                     
ttcttcCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTg                                                                                                                                                     <  1:1725050/150‑1 (MQ=255)
  cttcCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCt                                                                                                                                                   >  2:1025555/1‑150 (MQ=255)
  cttcCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCt                                                                                                                                                   >  1:1706834/1‑150 (MQ=255)
   ttcCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCTATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCtt                                                                                                                                                  <  1:1249412/150‑1 (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:386385/150‑1 (MQ=255)
     ccGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTAc                                                                                                                                                >  2:1541901/1‑150 (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)
          gTGGCTTTCATTTTTGCTCATGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAAc                                                                                                                                           <  2:135075/150‑1 (MQ=255)
            ggCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACaa                                                                                                                                         <  1:1170822/150‑1 (MQ=255)
             gCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGAACTTGCTTACATAACAat                                                                                                                                        >  2:613271/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                                                                                                                        >  2:1783797/1‑150 (MQ=255)
                             cTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTg                                                                                                                        <  1:214460/150‑1 (MQ=255)
                             cTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTg                                                                                                                        >  1:684749/1‑150 (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)
                                        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)
                                                 tatTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTAGTTGAAAGAGTGAGCTAAAATccc                                                                                                    >  1:651416/1‑150 (MQ=255)
                                                  atTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCt                                                                                                   <  2:184483/150‑1 (MQ=255)
                                                   ttttGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCta                                                                                                  <  1:136778/150‑1 (MQ=255)
                                                    tttGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCtat                                                                                                 >  2:1020283/1‑150 (MQ=255)
                                                    tttGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCtat                                                                                                 >  2:1430824/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)
                                                             ttttGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGa                                                                                        <  1:1207561/149‑1 (MQ=255)
                                                             atttGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGa                                                                                        >  2:66027/1‑150 (MQ=255)
                                                                    aTTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATAAAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTcc                                                                                 >  1:1147989/1‑150 (MQ=255)
                                                                      tAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCg                                                                               <  1:391031/150‑1 (MQ=255)
                                                                      tAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCg                                                                               <  1:180575/150‑1 (MQ=255)
                                                                           cATCAGGAAGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTTAGCTAAAATCCCTATAACAGTAGAACACTCCCGAGTGc                                                                          >  2:153392/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)
                                                                                        gATATGTCATGCGATTTAATGTTCTCCCTAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACg                                                             <  2:1327408/150‑1 (MQ=255)
                                                                                        gATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACg                                                             <  2:1501893/150‑1 (MQ=255)
                                                                                                 tGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAgg                                                    <  1:1816017/150‑1 (MQ=255)
                                                                                                  gCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGAAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGt                                                   >  2:17744/1‑150 (MQ=255)
                                                                                                   cGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGtt                                                  <  1:378816/150‑1 (MQ=255)
                                                                                                            gTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACATAATAGAGGTTTCAAAGTCa                                         >  2:1560284/1‑150 (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)
                                                                                                                          gCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaaaaaa                           >  2:590756/1‑150 (MQ=255)
                                                                                                                          gCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaaaaaa                           >  1:324468/1‑150 (MQ=255)
                                                                                                                           cAAAATTCTGACCGGTGTAAGCACTTGCTTACATAATAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAac                          >  2:113821/1‑150 (MQ=255)
                                                                                                                            aaaaTTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAaca                         <  1:992017/150‑1 (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:444598/150‑1 (MQ=255)
                                                                                                                                   tGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGc                  >  2:1754238/1‑150 (MQ=255)
                                                                                                                                      gttatgtAAGCAAGTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCAt               >  2:40694/5‑150 (MQ=255)
                                                                                                                                      gttatgtAAGCAAGTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCAt               >  1:991859/5‑150 (MQ=255)
                                                                                                                                       cGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATg              >  1:18379/1‑150 (MQ=255)
                                                                                                                                         gtgtAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATgcg            <  1:1095790/150‑1 (MQ=255)
                                                                                                                                           gtAAGCACTTGATTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAACTCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCt          <  2:248990/150‑1 (MQ=255)
                                                                                                                                            tAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGCTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTa         <  2:1486631/150‑1 (MQ=255)
                                                                                                                                             aaGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGAAGAAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTAt        >  1:1307186/1‑150 (MQ=255)
                                                                                                                                              aGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGAAGAACCCTCCCGAGTGCGGATGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTATg       >  1:1660839/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)
                                                                                                                                                   ttGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTATGTTGcc  <  1:66027/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTCTTCCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTATGTTGCC  >  NC_000913/230861‑231157

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: