Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,525,908 C→G 5.9% intergenic (+17/+22) 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,9080CG 5.9% 348.0 / 9.8 114intergenic (+17/+22)yfeN/yfeRputative outer membrane protein/transcriptional regulator of yefH
Reads supporting (aligned to +/- strand):  ref base C (49/56);  new base G (3/4);  total (53/61)
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.21e-01

Only 100 of 114 total aligned reads displayed.

GCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAC  >  NC_000913/2525761‑2526056
                                                                                                                                                   |                                                                                                                                                     
gcgcAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTc                                                                                                                                                     >  2:1453547/1‑150 (MQ=255)
   cAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCg                                                                                                                                                  <  1:1663591/150‑1 (MQ=255)
   cAGATGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATGCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCg                                                                                                                                                  <  1:109561/150‑1 (MQ=255)
     gTTGCACTGTGGTGGAATGCAACATCACACATTACTCCGGGGAATCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGc                                                                                                                                                <  1:1593110/150‑1 (MQ=255)
        gCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGt                                                                                                                                             <  1:662922/150‑1 (MQ=255)
        gCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTATTACCAGGATGCGATAATCTACTCCCTCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGt                                                                                                                                             <  1:1067547/150‑1 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGTCTAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGTaat                                                                                                                                          <  1:998035/150‑3 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATg                                                                                                                                          <  1:229379/150‑1 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATg                                                                                                                                          <  1:1475852/150‑1 (MQ=255)
           cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTAAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGTaat                                                                                                                                          <  1:62029/150‑3 (MQ=255)
             gtggtgGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGaa                                                                                                                                        >  1:1015346/1‑150 (MQ=255)
             gtggtgGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGTaatca                                                                                                                                        <  1:688035/150‑5 (MQ=255)
              tggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAAt                                                                                                                                       >  1:816477/1‑150 (MQ=255)
              gggtggAATGCAAAATCAGACATTACTACGGCGATTCAGTATCGATATTCGGATGACTAACTGGGCGTAGTATTTTACCAGTATGCGTTATTCTTCTCCATCAAATTTAACTTCAGAGAGATTACTGCCGGGGCTCCCGGCAGTATGAAt                                                                                                                                       <  2:1151034/149‑1 (MQ=255)
               ggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAAtt                                                                                                                                      <  2:478661/150‑1 (MQ=255)
                gtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTc                                                                                                                                     <  2:468550/150‑1 (MQ=255)
                gtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTc                                                                                                                                     >  1:1591817/1‑150 (MQ=255)
                   gAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCCTCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCTGCAGTATGAAATCAct                                                                                                                                  >  1:1570545/1‑150 (MQ=255)
                   gAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCAct                                                                                                                                  <  2:491043/150‑1 (MQ=255)
                        taaCATCACACATTACTTCGGTGATTCAGTATGGATATGCGGATGACTAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCTTCAAATTTAACTTCTGATTGATTACTGCCGGGGCTGCCGGCAGTATGAATTCACTCtata                                                                                                                             <  2:1027849/149‑1 (MQ=255)
                        cAACATCACACATTACTACGGGGATTCAGTCTCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCtata                                                                                                                             <  2:691698/150‑1 (MQ=255)
                        cAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCtata                                                                                                                             >  2:571811/1‑150 (MQ=255)
                        cAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCtata                                                                                                                             >  2:1274719/1‑150 (MQ=255)
                          aCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGAGGCTCCCGGCAGTATGAATTCACTCTatata                                                                                                                           >  1:1049574/1‑149 (MQ=255)
                           cATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCt                                                                                                                          <  2:1689729/150‑1 (MQ=255)
                                cacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATAc                                                                                                                     >  2:342834/1‑150 (MQ=255)
                                cacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATAc                                                                                                                     <  1:1126837/150‑1 (MQ=255)
                                   aTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAAc                                                                                                                  <  2:547194/150‑1 (MQ=255)
                                    ttACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCAATCAAATATAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTAACTCTATATCTGATACAACg                                                                                                                 >  1:700328/1‑150 (MQ=255)
                                       ctacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGAt                                                                                                              <  1:518465/150‑1 (MQ=255)
                                        tacttggATTCAGTATCGATATGCGGATGACAAACTGTGCTAAGATTTTTACCAGGATGCTATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGTGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATc                                                                                                            <  2:600149/145‑1 (MQ=255)
                                        tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATc                                                                                                            <  1:176841/150‑1 (MQ=255)
                                           ggggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATcccc                                                                                                         >  2:1571776/1‑150 (MQ=255)
                                            gggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCt                                                                                                        >  2:1492317/1‑150 (MQ=255)
                                            gggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCAGGGGATCCCGGAAGTATGAATTCACTCTATATATGATTCAACGTATcccca                                                                                                        >  1:222238/1‑149 (MQ=255)
                                              gATTCAGTATCGATATGCTGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTc                                                                                                      <  1:1282772/150‑1 (MQ=255)
                                              gATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTc                                                                                                      >  1:1046335/1‑150 (MQ=255)
                                                 tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCg                                                                                                   <  2:380583/150‑1 (MQ=255)
                                                 tCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCg                                                                                                   <  1:692594/150‑1 (MQ=255)
                                                   aGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGcc                                                                                                 >  2:11064/1‑150 (MQ=255)
                                                   aGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGcc                                                                                                 >  1:1530264/1‑150 (MQ=255)
                                                      aTCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAg                                                                                              <  1:994201/150‑1 (MQ=255)
                                                       tCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGc                                                                                             >  2:147751/1‑150 (MQ=255)
                                                        cGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGcc                                                                                            <  1:1183993/150‑1 (MQ=255)
                                                         gaaatGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCa                                                                                           <  2:1249094/147‑1 (MQ=255)
                                                         gATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCa                                                                                           >  2:909134/1‑150 (MQ=255)
                                                         gATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCa                                                                                           <  1:490456/150‑1 (MQ=255)
                                                           tatGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGATCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCAtt                                                                                         <  2:1203714/150‑1 (MQ=255)
                                                                ggATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTTCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGTAGCCCCTTCCCGCCCAGCCATTAgcgc                                                                                    >  1:90381/1‑150 (MQ=255)
                                                                    gACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTg                                                                                <  2:1687139/150‑1 (MQ=255)
                                                                    gACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTg                                                                                >  1:1693619/1‑150 (MQ=255)
                                                                     aCAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCACTTCCCGCCCAGCCATTCGCGCATTGc                                                                               >  2:174505/1‑150 (MQ=255)
                                                                       aaaCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                                                             <  2:900200/150‑1 (MQ=255)
                                                                            gggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAACTGTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtc                                                                        <  1:865219/150‑1 (MQ=255)
                                                                             ggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtct                                                                       <  1:681911/150‑1 (MQ=255)
                                                                             ggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtct                                                                       >  2:773115/1‑150 (MQ=255)
                                                                             ggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtct                                                                       >  1:521055/1‑150 (MQ=255)
                                                                              gCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCATATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtct                                                                       <  2:1681487/150‑1 (MQ=255)
                                                                              gCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGtctc                                                                      >  2:839798/1‑150 (MQ=255)
                                                                               cGAAGATTTTTACGAGGATGCGATAATCTACTCCATCAATTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCGTTGCCCTGGTCTCg                                                                     >  2:213175/1‑150 (MQ=255)
                                                                                 aaGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGAATGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTa                                                                   <  2:1094234/150‑1 (MQ=255)
                                                                                   gATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTAcc                                                                 >  2:1676446/1‑150 (MQ=255)
                                                                                   gATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTAcc                                                                 >  2:1533521/1‑150 (MQ=255)
                                                                                   gATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTACCGCCCAGCCATTAGCGCATTGACCTGGTCTCGTAcc                                                                 >  1:297969/1‑150 (MQ=255)
                                                                                      ttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACg                                                              <  1:1441266/150‑1 (MQ=255)
                                                                                         tACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTcc                                                           >  2:33776/1‑150 (MQ=255)
                                                                                          aCCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTccc                                                          <  2:557293/150‑1 (MQ=255)
                                                                                           ccAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCa                                                         <  1:1584052/150‑1 (MQ=255)
                                                                                           ccAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCa                                                         <  1:773115/150‑1 (MQ=255)
                                                                                           ccAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCa                                                         <  1:1345467/150‑1 (MQ=255)
                                                                                           ccAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCa                                                         <  1:405056/150‑1 (MQ=255)
                                                                                           ccAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCa                                                         <  2:1121232/150‑1 (MQ=255)
                                                                                                       aaTCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgc                                             <  1:1492317/150‑1 (MQ=255)
                                                                                                               ccATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTgg                                     <  2:1312651/150‑1 (MQ=255)
                                                                                                               ccATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTgg                                     >  2:1187149/1‑150 (MQ=255)
                                                                                                               ccATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTgg                                     >  1:1580464/1‑150 (MQ=255)
                                                                                                                 aTCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaa                                   >  2:514532/1‑150 (MQ=255)
                                                                                                                  tCAAAATTAACTTCAGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  <  1:880666/150‑1 (MQ=255)
                                                                                                                     aaTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTg                               <  2:611512/150‑1 (MQ=255)
                                                                                                                      aTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCCTTGCCCTGGTCTCGTACCACGTCCCATAAAGCGTCAGCGGCGGTGGAAAGTGa                              >  2:1040305/1‑150 (MQ=255)
                                                                                                                      aTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCACCTCCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGa                              <  1:56740/150‑1 (MQ=255)
                                                                                                                       tttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGaa                             >  2:343462/1‑150 (MQ=255)
                                                                                                                        ttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAAc                            <  2:1591817/150‑1 (MQ=255)
                                                                                                                                 tcttgattACTGCCGGGGCTCCATGCAGTATTAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCTGTGGAAAGTGAACGGTTCTTcc                   <  1:1142378/148‑1 (MQ=255)
                                                                                                                                  attgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCt                  >  1:789996/1‑150 (MQ=255)
                                                                                                                                  attgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCAAGCCATTAGCGCATTGCCCTGGTCTCGTACAACTTCCCATAAAGCTTCAGCGGCGGGGGAAAGTGAACGGTTCTTCCt                  >  2:862165/1‑150 (MQ=255)
                                                                                                                                     tcatACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGGCCCAGAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg               >  1:718373/4‑150 (MQ=255)
                                                                                                                                     gattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg               <  1:489503/150‑1 (MQ=255)
                                                                                                                                     gattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg               >  2:1676703/1‑150 (MQ=255)
                                                                                                                                     gattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGGGGAAAGTGAACGGTTCTTCCTgcg               >  1:812999/1‑150 (MQ=255)
                                                                                                                                     gattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCAAAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg               >  2:1535829/1‑150 (MQ=255)
                                                                                                                                       atACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATTCAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCACTGCAAAGTGCTCGGTTCTGCCTGCGCa             >  1:230996/2‑150 (MQ=255)
                                                                                                                                       atACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCCGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcgcc             >  2:1009359/2‑149 (MQ=255)
                                                                                                                                         aCTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCAcc           <  1:11064/150‑1 (MQ=255)
                                                                                                                                          cTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGGACGGTTCTTCCTGCGCACCa          >  1:475116/1‑150 (MQ=255)
                                                                                                                                            gCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGc        >  1:171079/1‑150 (MQ=255)
                                                                                                                                             ccGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGTATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCa       <  2:911358/150‑1 (MQ=255)
                                                                                                                                              cGGGGCTCCCGGCTGTATGAATTCACTCTATATCTGATACAACGGATCGCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCAt      <  2:586422/150‑1 (MQ=255)
                                                                                                                                                gggCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCa    >  2:1155113/1‑150 (MQ=255)
                                                                                                                                                  gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAc  <  1:142689/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAC  >  NC_000913/2525761‑2526056

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: