Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 662,686 A→T 6.2% intergenic (‑44/+66) ybeD ← / ← dacA DUF493 domain‑containing protein YbeD/D‑alanyl‑D‑alanine carboxypeptidase DacA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913662,6860AT 6.2% 324.2 / 6.1 103intergenic (‑44/+66)ybeD/dacADUF493 domain‑containing protein YbeD/D‑alanyl‑D‑alanine carboxypeptidase DacA
Reads supporting (aligned to +/- strand):  ref base A (48/48);  new base T (3/3);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

TACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGCAACAC  >  NC_000913/662538‑662835
                                                                                                                                                    |                                                                                                                                                     
tACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAg                                                                                                                                                      >  1:997580/1‑150 (MQ=255)
  ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGc                                                                                                                                                    >  1:1483318/1‑150 (MQ=255)
       tCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCTTCTGGCGGGAg                                                                                                                                               <  2:67824/150‑1 (MQ=255)
       tCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                                                                                                                               <  1:1344259/150‑1 (MQ=255)
       tCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                                                                                                                               <  1:486014/150‑1 (MQ=255)
       tCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                                                                                                                               >  1:87911/1‑150 (MQ=255)
        cAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTAGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGt                                                                                                                                              >  2:729231/1‑150 (MQ=255)
        cAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAACTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCTTCTGGCGGGAGt                                                                                                                                              <  1:1048046/150‑1 (MQ=255)
             accTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttttt                                                                                                                                         >  1:77420/1‑150 (MQ=255)
              cctcaccAACCACCTCACGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTCTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGTAGTTTTTTa                                                                                                                                        <  2:289531/144‑1 (MQ=255)
               cTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCTTCTGGCGGGAGttgctatt                                                                                                                                       <  1:1065399/150‑7 (MQ=255)
                 gATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTg                                                                                                                                     >  1:594961/1‑150 (MQ=255)
                 gATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTg                                                                                                                                     >  2:341402/1‑150 (MQ=255)
                  aTCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGa                                                                                                                                    <  1:86334/150‑1 (MQ=255)
                   tCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                                                                                                                                   >  1:1312886/1‑150 (MQ=255)
                   tCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGGTATTTACTGAt                                                                                                                                   >  2:1323217/1‑150 (MQ=255)
                     aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGc                                                                                                                                 >  1:1399712/1‑150 (MQ=255)
                      aCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCt                                                                                                                                >  1:205855/1‑150 (MQ=255)
                         aGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAg                                                                                                                             <  1:1535081/150‑1 (MQ=255)
                          gCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGt                                                                                                                            <  2:1046994/150‑1 (MQ=255)
                           cTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                                                                                                                           <  1:1299914/150‑1 (MQ=255)
                            tCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGtat                                                                                                                          <  1:776436/150‑1 (MQ=255)
                              aGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGt                                                                                                                            <  2:29008/146‑1 (MQ=255)
                              aGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGt                                                                                                                            >  1:29008/1‑146 (MQ=255)
                               ggTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATg                                                                                                                       >  2:242695/1‑150 (MQ=255)
                                  aaCGCCTGCCCTATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgggg                                                                                                                    <  2:1088771/150‑1 (MQ=255)
                                    cGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc                                                                                                                  >  1:1276200/1‑150 (MQ=255)
                                      ccTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACgg                                                                                                                >  1:1098517/1‑150 (MQ=255)
                                          ccccAAAACTTTGTAAGTAAAAGCAGTAGGGAATCCAAGCAGTTCGTTAAGTTTGCGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCGCGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAAt                                                                                                            <  2:349380/150‑1 (MQ=255)
                                             cttAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTac                                                                                                         <  1:454084/148‑1 (MQ=255)
                                                  cTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACttt                                                                                                    >  2:208145/1‑150 (MQ=255)
                                                   tttGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTc                                                                                                   <  2:1483318/150‑1 (MQ=255)
                                                    ttGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCa                                                                                                  >  2:1414759/1‑150 (MQ=255)
                                                      gtaagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAg                                                                                                >  2:302434/1‑150 (MQ=255)
                                                        aagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTATTATATGGGGACGGAAATTACACTTTCTAgtg                                                                                              <  2:3064/150‑1 (MQ=255)
                                                        aagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtg                                                                                              <  2:1276200/150‑1 (MQ=255)
                                                        aagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtg                                                                                              >  2:1317215/1‑150 (MQ=255)
                                                        aagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtg                                                                                              <  2:77420/150‑1 (MQ=255)
                                                        aagtaaAAGGAGTAGGGAATTCAAGCAGATCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTAATGCTGATTAGTATATGGGGACGGAAATAACACTTTCAAgtg                                                                                              >  2:1124237/1‑150 (MQ=255)
                                                         agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGATATTACACTTTCAAgtgt                                                                                             >  2:642377/1‑150 (MQ=255)
                                                             aaaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg                                                                                                                      <  1:530531/121‑1 (MQ=255)
                                                             aaaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg                                                                                                                      >  2:530531/1‑121 (MQ=255)
                                                              aaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCTGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAt                                                                                        <  1:673343/150‑1 (MQ=255)
                                                              aaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAt                                                                                        <  2:183076/150‑1 (MQ=255)
                                                                ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACGCTTTCAAGTGTTTAAttt                                                                                      >  1:845680/1‑150 (MQ=255)
                                                                ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAttt                                                                                      >  2:1129589/1‑150 (MQ=255)
                                                                     aGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaac                                                                                 <  2:235520/150‑1 (MQ=255)
                                                                      gggAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaacc                                                                                >  1:1188640/1‑150 (MQ=255)
                                                                        gAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaa                                                                              >  2:1241368/1‑150 (MQ=255)
                                                                           ttCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCATGTGTTTAATTTTTaaccaaacc                                                                           <  2:609386/150‑1 (MQ=255)
                                                                           ttCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaaacc                                                                           <  2:205855/150‑1 (MQ=255)
                                                                            tCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaaacca                                                                          >  2:863269/1‑150 (MQ=255)
                                                                             cAAGCAGTTCGTTAAGTTTGGTTTTCAAGTCAGCTCCGGAGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCCCGAGATTTTTACTGATGCTTAGAATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAg                                                                         >  2:648938/1‑150 (MQ=255)
                                                                              aaGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGt                                                                        >  1:1033366/1‑150 (MQ=255)
                                                                                 cAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAg                                                                                                <  1:383732/123‑1 (MQ=255)
                                                                                 cAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAg                                                                                                >  2:383732/1‑123 (MQ=255)
                                                                                  aGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATg                                                                    <  2:1432330/150‑1 (MQ=255)
                                                                                      ccttAAGTTTGGATATCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAAc                                                                >  2:471147/3‑150 (MQ=255)
                                                                                        ttAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAt                                                              <  2:171364/150‑1 (MQ=255)
                                                                                              ttGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaattt                                                        <  2:1312886/150‑1 (MQ=255)
                                                                                              ttGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaattt                                                        <  2:952517/150‑1 (MQ=255)
                                                                                               tGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaattta                                                       >  2:443258/1‑150 (MQ=255)
                                                                                                ggTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatttaa                                                      >  2:1057700/1‑150 (MQ=255)
                                                                                                ggTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatttaa                                                      <  1:30524/150‑1 (MQ=255)
                                                                                                ggTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatttaa                                                      >  2:1243521/1‑150 (MQ=255)
                                                                                                  ttttCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATg                                                    <  2:1457770/150‑1 (MQ=255)
                                                                                                  ttttCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATg                                                    <  2:895070/150‑1 (MQ=255)
                                                                                                   tttCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGt                                                   <  2:1399712/150‑1 (MQ=255)
                                                                                                      cATGACAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGATTAGTATATGGGGACGGAAATTACAATTTAAAGTGTTTAATTTTTAACAAAACCAGTGATGGAAAATTAATTTAATGTaaa                                                >  1:669917/1‑149 (MQ=255)
                                                                                                       aTGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATc                                               <  2:817990/150‑1 (MQ=255)
                                                                                                       aTGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATc                                               <  2:969314/150‑1 (MQ=255)
                                                                                                        tGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCa                                              <  1:1183492/150‑1 (MQ=255)
                                                                                                         gTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAAGATTAATTTAATGTAATCaa                                             >  1:1134512/1‑150 (MQ=255)
                                                                                                           cAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATg                                           <  2:475548/150‑1 (MQ=255)
                                                                                                                 cGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACGCTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTg                                     <  2:845680/150‑1 (MQ=255)
                                                                                                                  ggCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCTAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGc                                    <  1:1414759/150‑1 (MQ=255)
                                                                                                                  ggCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGc                                    >  1:84562/1‑150 (MQ=255)
                                                                                                                  gccgtaacgtaaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGc                                    >  1:442770/3‑150 (MQ=255)
                                                                                                                   gCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGcc                                   >  1:232670/1‑150 (MQ=255)
                                                                                                                     gtaacgtaaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATTCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCga                                 <  2:1409084/150‑1 (MQ=255)
                                                                                                                       aacgtaaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCgaag                               >  1:1085554/1‑150 (MQ=255)
                                                                                                                           taaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGt                           >  1:466068/1‑150 (MQ=255)
                                                                                                                             aTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTa                         >  2:1328117/1‑150 (MQ=255)
                                                                                                                                  taaaaaaCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTc                    >  2:698227/6‑150 (MQ=255)
                                                                                                                                      cAACTCCCGCCATAATGCGGGATTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCggg                <  2:669917/150‑1 (MQ=255)
                                                                                                                                      cAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCggg                <  1:242695/150‑1 (MQ=255)
                                                                                                                                       aaCTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGa               >  1:1190673/1‑150 (MQ=255)
                                                                                                                                        aCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATGTGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGAt              >  2:691745/1‑150 (MQ=255)
                                                                                                                                        aCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGAt              >  1:1497148/1‑150 (MQ=255)
                                                                                                                                         cTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGAtt             >  1:1263636/1‑150 (MQ=255)
                                                                                                                                         cTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGAtt             <  1:1243521/150‑1 (MQ=255)
                                                                                                                                           cccGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTc           >  1:1108523/1‑150 (MQ=255)
                                                                                                                                             cGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCGGAAGAAGTTACCTTCCGGGATTTCtt         >  1:1280810/1‑150 (MQ=255)
                                                                                                                                             cGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCtt         >  1:580351/1‑150 (MQ=255)
                                                                                                                                               ccAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGc       <  1:943571/150‑1 (MQ=255)
                                                                                                                                               ccAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGc       <  1:1020787/150‑1 (MQ=255)
                                                                                                                                               ccAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGc       <  1:302434/150‑1 (MQ=255)
                                                                                                                                                  gAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGCAac    <  1:443258/150‑1 (MQ=255)
                                                                                                                                                    aGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGCAacac  <  2:883688/150‑1 (MQ=255)
                                                                                                                                                    aGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGCAacac  <  2:585980/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTCTTGCAACAC  >  NC_000913/662538‑662835

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