Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 146,734 A→T 8.8% intergenic (‑40/‑234) 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,7340AT 8.8% 298.7 / 15.2 113intergenic (‑40/‑234)panD/yadDaspartate 1‑decarboxylase/transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base A (50/52);  new base T (6/4);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 113 total aligned reads displayed.

TCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAACGC  >  NC_000913/146594‑146879
                                                                                                                                            |                                                                                                                                                 
tCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCCCGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCtt                                                                                                                                          <  2:87483/150‑1 (MQ=255)
 catgAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCATCGCCTGCttc                                                                                                                                         <  1:1569435/147‑2 (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:1844362/150‑5 (MQ=255)
    gAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCATCGCCTGCttcgtt                                                                                                                                      <  1:348678/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)
            tGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATtct                                                                                                                              <  2:1062068/150‑1 (MQ=255)
            tGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTAACTTTCTACCCTGTCGTTAACTACGCAGGCGATGCCTGCGTTGAGGAAATtca                                                                                                                              >  1:182895/1‑149 (MQ=255)
            tGGTCAATGGCGCAAGAACCTTCATAGTGCAGCTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTCCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGAGATGCCTGCTTTGAGGAAATtca                                                                                                                              >  1:915663/1‑149 (MQ=255)
                 aaTGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGc                                                                                                                         >  2:1604719/1‑150 (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)
                        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)
                              ccTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGAt                                                                                                            >  2:1294162/1‑150 (MQ=255)
                                   aTAGTGCAGGGCCGCAAGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCGTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTg                                                                                                       <  2:1125905/150‑1 (MQ=255)
                                      gTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGttt                                                                                                    <  2:1770881/150‑1 (MQ=255)
                                      gTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGttt                                                                                                    <  2:1270514/150‑1 (MQ=255)
                                          aGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGCTGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAgtg                                                                                                >  2:437239/1‑150 (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)
                                                         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)
                                                                  cgGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTGACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCg                                                                        >  1:329271/1‑150 (MQ=255)
                                                                  cgGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGGCGTTAACGAAGCAGGCGTTGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCg                                                                        <  2:97698/150‑1 (MQ=255)
                                                                    gTGGATGTTGCCCTGCAGCATCGTGCGTATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGAATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCACTTGCggc                                                                      <  1:725415/150‑3 (MQ=255)
                                                                    gTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTc                                                                      >  2:233368/1‑150 (MQ=255)
                                                                      ggAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGGCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTcag                                                                    <  2:551092/150‑1 (MQ=255)
                                                                        aGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTcagca                                                                  >  1:416576/1‑150 (MQ=255)
                                                                        aGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTATACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGTCATTTGCGTcagca                                                                  >  1:67859/1‑150 (MQ=255)
                                                                            tGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTAAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGttt                                                              <  1:645594/150‑1 (MQ=255)
                                                                              cccTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATACTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGc                                                            >  2:1634966/1‑150 (MQ=255)
                                                                                cTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGTCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGGCAGCAGTTTGCGt                                                          <  1:1010133/150‑1 (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                                         >  2:1466368/1‑150 (MQ=255)
                                                                                                 tCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGc                                         <  1:1285884/150‑1 (MQ=255)
                                                                                                   aTAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGTAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCa                                       >  2:90924/1‑150 (MQ=255)
                                                                                                   aTAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCTTACCGCGCCGAAGCGTGCCa                                       >  1:1664112/1‑150 (MQ=255)
                                                                                                   aTAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTTCCGCGCCGAAGCGTGCCa                                       >  1:985991/1‑150 (MQ=255)
                                                                                                    aaaCTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCAt                                      <  1:927274/149‑1 (MQ=255)
                                                                                                     aaCTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATg                                     <  1:1534225/150‑1 (MQ=255)
                                                                                                     aaCTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATg                                     <  1:1473503/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)
                                                                                                            aCCTTTCTACCCTGTCGTTAACGAAGCAGGCGATTCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGcc                              <  2:1114694/150‑1 (MQ=255)
                                                                                                               tttCTACCCTGTCGTTACCGCCGCCGGCGATGCCTGCTTTGAGGACTTCCGCCCGCTCTATAGCACGATTTTTGTTTAGTGTCTACTCATCTGAAGGCATTTGCGTCACCAGTTTGCGTACCGCGCCGAAGCGTGCCATGAAGCGCCTGt                           <  2:294190/150‑1 (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)
                                                                                                                      ccTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAACTTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCAACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAAt                    <  1:1663797/150‑1 (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)
                                                                                                                                      cAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGAGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAAc    <  1:247927/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
TCAAGAAAATCCTGGTCAATGGCGCAAGAACCTTCATAGTGCAGGTCCGCATGAGTCACTTTCACGCGGTGGAGTTTGCCCTGCAGCATCGTGCGAATCATAACTTCTACCTTTCTACCCTGTCGTTAACGAAGCAGGCGATGCCTGCTTTGAGGAAATTCTCACGCAGTATTGCCCGATTTTTGTTTAGTGTCTACTCATCTGACGGCATTTGCGTCAGCAGTTTGCGTACCGCGCCGAAGCGTGCCATGTATCGCCTGTAACGAATTGATATTTCTCTGAACGC  >  NC_000913/146594‑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: