Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 458,738 G→A 6.8% intergenic (+38/‑150) clpX → / → lon ATP‑dependent Clp protease ATP‑binding subunit ClpX/Lon protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913458,7380GA 6.8% 312.2 / 8.5 102intergenic (+38/‑150)clpX/lonATP‑dependent Clp protease ATP‑binding subunit ClpX/Lon protease
Reads supporting (aligned to +/- strand):  ref base G (41/54);  new base A (2/5);  total (43/59)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

Only 100 of 110 total aligned reads displayed.

GTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCT  >  NC_000913/458589‑458887
                                                                                                                                                     |                                                                                                                                                     
gTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAggggg                                                                                                                                                       >  1:844661/1‑150 (MQ=255)
 tCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAgggggg                                                                                                                                                      >  1:785693/1‑150 (MQ=255)
       gAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTa                                                                                                                                                <  2:673954/150‑1 (MQ=255)
       gAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTa                                                                                                                                                >  1:335266/1‑150 (MQ=255)
        aaGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGAtaaaat                                                                                                                                               >  2:1419186/1‑145 (MQ=255)
         aGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtc                                                                                                                                              <  1:42044/150‑1 (MQ=255)
           aCGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtctc                                                                                                                                            >  1:1088315/1‑150 (MQ=255)
              tCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTcccc                                                                                                                                         >  1:737837/1‑150 (MQ=255)
               cGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACCGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTAAACCAAAAAGGCGGGATTTTAtctaacct                                                                                                                                        >  2:1498117/1‑145 (MQ=255)
               cGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCt                                                                                                                                        <  1:1099504/150‑1 (MQ=255)
               cGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCt                                                                                                                                        >  2:405204/1‑150 (MQ=255)
               cGAAAAAGCGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCt                                                                                                                                        <  1:184013/150‑1 (MQ=255)
                gAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtt                                                                                                                                       <  2:432247/150‑1 (MQ=255)
                  aaaaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTa                                                                                                                                     >  1:1357628/1‑150 (MQ=255)
                   aaaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCttttt                                                                                                                                    <  2:1325729/150‑3 (MQ=255)
                   aaaGTGGTTATCGACGAGTCGGAAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTTTTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCttttt                                                                                                                                    <  2:592000/150‑3 (MQ=255)
                    aaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttg                                                                                                                                   <  1:1419186/150‑4 (MQ=255)
                    aaGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGCGGGATTTTATCTCACCTTTAAt                                                                                                                                   >  2:720259/1‑150 (MQ=255)
                         gTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTc                                                                                                                              <  2:892900/150‑1 (MQ=255)
                              cgacgaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTa                                                                                                                         >  2:391584/1‑150 (MQ=255)
                               gacgaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTAt                                                                                                                        <  1:1086234/150‑1 (MQ=255)
                               gacgaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTAt                                                                                                                        >  2:1030202/1‑150 (MQ=255)
                                 cgaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtc                                                                                                                      >  1:56335/1‑150 (MQ=255)
                                 cgaGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtc                                                                                                                      >  1:61891/1‑150 (MQ=255)
                                    gTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCg                                                                                                                   >  1:1013791/1‑150 (MQ=255)
                                    gTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCg                                                                                                                   <  1:391584/150‑1 (MQ=255)
                                       ggTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCg                                                                                                                >  1:1455734/1‑150 (MQ=255)
                                       ggTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCGCTATTCTCGGCg                                                                                                                >  2:480472/1‑150 (MQ=255)
                                       ggTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATAGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCg                                                                                                                >  2:601962/1‑150 (MQ=255)
                                          aaTTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTg                                                                                                             <  1:1058593/150‑1 (MQ=255)
                                           aTTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGa                                                                                                            <  2:1088315/150‑1 (MQ=255)
                                              gATGGTCAAAGCAAACCGTTGGTGATTTATGGCAAGCCGGTAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCATTGAAtg                                                                                                         <  2:1435095/150‑1 (MQ=255)
                                                tGGTCAAAGCAAAGCGTTGCTGATTTATGGCAAGCGGGAAGCGCAACAGGCATCTTGTGAATATTTAACCATTCCCAAACAATTAGTTAAGCAACAAGGGGGGATTTTATCTCCGCTTTTATTTTTCCTCTATTCTCGGCGTTGAAtgtg                                                                                                       <  2:860756/150‑1 (MQ=255)
                                                    cAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTggggg                                                                                                   >  2:1068198/1‑150 (MQ=255)
                                                       aGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGaaa                                                                                                <  2:1357628/150‑1 (MQ=255)
                                                       aGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGaaa                                                                                                >  1:1096057/1‑150 (MQ=255)
                                                        gCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAAc                                                                                               <  1:944058/150‑1 (MQ=255)
                                                        gCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAAc                                                                                               >  1:557146/1‑150 (MQ=255)
                                                                 tGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCata                                                                                      <  1:490199/150‑1 (MQ=255)
                                                                    tGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATAc                                                                                   >  1:430812/1‑150 (MQ=255)
                                                                    tGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCGCGGCGTTGAATGTGGGGGAAACATCCCCATATAc                                                                                   <  2:275260/150‑1 (MQ=255)
                                                                     gATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACt                                                                                  <  2:1455734/150‑1 (MQ=255)
                                                                         tATGGCACGCCCGAAGCGCAACAGGCATCTGGTTAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACg                                                                              <  2:1302784/150‑1 (MQ=255)
                                                                         tATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACg                                                                              <  2:737837/150‑1 (MQ=255)
                                                                         tATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACg                                                                              <  1:1030202/150‑1 (MQ=255)
                                                                           tGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTa                                                                            <  2:892625/150‑1 (MQ=255)
                                                                            ggCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtct                                                                                                                     <  2:537891/108‑1 (MQ=255)
                                                                            ggCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtct                                                                                                                     >  1:537891/1‑108 (MQ=255)
                                                                             gCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACa                                                                          <  2:230268/150‑1 (MQ=255)
                                                                              cAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACAt                                                                         >  2:646786/1‑150 (MQ=255)
                                                                                 gCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGtt                                                                      <  2:393611/150‑1 (MQ=255)
                                                                                 gCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGtt                                                                      <  1:916537/150‑1 (MQ=255)
                                                                                 gCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGtt                                                                      <  2:212292/150‑1 (MQ=255)
                                                                                 gCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGtt                                                                      <  1:740242/150‑1 (MQ=255)
                                                                                     gAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCGCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATa                                                                  >  2:972676/1‑150 (MQ=255)
                                                                                          gcAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATgg                                                             <  1:801597/150‑1 (MQ=255)
                                                                                            aaCAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCg                                                           <  1:1038097/150‑1 (MQ=255)
                                                                                             aCAGGCATCTGGTGAATAATTAACCATTCCAATACAATTAGTTAACAAAAAAGGGGGGATTTTATATCCCATTTAATTTTTCATCTATAAACTCATTTCAATATGGGGGAAACATACCCATATACTAACTAAAATTTTAAAAGAAAGCGt                                                          >  2:661434/1‑150 (MQ=255)
                                                                                               aGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                                                        <  2:844661/150‑1 (MQ=255)
                                                                                               aGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                                                        <  2:236343/150‑1 (MQ=255)
                                                                                                 gCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAg                                                      <  2:343634/150‑1 (MQ=255)
                                                                                                 gCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAg                                                      <  2:156519/150‑1 (MQ=255)
                                                                                                       ggTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGt                                                >  2:561116/1‑150 (MQ=255)
                                                                                                       ggTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGt                                                >  1:647717/1‑150 (MQ=255)
                                                                                                        gTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTc                                               >  1:1241024/1‑150 (MQ=255)
                                                                                                        gTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTc                                               >  1:1515848/1‑150 (MQ=255)
                                                                                                         tGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCg                                              >  1:339790/1‑150 (MQ=255)
                                                                                                            ataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtgt                                           >  1:1114391/1‑150 (MQ=255)
                                                                                                             taatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGGTGGCGTGAAGCACAGTCGTGTc                                          >  2:1365508/1‑150 (MQ=255)
                                                                                                             taatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTc                                          <  2:785693/150‑1 (MQ=255)
                                                                                                             taatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGACGTGAAGCACAGTCGTGTc                                          >  2:1173401/1‑150 (MQ=255)
                                                                                                             taatTAACCATTCCCATACAATTAGTTAACAAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTc                                          >  2:458113/1‑150 (MQ=255)
                                                                                                              aatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGTGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCa                                         >  1:1034845/1‑150 (MQ=255)
                                                                                                               atTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCACCTTTAATTTTTCCTCTATTGTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCAt                                        >  2:1249506/1‑150 (MQ=255)
                                                                                                                ttAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATc                                       <  2:61891/150‑1 (MQ=255)
                                                                                                                ttAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATc                                       <  2:56335/150‑1 (MQ=255)
                                                                                                                 tAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCt                                      <  2:755387/150‑1 (MQ=255)
                                                                                                                   aCCATTCCCATACAATTAGTTAACCAGAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGa                                    <  1:1536705/150‑1 (MQ=255)
                                                                                                                     cATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGAtt                                  >  1:755319/1‑150 (MQ=255)
                                                                                                                     cATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGAtt                                  <  2:1096057/150‑1 (MQ=255)
                                                                                                                      aTTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTa                                 <  2:335266/150‑1 (MQ=255)
                                                                                                                      aTTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTa                                 <  1:203956/150‑1 (MQ=255)
                                                                                                                        tCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAcc                               <  2:1013791/150‑1 (MQ=255)
                                                                                                                          ccATACAATTANTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTg                             >  1:1373760/1‑150 (MQ=255)
                                                                                                                          ccATACAATTAGTTAACCAAAAAGGCGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTGATCTGATTACCTg                             <  2:170199/150‑1 (MQ=255)
                                                                                                                             tACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCg                          <  1:369801/150‑1 (MQ=255)
                                                                                                                              aCAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCgg                         <  1:839397/150‑1 (MQ=255)
                                                                                                                               cAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGa                        >  2:750234/1‑150 (MQ=255)
                                                                                                                                 aTTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGaaa                      <  1:24136/150‑1 (MQ=255)
                                                                                                                                    aGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTa                   <  1:894285/150‑1 (MQ=255)
                                                                                                                                         aaattaaaGGGGAGATAAAATCCCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTa              >  2:303531/6‑150 (MQ=255)
                                                                                                                                           cAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAag            <  2:1373760/150‑1 (MQ=255)
                                                                                                                                            aaaaaGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAaga           <  1:561116/150‑1 (MQ=255)
                                                                                                                                             aaaaGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAagag          >  2:650381/1‑150 (MQ=255)
                                                                                                                                             aaaaGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAagag          <  1:567774/150‑1 (MQ=255)
                                                                                                                                             aaaaGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAagag          <  2:211855/150‑1 (MQ=255)
                                                                                                                                                  gggggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGc     <  1:750234/150‑1 (MQ=255)
                                                                                                                                                   ggggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGct    >  1:300156/1‑150 (MQ=255)
                                                                                                                                                    gggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGctc   <  1:405204/150‑1 (MQ=255)
                                                                                                                                                     ggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGctct  >  1:980931/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCT  >  NC_000913/458589‑458887

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: