Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,633,634 C→T 17.9% E38K (GAG→AAG)  yjjX ← non‑canonical purine NTP phosphatase, ITPase/XTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,633,6340CT17.9% 1264.4 / 231.5 502E38K (GAG→AAG) yjjXnon‑canonical purine NTP phosphatase, ITPase/XTPase
Reads supporting (aligned to +/- strand):  ref base C (199/212);  new base T (38/52);  total (238/264)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 524 total aligned reads displayed.

CCAACTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTAT  >  NC_000913/4633485‑4633783
                                                                                                                                                     |                                                                                                                                                     
ccAACTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCtctc                                                                                                                                                       >  2:2740159/1‑150 (MQ=255)
 cAACTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCg                                                                                                                                                      >  2:4072719/1‑150 (MQ=255)
 cAACTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCg                                                                                                                                                      >  2:726414/1‑150 (MQ=255)
  aaCTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGa                                                                                                                                                     <  1:1734037/150‑1 (MQ=255)
    cTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACg                                                                                                                                                   <  1:5490575/150‑1 (MQ=255)
       aaaGTGCTGTCACCATCGATGCCAGCTACAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATGTCGTGCGCCAGCGCGCGATTCCTCACTGCCAAAGGGCTGTTCCGGTCCACCGCTCTCGACGGCa                                                                                                                                                >  1:4645527/1‑150 (MQ=255)
           tGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACgg                                                                                                                                            >  2:1793710/1‑150 (MQ=255)
           tGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACgg                                                                                                                                            >  1:4017508/1‑150 (MQ=255)
           tGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACgg                                                                                                                                            <  1:3155687/150‑1 (MQ=255)
              tGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATg                                                                                                                                         >  1:5551334/1‑150 (MQ=255)
                tCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAATGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCa                                                                                                                                       >  1:158086/1‑150 (MQ=255)
                tCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCa                                                                                                                                       <  2:5072649/150‑1 (MQ=255)
                  aCCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAat                                                                                                                                     >  1:2862470/1‑150 (MQ=255)
                   ccATCGATGCCAGCTTCTATCGCCACCCAAAAATCAGCCTCTGGAAGTAACCGGAGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCTACGGATGCAata                                                                                                                                    <  2:5003292/150‑1 (MQ=255)
                        gATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCa                                                                                                                               <  1:3933266/150‑1 (MQ=255)
                           gCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGTAACGGATGCAATATGGTAGGa                                                                                                                            >  2:4712342/1‑150 (MQ=255)
                             cAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATc                                                                                                                          >  1:593051/1‑150 (MQ=255)
                             cAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATc                                                                                                                          >  1:5179080/1‑150 (MQ=255)
                             cAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATc                                                                                                                          >  1:4995942/1‑150 (MQ=255)
                              aGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATcc                                                                                                                         <  1:2662380/150‑1 (MQ=255)
                                cTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCtt                                                                                                                       >  2:5933271/1‑150 (MQ=255)
                                cTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCtt                                                                                                                       <  2:2311247/150‑1 (MQ=255)
                                   cAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGc                                                                                                                    >  2:1461995/1‑150 (MQ=255)
                                   cAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGc                                                                                                                    >  2:3095297/1‑150 (MQ=255)
                                    aaTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGcc                                                                                                                   <  2:246421/150‑1 (MQ=255)
                                       cGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGaa                                                                                                                <  1:2993355/150‑1 (MQ=255)
                                        gCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAg                                                                                                               >  2:583490/1‑150 (MQ=255)
                                           aCCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtc                                                                                                            >  1:1305937/1‑150 (MQ=255)
                                            cccAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGAACACCGCTCTCGACGGCAACGGATGCACTATGGCAGGATCCTTCGCCGACAAtct                                                                                                           >  2:3950450/1‑150 (MQ=255)
                                              cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCg                                                                                                         >  2:251836/1‑150 (MQ=255)
                                              cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCg                                                                                                         >  2:4468383/1‑150 (MQ=255)
                                              cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCg                                                                                                         <  2:6135096/150‑1 (MQ=255)
                                              cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCg                                                                                                         >  2:6214445/1‑150 (MQ=255)
                                               aaaaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGt                                                                                                        <  1:2564795/150‑1 (MQ=255)
                                               aaaaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGt                                                                                                        >  1:4880275/1‑150 (MQ=255)
                                               aaaaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGt                                                                                                        <  1:1006605/150‑1 (MQ=255)
                                                         ctctGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTg                                                                                              <  1:600867/150‑1 (MQ=255)
                                                          tctGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGc                                                                                             >  2:1195421/1‑150 (MQ=255)
                                                          tctGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGc                                                                                             >  1:360696/1‑150 (MQ=255)
                                                            tGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAg                                                                                           >  2:100321/1‑150 (MQ=255)
                                                                aGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATg                                                                                       <  2:2147270/150‑1 (MQ=255)
                                                                   aaaCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGcc                                                                                    >  2:5218489/1‑150 (MQ=255)
                                                                       ggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGaa                                                                                <  1:1408530/150‑1 (MQ=255)
                                                                       ggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGaa                                                                                <  2:4327272/150‑1 (MQ=255)
                                                                            gCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTa                                                                           >  2:3540113/1‑150 (MQ=255)
                                                                            gCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTa                                                                           <  2:1867186/150‑1 (MQ=255)
                                                                            gCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTa                                                                           >  1:5831288/1‑150 (MQ=255)
                                                                              aTTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGc                                                                         <  1:3748987/150‑1 (MQ=255)
                                                                               ttGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCg                                                                        <  1:4186313/150‑1 (MQ=255)
                                                                                     aCCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtg                                                                  <  1:4455345/150‑1 (MQ=255)
                                                                                     aCCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtg                                                                  >  2:1631061/1‑150 (MQ=255)
                                                                                      cccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtgg                                                                 >  1:2183346/1‑150 (MQ=255)
                                                                                         gATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggtgg                                                              >  1:2468632/1‑150 (MQ=255)
                                                                                           tttCGTGCGCCAGCGCGCGTTTCCTCACTGCCACAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTc                                                            <  2:108347/150‑1 (MQ=255)
                                                                                           tttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTc                                                            <  2:6372130/150‑1 (MQ=255)
                                                                                           tttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTc                                                            <  2:3793054/150‑1 (MQ=255)
                                                                                           tttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTc                                                            <  2:6486915/150‑1 (MQ=255)
                                                                                           tttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTc                                                            <  2:5436933/150‑1 (MQ=255)
                                                                                           tgtCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAAGGGCCTGAATTTTAGCGGGATTGGTGGTc                                                            <  1:2008437/148‑1 (MQ=255)
                                                                                            ttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCg                                                           <  2:1991730/150‑1 (MQ=255)
                                                                                                tGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcaca                                                       >  1:1264073/1‑150 (MQ=255)
                                                                                                   gcCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGAc                                                    <  2:4789932/150‑1 (MQ=255)
                                                                                                    ccAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACa                                                   >  2:2783797/1‑150 (MQ=255)
                                                                                                        cgcgcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACtt                                               <  1:3069677/150‑1 (MQ=255)
                                                                                                        cgcgcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACtt                                               <  2:1550281/150‑1 (MQ=255)
                                                                                                        cgcgcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACtt                                               <  2:3316180/150‑1 (MQ=255)
                                                                                                             gTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGc                                          <  2:2939716/150‑1 (MQ=255)
                                                                                                                tCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATa                                       >  1:3704108/1‑150 (MQ=255)
                                                                                                                  cTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAAt                                     <  1:4920482/150‑1 (MQ=255)
                                                                                                                  cTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAAt                                     <  1:5419647/150‑1 (MQ=255)
                                                                                                                   tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                    >  2:3406129/1‑150 (MQ=255)
                                                                                                                   tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATTCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                    <  1:1873151/150‑1 (MQ=255)
                                                                                                                   tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                    >  1:3434477/1‑150 (MQ=255)
                                                                                                                   tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                    >  2:3489448/1‑150 (MQ=255)
                                                                                                                     aCTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAg                                  <  1:233298/150‑1 (MQ=255)
                                                                                                                     aCTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAg                                  >  2:3213054/1‑150 (MQ=255)
                                                                                                                      cTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGc                                 <  2:3721769/150‑1 (MQ=255)
                                                                                                                       tGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCa                                >  1:4792147/1‑150 (MQ=255)
                                                                                                                        gCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCAt                               >  1:6390223/1‑150 (MQ=255)
                                                                                                                          cAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTa                             >  1:5325432/1‑150 (MQ=255)
                                                                                                                          cAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTa                             >  2:1403795/1‑150 (MQ=255)
                                                                                                                            aaGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACt                           <  2:6580713/150‑1 (MQ=255)
                                                                                                                              gggCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCa                         <  1:1045712/150‑1 (MQ=255)
                                                                                                                              gggCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCa                         <  1:5638763/150‑1 (MQ=255)
                                                                                                                              gggCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCa                         <  1:4321353/150‑1 (MQ=255)
                                                                                                                                gCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGa                       >  2:1286003/1‑150 (MQ=255)
                                                                                                                                gCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGa                       <  1:2780683/150‑1 (MQ=255)
                                                                                                                                 cTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGaa                      >  2:6263426/1‑150 (MQ=255)
                                                                                                                                    ttCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAAt                   >  1:6583346/1‑150 (MQ=255)
                                                                                                                                         gTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACg              >  1:874205/1‑150 (MQ=255)
                                                                                                                                            cacCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTa           <  2:3607762/150‑1 (MQ=255)
                                                                                                                                              ccGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAca         <  1:3355/150‑1 (MQ=255)
                                                                                                                                              ccGCTCTCGACGGCAACGGATTCATTTTGGTAGGATCATTCGCCGAAGATCTCGTGATATGCCTGCAGAATGGCCTGATTTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGTAAATTATCGTTAca         <  1:457155/150‑1 (MQ=255)
                                                                                                                                              ccGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAca         <  1:6214445/150‑1 (MQ=255)
                                                                                                                                              ccGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAca         <  1:1443951/150‑1 (MQ=255)
                                                                                                                                               cGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcag        >  1:240360/1‑150 (MQ=255)
                                                                                                                                                 ctctCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcagca      <  1:5719009/150‑1 (MQ=255)
                                                                                                                                                 ctctCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcagca      <  2:5778248/150‑1 (MQ=255)
                                                                                                                                                   ctcGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGt    <  2:5626950/150‑1 (MQ=255)
                                                                                                                                                   ctcGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGt    <  1:1174234/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAACTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTAT  >  NC_000913/4633485‑4633783

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