Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,633,657 C→T 100% C30Y (TGC→TAC)  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,6570CT96.9% 2911.9 / 44.1 819C30Y (TGC→TAC) yjjXnon‑canonical purine NTP phosphatase, ITPase/XTPase
Reads supporting (aligned to +/- strand):  ref base C (14/12);  new base T (381/410);  total (396/423)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 820 total aligned reads displayed.

CGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTAC  >  NC_000913/4633508‑4633803
                                                                                                                                                     |                                                                                                                                                  
  aTGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCACCGGATGCAATATGGTAg                                                                                                                                                  >  1:8468013/1‑150 (MQ=255)
    gCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGa                                                                                                                                                >  2:4062521/1‑150 (MQ=255)
      cAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATc                                                                                                                                              >  1:7342618/1‑150 (MQ=255)
       aGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATcc                                                                                                                                             >  2:4331278/1‑150 (MQ=255)
       aGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATcc                                                                                                                                             >  1:327231/1‑150 (MQ=255)
       aGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATcc                                                                                                                                             >  1:11175789/1‑150 (MQ=255)
          ttCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTc                                                                                                                                          >  2:10820697/1‑150 (MQ=255)
            cAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGc                                                                                                                                        <  1:9154272/150‑1 (MQ=255)
             aaTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGcc                                                                                                                                       <  1:8192347/150‑1 (MQ=255)
             aaTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGcc                                                                                                                                       >  1:6746658/1‑150 (MQ=255)
                  ccACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGa                                                                                                                                  <  2:4596424/150‑1 (MQ=255)
                  ccACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGa                                                                                                                                  >  2:8059343/1‑150 (MQ=255)
                      ccAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtctc                                                                                                                              <  2:5157891/150‑1 (MQ=255)
                      ccAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtctc                                                                                                                              >  2:3722755/1‑150 (MQ=255)
                      ccAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtctc                                                                                                                              <  2:1703694/150‑1 (MQ=255)
                      ccAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtctc                                                                                                                              <  1:2288007/150‑1 (MQ=255)
                      ccAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTACCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtctc                                                                                                                              >  2:551128/1‑150 (MQ=255)
                       cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCg                                                                                                                             <  1:4085899/150‑1 (MQ=255)
                       cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCg                                                                                                                             <  1:5374001/150‑1 (MQ=255)
                        aaaaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGt                                                                                                                            >  1:2690955/1‑150 (MQ=255)
                          aaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGa                                                                                                                          <  2:11586691/150‑1 (MQ=255)
                           aaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGaa                                                                                                                         <  1:686302/150‑1 (MQ=255)
                           aaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGaa                                                                                                                         <  2:981398/150‑1 (MQ=255)
                                 ccTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCt                                                                                                                   >  2:5765465/1‑150 (MQ=255)
                                       gAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGaa                                                                                                             >  1:8301535/1‑150 (MQ=255)
                                        aaGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAAt                                                                                                            <  2:8178543/150‑1 (MQ=255)
                                          gTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATgg                                                                                                          >  1:10262648/1‑150 (MQ=255)
                                           tAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGc                                                                                                         <  1:11918114/150‑1 (MQ=255)
                                           tAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGc                                                                                                         <  2:301935/150‑1 (MQ=255)
                                             aaCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCt                                                                                                       >  2:8742494/1‑150 (MQ=255)
                                             aaCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCt                                                                                                       <  1:10407059/150‑1 (MQ=255)
                                             aaCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCt                                                                                                       >  2:1415941/1‑150 (MQ=255)
                                               cggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGa                                                                                                     >  1:1784364/1‑150 (MQ=255)
                                               cggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGa                                                                                                     >  2:12286008/1‑150 (MQ=255)
                                                ggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGaa                                                                                                    <  1:10105786/150‑1 (MQ=255)
                                                  cggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAAtt                                                                                                  <  1:3734599/150‑1 (MQ=255)
                                                     gCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTa                                                                                               <  1:11278640/150‑1 (MQ=255)
                                                     gCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTa                                                                                               >  1:2571117/1‑150 (MQ=255)
                                                       aTTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGc                                                                                             >  1:5149560/1‑150 (MQ=255)
                                                       aTTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGc                                                                                             <  1:7907020/150‑1 (MQ=255)
                                                         tGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCgg                                                                                           <  2:8301535/150‑1 (MQ=255)
                                                           gCTACCCGCTTTCGTGCGCCTGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCTTTATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGACTGCAGAATGGCCTGAATTTTAGCTGGa                                                                                         <  2:11047233/150‑1 (MQ=255)
                                                             tACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGAtt                                                                                       >  1:9983449/1‑150 (MQ=255)
                                                               cccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtgg                                                                                     >  1:8714091/1‑150 (MQ=255)
                                                               cccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtgg                                                                                     <  2:10225861/150‑1 (MQ=255)
                                                                ccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggt                                                                                    >  2:733080/1‑150 (MQ=255)
                                                                ccGATTTCATGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggt                                                                                    >  1:11652507/1‑150 (MQ=255)
                                                                 cGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggtg                                                                                   <  2:9371666/150‑1 (MQ=255)
                                                                   aTTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATGCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTggtggt                                                                                 >  1:4539380/1‑150 (MQ=255)
                                                                   aTTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTggtggt                                                                                 >  1:4345306/1‑150 (MQ=255)
                                                                       cGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGca                                                                             >  2:6576513/1‑150 (MQ=255)
                                                                        gTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcac                                                                            <  1:9560615/150‑1 (MQ=255)
                                                                        gTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcac                                                                            <  1:3849321/150‑1 (MQ=255)
                                                                        gTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcac                                                                            <  1:1981501/150‑1 (MQ=255)
                                                                         tGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcaca                                                                           <  1:6319076/150‑1 (MQ=255)
                                                                         tGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcaca                                                                           >  2:2636317/1‑150 (MQ=255)
                                                                           cgcCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                                                         >  2:281405/1‑150 (MQ=255)
                                                                           cgcCAGCGCGAGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                                                         <  1:5823125/150‑1 (MQ=255)
                                                                            gcCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGAc                                                                        <  2:7072439/150‑1 (MQ=255)
                                                                            gcCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGAc                                                                        <  1:10750416/150‑1 (MQ=255)
                                                                              cAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACaa                                                                      >  1:10600613/1‑150 (MQ=255)
                                                                                  gcgcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTg                                                                  <  1:6648663/150‑1 (MQ=255)
                                                                                   cgcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTgg                                                                 <  1:4622794/150‑1 (MQ=255)
                                                                                    gcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGt                                                                <  2:11794884/150‑1 (MQ=255)
                                                                                     cgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTg                                                               <  2:6662868/150‑1 (MQ=255)
                                                                                         tCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATa                                                           <  2:592845/150‑1 (MQ=255)
                                                                                         tCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATa                                                           >  2:4851942/1‑150 (MQ=255)
                                                                                         tCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATa                                                           >  2:8482307/1‑150 (MQ=255)
                                                                                          ccTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATaa                                                          <  1:11363641/150‑1 (MQ=255)
                                                                                            tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                                        >  2:3409564/1‑150 (MQ=255)
                                                                                            tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                                        >  1:7675283/1‑150 (MQ=255)
                                                                                              aCTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAg                                                      >  2:10338720/1‑150 (MQ=255)
                                                                                                  ccAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCAtt                                                  >  2:9604338/1‑150 (MQ=255)
                                                                                                         gCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGa                                           >  2:6121073/1‑150 (MQ=255)
                                                                                                            gTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGaaaa                                        <  2:7689088/150‑1 (MQ=255)
                                                                                                             ttCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAGCTTGGTGCATAATCAGCATTACTCAGAAAAt                                       >  2:716626/1‑150 (MQ=255)
                                                                                                               ccGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTa                                     <  1:12091195/150‑1 (MQ=255)
                                                                                                                cGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTaa                                    <  2:2158743/150‑1 (MQ=255)
                                                                                                                 ggTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAAc                                   <  1:1626179/150‑1 (MQ=255)
                                                                                                                 ggTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAAc                                   >  2:6375758/1‑150 (MQ=255)
                                                                                                                 ggTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAAc                                   <  1:5005781/150‑1 (MQ=255)
                                                                                                                   tACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGt                                 >  1:3495934/1‑150 (MQ=255)
                                                                                                                       ccGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAca                             <  2:10516030/150‑1 (MQ=255)
                                                                                                                        cGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcag                            >  2:11199916/1‑150 (MQ=255)
                                                                                                                        cGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcag                            >  1:9030405/1‑150 (MQ=255)
                                                                                                                        cGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcag                            >  2:1966738/1‑150 (MQ=255)
                                                                                                                          ctctCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcagca                          <  1:9849561/150‑1 (MQ=255)
                                                                                                                             tcGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGta                       >  2:161447/1‑150 (MQ=255)
                                                                                                                               gACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGtata                     <  1:10356076/150‑1 (MQ=255)
                                                                                                                                aCGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAc                    <  2:4994582/150‑1 (MQ=255)
                                                                                                                                 cGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACg                   <  1:2162721/150‑1 (MQ=255)
                                                                                                                                        ggATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaa            <  2:3915928/150‑1 (MQ=255)
                                                                                                                                        ggATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaa            >  1:2840329/1‑150 (MQ=255)
                                                                                                                                        ggATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGaaaaaa            >  1:11440165/1‑150 (MQ=255)
                                                                                                                                        ggATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGaaaaaa            >  2:904167/1‑150 (MQ=255)
                                                                                                                                         gATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaaa           <  1:12286008/150‑1 (MQ=255)
                                                                                                                                             cAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCAt       <  2:12074348/150‑1 (MQ=255)
                                                                                                                                              aaTATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATg      >  2:6766517/1‑150 (MQ=255)
                                                                                                                                              aaTATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATg      >  1:302158/1‑150 (MQ=255)
                                                                                                                                                 atGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTa   <  1:7071640/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTAC  >  NC_000913/4633508‑4633803

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: