Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,525,907 G→A 5.6% intergenic (+16/+23) yfeN → / ← yfeR putative outer membrane protein/transcriptional regulator of yefH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,525,9070GA 5.6% 345.7 / 4.2 107intergenic (+16/+23)yfeN/yfeRputative outer membrane protein/transcriptional regulator of yefH
Reads supporting (aligned to +/- strand):  ref base G (50/50);  new base A (3/3);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

Only 100 of 108 total aligned reads displayed.

GGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATC  >  NC_000913/2525760‑2526053
                                                                                                                                                   |                                                                                                                                                  
ggCGCAGTTGCACTGTTGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTAACAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTTCCGGGGCt                                                                                                                                                  >  1:1215921/1‑150 (MQ=255)
ggCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCt                                                                                                                                                  >  1:902376/1‑150 (MQ=255)
ggCGCAGTTGCACTGTGGTGGAATGCAACATCACACAATACTACGGGGATTCAGAATCGATATGCGGATGACAAACTGGGCGAAGATTATTAACAGGATGCGCTAATCAACACCATCAAATTTAACTTCAGATTGATTACTGCCGGGGCt                                                                                                                                                  >  1:1215925/1‑150 (MQ=255)
    cAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCg                                                                                                                                              <  1:594543/150‑1 (MQ=255)
      gTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGc                                                                                                                                            <  2:1374196/150‑1 (MQ=255)
        tGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAg                                                                                                                                          >  2:1111576/1‑150 (MQ=255)
        tGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAg                                                                                                                                          >  1:1593853/1‑150 (MQ=255)
         gCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGt                                                                                                                                         >  1:836497/1‑150 (MQ=255)
         gCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGt                                                                                                                                         <  2:189584/150‑1 (MQ=255)
         gCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGt                                                                                                                                         <  2:53637/150‑1 (MQ=255)
          cACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGTa                                                                                                                                        <  1:1245584/150‑1 (MQ=255)
            cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGAGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATg                                                                                                                                      >  1:757878/1‑150 (MQ=255)
             tgtgTTTGAATGCAACATCACACATTACTACGGTTATTCAGTATCGATATGCGGATGACAAACTGTTCTAAGATTATTACAAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGTaatc                                                                                                                                     <  1:695204/150‑4 (MQ=255)
               tggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAAt                                                                                                                                   <  1:207861/150‑1 (MQ=255)
                ggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAAtt                                                                                                                                  >  2:501520/1‑150 (MQ=255)
                 gtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTc                                                                                                                                 <  2:1678032/150‑1 (MQ=255)
                    gAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCAct                                                                                                                              <  2:1575766/150‑1 (MQ=255)
                        gCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCtat                                                                                                                          >  1:780707/1‑150 (MQ=255)
                        gCAACATCACACATTACTACCGGGATTCAGTATCGATAAGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATAAAAATTAACTTCTGATTGATTACTTCCGGGGCTCCCGGCAGTATGAATTAAATCtaa                                                                                                                          >  2:724610/1‑149 (MQ=255)
                         cAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCtata                                                                                                                         >  1:95332/1‑150 (MQ=255)
                           aCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATc                                                                                                                       <  1:1716530/150‑1 (MQ=255)
                            cATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCt                                                                                                                      <  1:983942/150‑1 (MQ=255)
                            aaTCACACATTACTACGGTGATTCAGTATCTATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCt                                                                                                                      <  2:129038/149‑1 (MQ=255)
                             aTCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTg                                                                                                                     >  2:924921/1‑150 (MQ=255)
                               cacacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGAt                                                                                                                   >  1:1221964/1‑150 (MQ=255)
                                acacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCTGGCAGTATGAATTCACTCTATATCttcta                                                                                                                  >  2:684042/1‑146 (MQ=255)
                                acacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATa                                                                                                                  >  2:605606/1‑150 (MQ=255)
                                   caTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATTAATTCACTCTATATCTGATACaa                                                                                                               >  1:634371/1‑150 (MQ=255)
                                      tactacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACgg                                                                                                            <  1:1529491/150‑1 (MQ=255)
                                      tactacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACgg                                                                                                            <  1:787117/150‑1 (MQ=255)
                                        ctacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGAt                                                                                                          >  1:586144/1‑150 (MQ=255)
                                          acGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATcc                                                                                                        <  2:1456434/150‑1 (MQ=255)
                                            ggggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATcccc                                                                                                      <  2:925574/150‑1 (MQ=255)
                                             gggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCt                                                                                                     <  2:1270992/150‑1 (MQ=255)
                                              ggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCTATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCtt                                                                                                    <  1:856030/150‑1 (MQ=255)
                                              ggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCtt                                                                                                    <  2:289752/150‑1 (MQ=255)
                                               gATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTc                                                                                                   <  2:1592824/150‑1 (MQ=255)
                                                aTTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTcc                                                                                                  <  2:1130537/150‑1 (MQ=255)
                                                  tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATGCAACGGATCCCCTTCCCg                                                                                                >  2:423385/1‑150 (MQ=255)
                                                   cAGTATCAATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGc                                                                                               >  2:331944/1‑150 (MQ=255)
                                                   cAGAATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCAGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGc                                                                                               >  2:1152021/1‑150 (MQ=255)
                                                    aGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGcc                                                                                              >  2:1214688/1‑150 (MQ=255)
                                                     gTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGccc                                                                                             <  1:417456/150‑1 (MQ=255)
                                                      tATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCa                                                                                            <  1:501520/150‑1 (MQ=255)
                                                       aTCGATATGCGGATGACAAACTGGGAGAAGATTATTACCAGGATGCGATAATCTACTCCATCAAATATAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAg                                                                                           <  2:255671/150‑1 (MQ=255)
                                                          gATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCa                                                                                        >  1:1433392/1‑150 (MQ=255)
                                                            tatGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCAtt                                                                                      <  2:783481/150‑1 (MQ=255)
                                                            tatGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCAtt                                                                                      >  1:1652558/1‑150 (MQ=255)
                                                                 ggATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAgcgc                                                                                 >  2:1401411/1‑150 (MQ=255)
                                                                 ggATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAgcgc                                                                                 >  2:1597816/1‑150 (MQ=255)
                                                                   aTGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCAt                                                                               <  1:331944/150‑1 (MQ=255)
                                                                          aCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTg                                                                        <  2:1690682/150‑1 (MQ=255)
                                                                               gCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACGCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtctc                                                                   >  2:657597/1‑150 (MQ=255)
                                                                                cGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCg                                                                  >  2:470626/1‑150 (MQ=255)
                                                                                cGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCg                                                                  >  2:1634830/1‑150 (MQ=255)
                                                                                  aaGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGTTTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTa                                                                <  2:609323/150‑1 (MQ=255)
                                                                                    gATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTAcc                                                              >  2:1492539/1‑150 (MQ=255)
                                                                                     aTTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCTGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCAATAGCGCATTGCCATGGTCTCGTACCa                                                             >  1:1084578/1‑150 (MQ=255)
                                                                                         ttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTc                                                         <  2:505462/150‑1 (MQ=255)
                                                                                          tACCAGGATGCGATAATCTACTCCATCAAATTTACCTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTcc                                                        <  1:727490/150‑1 (MQ=255)
                                                                                           aCCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTccc                                                       >  1:442803/1‑150 (MQ=255)
                                                                                            ccAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCa                                                      >  2:71832/1‑150 (MQ=255)
                                                                                              aGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCATCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATa                                                    <  1:1152021/150‑1 (MQ=255)
                                                                                               ggATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCAGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATaa                                                   <  1:1325888/150‑1 (MQ=255)
                                                                                                 aTGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAg                                                 <  2:1380948/150‑1 (MQ=255)
                                                                                                     gATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTc                                             <  1:605606/150‑1 (MQ=255)
                                                                                                      aTAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCGCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCa                                            <  2:149142/150‑1 (MQ=255)
                                                                                                         aTCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcg                                         >  2:1625999/1‑150 (MQ=255)
                                                                                                              cTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGt                                    <  2:787743/150‑1 (MQ=255)
                                                                                                               tCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTg                                   <  1:71832/150‑1 (MQ=255)
                                                                                                                ccATCAAATTTAACTTCTGATTGATTTCTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACAGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTgg                                  >  2:1311710/1‑150 (MQ=255)
                                                                                                                 cATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGCTCTCGTACCACGTCCCAAAAAGCTTCAGCGGCggagga                                 >  1:1104700/1‑150 (MQ=255)
                                                                                                                   tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                               >  1:478154/1‑150 (MQ=255)
                                                                                                                    aaaaTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATTAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAg                              <  1:1570682/149‑1 (MQ=255)
                                                                                                                     aaaTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGt                             >  2:1383457/1‑150 (MQ=255)
                                                                                                                       aTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGa                           <  2:1593853/150‑1 (MQ=255)
                                                                                                                         ttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATGCAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAAc                         <  1:423385/150‑1 (MQ=255)
                                                                                                                          tAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACg                        <  1:470626/150‑1 (MQ=255)
                                                                                                                           aaCTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACgg                       >  2:859093/1‑150 (MQ=255)
                                                                                                                            aCTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGt                      <  1:924921/150‑1 (MQ=255)
                                                                                                                            aCTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGt                      >  1:311979/1‑150 (MQ=255)
                                                                                                                             cTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                     >  1:521027/1‑150 (MQ=255)
                                                                                                                               tCTGTTTGTTTACTGCCGGGGCTCCCGGCAGTATGTATTCACTCTATATGTGATACTACGGATCCCCTTCCCGCCAAGAAATTATCGCATTTCCCTGGTCTCTTACCACGTCCCATAAAGCTTCTGCGGCTGTGGAAAGTGAACTGttct                   <  2:1215925/150‑1 (MQ=255)
                                                                                                                               tCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTTGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtgct                   <  2:1215921/150‑1 (MQ=255)
                                                                                                                                  gattgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTcc                <  1:1111576/150‑1 (MQ=255)
                                                                                                                                   attgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCt               <  1:1625999/150‑1 (MQ=255)
                                                                                                                                   attgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCACTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGGACCACGTCCAATAAAGCATCAGCGGCGGTGGAAAGTGAACGGTTCTTCCt               >  1:956560/1‑150 (MQ=255)
                                                                                                                                     tgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCCGCCATTAGCGCATTGCCCTGGGCTCGTACCCCGTCCCAGAACGCTTCAGCGGCGGTGGAAAGGGCACGGTTCTTCCTgc             >  1:174759/1‑150 (MQ=255)
                                                                                                                                     tgattACTGCCGGGGCTCCCGGCAGTATGAATTCAATCTATATCTGATAAAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCGTGGTCTCGTACCACGACCCATATAGCTTCAGCGGCGGTGGAAAGTGAAAGGTTCTTCCTgc             >  1:189526/1‑150 (MQ=255)
                                                                                                                                      gattACTGCCGTGGCTCCCGGCAGTATGAATTCACTCTTTATCTGATACAACGGATCCCCTTCCCGCCCATCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgct            <  2:757878/150‑2 (MQ=255)
                                                                                                                                      gattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg            >  2:600151/1‑150 (MQ=255)
                                                                                                                                       catACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcgc           >  1:806199/3‑150 (MQ=255)
                                                                                                                                        atACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCa          >  1:1282533/2‑150 (MQ=255)
                                                                                                                                          aCTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCAcc        >  2:359431/1‑150 (MQ=255)
                                                                                                                                           cTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCAccc       >  2:513518/1‑149 (MQ=255)
                                                                                                                                             gCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGc     >  2:6617/1‑150 (MQ=255)
                                                                                                                                             gCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGc     >  2:834383/1‑150 (MQ=255)
                                                                                                                                              ccGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCa    <  1:1419452/150‑1 (MQ=255)
                                                                                                                                               cGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCCTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCAt   <  1:1214688/150‑1 (MQ=255)
                                                                                                                                                ggggCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATc  >  1:1649539/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATC  >  NC_000913/2525760‑2526053

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: