Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 146,733 G→A 8.1% intergenic (‑39/‑235) panD ← / → yadD aspartate 1‑decarboxylase/transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913146,7330GA 8.1% 343.0 / 12.1 113intergenic (‑39/‑235)panD/yadDaspartate 1‑decarboxylase/transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (50/53);  new base A (6/3);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

ACCGGCTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAACGC  >  NC_000913/146584‑146879
                                                                                                                                                     |                                                                                                                                                  
          tCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCCCGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCtt                                                                                                                                          <  2:87483/150‑1 (MQ=255)
           catgAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCATCGCCTGCttc                                                                                                                                         <  1:1569435/147‑2 (MQ=255)
           cAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCttt                                                                                                                                         <  2:1036614/150‑1 (MQ=255)
           cAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTATCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCttt                                                                                                                                         <  2:1447058/150‑1 (MQ=255)
            aaGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTg                                                                                                                                        <  2:277388/150‑1 (MQ=255)
              taaaaTCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAg                                                                                                                                      <  1:1529824/149‑1 (MQ=255)
              gAATATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAg                                                                                                                                      <  1:915471/150‑1 (MQ=255)
              gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAg                                                                                                                                      >  2:1534225/1‑150 (MQ=255)
              gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCATCGCCTGCttcgtt                                                                                                                                      <  1:348678/150‑5 (MQ=255)
              gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCATCGCCTGCttcgtt                                                                                                                                      <  1:1844362/150‑5 (MQ=255)
              gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGGTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAg                                                                                                                                      <  1:768871/150‑1 (MQ=255)
              caaaaTCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAg                                                                                                                                      <  1:1102899/149‑1 (MQ=255)
                   tCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATACCTTGTACCTTTGTACCCTGTCGGTAACGAAGCAGGCGATGCCTGCCTTGAGGAAAt                                                                                                                                 >  1:1234644/1‑150 (MQ=255)
                     cTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATtc                                                                                                                               <  1:636427/150‑1 (MQ=255)
                      tGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATtct                                                                                                                              <  2:1062068/150‑1 (MQ=255)
                      tGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTAACTTTCTACCCTGTCGTTAACTACGCAGGCGATGCCTGCGTTGAGGAAATtca                                                                                                                              >  1:182895/1‑149 (MQ=255)
                      tGGTCAATGGCGCAAGAACCTTCATAGTGCAGCTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTCCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGAGATGCCTGCTTTGAGGAAATtca                                                                                                                              >  1:915663/1‑149 (MQ=255)
                            aTGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCa                                                                                                                        >  2:1035198/1‑150 (MQ=255)
                              ggCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGt                                                                                                                      >  1:1602690/1‑150 (MQ=255)
                               gcgcAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGAAGCCTGCTTTGAGGAAATTCTCACGCAGTa                                                                                                                     >  1:1082211/1‑150 (MQ=255)
                                cgcAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTAt                                                                                                                    >  1:1331673/1‑150 (MQ=255)
                                 gcAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTAtt                                                                                                                   >  1:1506089/1‑150 (MQ=255)
                                  cAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTg                                                                                                                  <  2:1281289/150‑1 (MQ=255)
                                    aGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGcc                                                                                                                <  2:640725/150‑1 (MQ=255)
                                    aGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCGTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGcc                                                                                                                >  2:315725/1‑150 (MQ=255)
                                     gtaCCTTCATAGTGCAGGTCAGCATGAGTCAATTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGccc                                                                                                               <  1:1102586/148‑1 (MQ=255)
                                     gAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGccc                                                                                                               <  1:370981/150‑1 (MQ=255)
                                             aTAGTGCAGGGCCGCAAGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCGTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTg                                                                                                       <  2:1125905/150‑1 (MQ=255)
                                                gTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGttt                                                                                                    <  2:1270514/150‑1 (MQ=255)
                                                gTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGttt                                                                                                    <  2:1770881/150‑1 (MQ=255)
                                                    aGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAgtg                                                                                                <  2:933444/150‑1 (MQ=255)
                                                          gCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCGGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACt                                                                                          >  1:1076007/1‑150 (MQ=255)
                                                              gAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATc                                                                                      <  2:235531/150‑1 (MQ=255)
                                                               aGTCACTTTCACGCGGTGGAGTTTGCCCTGTAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCt                                                                                     >  2:960981/1‑150 (MQ=255)
                                                                gTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTg                                                                                    >  2:834172/1‑150 (MQ=255)
                                                                  cACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGAc                                                                                  >  2:1530150/1‑150 (MQ=255)
                                                                   aCTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACg                                                                                 <  1:800141/150‑1 (MQ=255)
                                                                      ttCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCa                                                                              >  1:1048147/1‑150 (MQ=255)
                                                                        cACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTTACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCAtt                                                                            >  1:522692/1‑150 (MQ=255)
                                                                        cACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCAtt                                                                            >  2:1080154/1‑150 (MQ=255)
                                                                          cgcgGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTg                                                                          <  2:1851841/150‑1 (MQ=255)
                                                                            cgGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTGACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCg                                                                        >  1:329271/1‑150 (MQ=255)
                                                                            cgGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGGCGTTAACGAAGCAGGCGTTGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCg                                                                        <  2:97698/150‑1 (MQ=255)
                                                                              gTGGATGTTGCCCTGCAGCATCGTGCGTATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGAATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCACTTGCggc                                                                      <  1:725415/150‑3 (MQ=255)
                                                                                ggAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGGCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTcag                                                                    <  2:551092/150‑1 (MQ=255)
                                                                                  aGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTcagca                                                                  >  1:416576/1‑150 (MQ=255)
                                                                                      tGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTAAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGttt                                                              <  1:645594/150‑1 (MQ=255)
                                                                                        cccTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGACCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGttttc                                                            >  1:291695/1‑148 (MQ=255)
                                                                                        cccTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATACTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGc                                                            >  2:1634966/1‑150 (MQ=255)
                                                                                          cTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGTCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGGCAGCAGTTTGCGt                                                          <  1:1010133/150‑1 (MQ=255)
                                                                                           tGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTCAGGATATTCTCACGCAGTAATGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTa                                                         >  1:892630/1‑150 (MQ=255)
                                                                                            gcagcaTCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCTCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTAc                                                        >  1:252965/1‑150 (MQ=255)
                                                                                              agcaTCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACcg                                                      >  2:1057030/1‑150 (MQ=255)
                                                                                                caTCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACcgcg                                                    >  2:448770/1‑150 (MQ=255)
                                                                                                  tCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGcc                                                  <  1:1572533/150‑1 (MQ=255)
                                                                                                    gTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGa                                                >  1:552297/1‑150 (MQ=255)
                                                                                                     tGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGaa                                               >  1:53363/1‑150 (MQ=255)
                                                                                                      gCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAg                                              <  1:646186/150‑1 (MQ=255)
                                                                                                       cGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGAAGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGc                                             >  1:645205/1‑150 (MQ=255)
                                                                                                         aaTCATAACTTCTACCTTTCTACCCTGTCGGTAAAGAAGCAGGCGATGCCTGATTTGAGGAAATTCTCACGCAGTATTGCCAGATTTTGGTTTAGTGTCTACTCAGCTGACGGCATTTGCGTCAGCAGTTGGCGTACCGCGCCGCAGCgg                                           >  2:1607587/1‑149 (MQ=255)
                                                                                                           tCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGc                                         <  1:358674/150‑1 (MQ=255)
                                                                                                           tCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGc                                         <  1:1285884/150‑1 (MQ=255)
                                                                                                           tCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGc                                         >  2:1466368/1‑150 (MQ=255)
                                                                                                             aTAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGTAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                       >  2:90924/1‑150 (MQ=255)
                                                                                                             aTAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCTTACCGCGCCGAAGCGTGCCa                                       >  1:1664112/1‑150 (MQ=255)
                                                                                                             aTAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTTCCGCGCCGAAGCGTGCCa                                       >  1:985991/1‑150 (MQ=255)
                                                                                                               aaCTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATg                                     <  1:1534225/150‑1 (MQ=255)
                                                                                                               aaCTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTGTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATg                                     <  2:372045/150‑1 (MQ=255)
                                                                                                                 cTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTa                                   >  1:1434707/1‑150 (MQ=255)
                                                                                                                   tCTACCTTTCTACCCTGTCTTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTATAGTGTCTACTCATCTGACGGCATTTGCGTAAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATc                                 <  2:1156371/150‑1 (MQ=255)
                                                                                                                   tCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGCCGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATc                                 >  1:1185996/1‑150 (MQ=255)
                                                                                                                     tACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGCAGCGGGCCATGTATCGc                               >  2:902128/1‑150 (MQ=255)
                                                                                                                     tACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGc                               >  1:987966/1‑150 (MQ=255)
                                                                                                                         tttCTACCCTGTCGTTACCGCCGCCGGCGATGCCTGCTTTGAGGACTTCCGCCCGCTCTATAGCACGATTTTTGTTTAGTGTCTACTCATCTGAAGGCATTTGCGTCACCAGTTTGCGTACCGCGCCGAAGCGTGCCATGAAGCGCCTGt                           <  2:294190/150‑1 (MQ=255)
                                                                                                                          ttCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTa                          >  1:1780946/1‑150 (MQ=255)
                                                                                                                            cTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAAc                        <  1:841838/150‑1 (MQ=255)
                                                                                                                             tACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACg                       <  2:529617/150‑1 (MQ=255)
                                                                                                                               cccTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGaa                     >  2:1835858/1‑150 (MQ=255)
                                                                                                                               cccTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTATCACGCAGTATTCCCCGATTATTGTTTAGTGACTACTCATCTGACGGCATTTGCGTCAGCAGTTTGAGTACCGCGCCGAAGAGTGCCATGAATCGCCTGTAACaaa                     >  1:1190477/1‑150 (MQ=255)
                                                                                                                                 cTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAAtt                   <  1:1152076/150‑1 (MQ=255)
                                                                                                                                    tCGTAAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGat                <  2:1518314/150‑1 (MQ=255)
                                                                                                                                     cGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGata               >  1:1781107/1‑150 (MQ=255)
                                                                                                                                        cctcaaaGCAGGCATCGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTAATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGCAGCGTGCCATGTATCGCCTGTAACGAATTGATtttt            >  2:477612/6‑150 (MQ=255)
                                                                                                                                        cctcaaaGCAGGCATCGCCTGCTTCGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCCTCTGACGGCATTTGCGTCCGCAGTTTGCGTCCCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATAttt            >  1:1654303/6‑150 (MQ=255)
                                                                                                                                         ctcaaaGCAGGCATCGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtc           >  1:1317091/5‑150 (MQ=255)
                                                                                                                                         aaCGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtc           <  2:1870570/150‑1 (MQ=255)
                                                                                                                                         aaCGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtc           <  1:1841915/150‑1 (MQ=255)
                                                                                                                                           caaaGCAGGCATCGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtctc         >  2:524800/3‑150 (MQ=255)
                                                                                                                                           cGAAGCACGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTtctc         <  2:567335/150‑1 (MQ=255)
                                                                                                                                            gAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCCCGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTctct        <  2:1138491/150‑1 (MQ=255)
                                                                                                                                            gAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTctct        <  1:960981/150‑1 (MQ=255)
                                                                                                                                            gAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTctct        >  2:27343/1‑150 (MQ=255)
                                                                                                                                            gAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTctct        <  1:1792963/150‑1 (MQ=255)
                                                                                                                                            gAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTctct        <  1:1031956/150‑1 (MQ=255)
                                                                                                                                            aaaGCAGGCATCGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTTCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTctct        >  2:1386169/2‑150 (MQ=255)
                                                                                                                                            aaaGCAGGCATCGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTctct        >  2:1802312/2‑150 (MQ=255)
                                                                                                                                               gCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGaa     >  2:1202234/1‑150 (MQ=255)
                                                                                                                                               gCAGGCGATGCCTGATTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGaa     <  2:329271/150‑1 (MQ=255)
                                                                                                                                                ccggCGATGCCTGCGGTGAGGAAAGGCTCACGCAGTATTGCCCGATTCTTGTTTAGAGTCTACTCATCTGACGGCATTTGCGGCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGACTTGATATTTCTCTGAAc    <  2:142602/148‑1 (MQ=255)
                                                                                                                                                  ggCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTTAACgc  >  2:1205485/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACCGGCTGCGTCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAACGC  >  NC_000913/146584‑146879

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