Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,633,657 C→T 5.1% 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,6570CT 5.1% 1521.2 / 40.1 434C30Y (TGC→TAC) yjjXnon‑canonical purine NTP phosphatase, ITPase/XTPase
Reads supporting (aligned to +/- strand):  ref base C (194/217);  new base T (10/12);  total (205/229)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 435 total aligned reads displayed.

GATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGG  >  NC_000913/4633509‑4633806
                                                                                                                                                    |                                                                                                                                                     
gATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCa                                                                                                                                                      <  2:1130501/150‑1 (MQ=255)
  tGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAgg                                                                                                                                                    >  1:4171602/1‑150 (MQ=255)
    ccAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGAt                                                                                                                                                  <  2:823267/150‑1 (MQ=255)
    ccAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGAt                                                                                                                                                  >  1:2722872/1‑150 (MQ=255)
      aGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATcc                                                                                                                                                >  2:5385331/1‑150 (MQ=255)
       gCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCt                                                                                                                                               >  1:4859411/1‑150 (MQ=255)
        cTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCtt                                                                                                                                              >  2:3301006/1‑150 (MQ=255)
           cAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAAGGGATGCAATATGGCAGGATCCTTCGc                                                                                                                                           <  2:5244553/150‑1 (MQ=255)
           cAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGc                                                                                                                                           <  2:1432894/150‑1 (MQ=255)
            aaTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGcc                                                                                                                                          <  1:3384946/150‑1 (MQ=255)
             aTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCg                                                                                                                                         >  2:5151810/1‑150 (MQ=255)
             aTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCg                                                                                                                                         <  1:3646553/150‑1 (MQ=255)
             aTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCg                                                                                                                                         >  2:3131763/1‑150 (MQ=255)
             aTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCg                                                                                                                                         <  1:708787/150‑1 (MQ=255)
              tCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGa                                                                                                                                        >  1:358896/1‑150 (MQ=255)
              tCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGa                                                                                                                                        >  1:5251143/1‑150 (MQ=255)
                     ccAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGAtctc                                                                                                                                 <  2:2092467/150‑1 (MQ=255)
                     ccAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGAtctc                                                                                                                                 <  1:2531281/150‑1 (MQ=255)
                     ccAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGAtctc                                                                                                                                 <  2:3885698/150‑1 (MQ=255)
                      cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCg                                                                                                                                >  2:4214035/1‑150 (MQ=255)
                      cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCg                                                                                                                                <  1:3803008/150‑1 (MQ=255)
                      cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCg                                                                                                                                >  1:4515531/1‑150 (MQ=255)
                       aaaaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGt                                                                                                                               >  2:5078111/1‑150 (MQ=255)
                           ctCAGCCTCTGGTAGTAACCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTGTCCTCACTGCCAACGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGaaa                                                                                                                           <  2:2991778/149‑1 (MQ=255)
                               gCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGcc                                                                                                                       >  2:4547017/1‑150 (MQ=255)
                                ccTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCt                                                                                                                      >  1:3544815/1‑150 (MQ=255)
                                    tGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAg                                                                                                                  <  2:2933844/150‑1 (MQ=255)
                                     ggAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGa                                                                                                                 >  2:5815300/1‑150 (MQ=255)
                                        aGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATg                                                                                                              >  1:2421870/1‑150 (MQ=255)
                                         gTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATgg                                                                                                             <  1:1108798/150‑1 (MQ=255)
                                          tAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGc                                                                                                            >  1:928474/1‑150 (MQ=255)
                                          tAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGc                                                                                                            >  1:534053/1‑150 (MQ=255)
                                          tAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGc                                                                                                            <  2:649698/150‑1 (MQ=255)
                                            aaCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCt                                                                                                          >  2:4063314/1‑150 (MQ=255)
                                            aaCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCt                                                                                                          >  1:4569207/1‑150 (MQ=255)
                                              cggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGa                                                                                                        >  1:1280795/1‑150 (MQ=255)
                                              cggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGa                                                                                                        <  2:1700717/150‑1 (MQ=255)
                                                      aTTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGc                                                                                                >  1:3941176/1‑150 (MQ=255)
                                                        tGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCgg                                                                                              <  2:4930956/150‑1 (MQ=255)
                                                         ggCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCggg                                                                                             <  1:5507450/150‑1 (MQ=255)
                                                           ntACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCACTATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGACTTTTAGCGGGAt                                                                                           >  1:2038458/2‑150 (MQ=255)
                                                              cccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtgg                                                                                        <  2:3755099/150‑1 (MQ=255)
                                                               ccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggt                                                                                       >  2:3540722/1‑150 (MQ=255)
                                                               ccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggt                                                                                       <  1:3307702/150‑1 (MQ=255)
                                                               ccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggt                                                                                       >  2:3646110/1‑150 (MQ=255)
                                                               ccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggt                                                                                       >  2:2517884/1‑150 (MQ=255)
                                                               ccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggt                                                                                       <  1:5193426/150‑1 (MQ=255)
                                                                 gATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggtgg                                                                                     <  2:1529713/150‑1 (MQ=255)
                                                                       gTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcac                                                                               >  2:1208852/1‑150 (MQ=255)
                                                                          cgcCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                                                            <  2:4988092/150‑1 (MQ=255)
                                                                          cgcCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                                                            <  2:5738968/150‑1 (MQ=255)
                                                                          cgcCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                                                            <  2:4859411/150‑1 (MQ=255)
                                                                              aGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAAc                                                                        >  2:2268186/1‑150 (MQ=255)
                                                                                   gcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGt                                                                   <  1:3190962/150‑1 (MQ=255)
                                                                                   gcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGt                                                                   >  1:5058641/1‑150 (MQ=255)
                                                                                       ttCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCAt                                                               <  2:609337/150‑1 (MQ=255)
                                                                                           tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                                           >  2:1058932/1‑150 (MQ=255)
                                                                                           tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                                           >  2:1913271/1‑150 (MQ=255)
                                                                                           tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                                           >  1:1737865/1‑150 (MQ=255)
                                                                                            cACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCa                                                          >  1:5538900/1‑150 (MQ=255)
                                                                                            cACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCa                                                          <  2:355026/150‑1 (MQ=255)
                                                                                                 ccAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCAtt                                                     >  1:1571477/1‑150 (MQ=255)
                                                                                                  cAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTa                                                    <  1:1182530/150‑1 (MQ=255)
                                                                                                  cAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTa                                                    >  2:4749957/1‑150 (MQ=255)
                                                                                                  cAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTa                                                    <  1:577983/150‑1 (MQ=255)
                                                                                                    aaGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACt                                                  >  1:10031/1‑150 (MQ=255)
                                                                                                    aaGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACt                                                  <  1:1419052/150‑1 (MQ=255)
                                                                                                            ttCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAAt                                          <  2:3209044/150‑1 (MQ=255)
                                                                                                            ttCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAAt                                          >  2:3415427/1‑150 (MQ=255)
                                                                                                             tCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAAtt                                         >  1:2980016/1‑150 (MQ=255)
                                                                                                                  tACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGt                                    <  1:4973438/150‑1 (MQ=255)
                                                                                                                     acCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAc                                 >  2:3777794/1‑150 (MQ=255)
                                                                                                                      ccGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAca                                <  1:3009407/150‑1 (MQ=255)
                                                                                                                       cGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcag                               >  2:2582970/1‑150 (MQ=255)
                                                                                                                        gCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcagc                              <  1:468765/150‑1 (MQ=255)
                                                                                                                        gCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcagc                              <  1:1933184/150‑1 (MQ=255)
                                                                                                                         ctctCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcagca                             >  2:545062/1‑150 (MQ=255)
                                                                                                                              gACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGtata                        >  2:5630552/1‑150 (MQ=255)
                                                                                                                               aCGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAc                       <  1:3249941/150‑1 (MQ=255)
                                                                                                                                cGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACg                      <  2:4463933/150‑1 (MQ=255)
                                                                                                                                cGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACg                      >  1:4059844/1‑150 (MQ=255)
                                                                                                                                 ggCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACgg                     <  1:1652093/150‑1 (MQ=255)
                                                                                                                                 ggCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACgg                     >  2:2857981/1‑150 (MQ=255)
                                                                                                                                 ggCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACgg                     <  2:5576439/150‑1 (MQ=255)
                                                                                                                                   cAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaa                   >  1:5714789/1‑150 (MQ=255)
                                                                                                                                   cAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaa                   <  1:3472573/150‑1 (MQ=255)
                                                                                                                                   cAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaa                   >  1:4793993/1‑150 (MQ=255)
                                                                                                                                       ggATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaa               >  2:4168793/1‑150 (MQ=255)
                                                                                                                                       ggATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaa               >  1:2164572/1‑150 (MQ=255)
                                                                                                                                        gATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaaa              <  2:3936580/150‑1 (MQ=255)
                                                                                                                                        gATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaaa              <  1:544314/150‑1 (MQ=255)
                                                                                                                                         aTGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAGAATTAACGTTACAGCAGTATACGGAAAAAAAg             <  1:4204625/150‑1 (MQ=255)
                                                                                                                                         aTGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAg             <  1:3742144/150‑1 (MQ=255)
                                                                                                                                          tGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGc            <  1:2833114/150‑1 (MQ=255)
                                                                                                                                          tGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGc            <  2:260744/150‑1 (MQ=255)
                                                                                                                                          tGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGc            <  1:4410377/150‑1 (MQ=255)
                                                                                                                                           gCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCa           <  2:731586/150‑1 (MQ=255)
                                                                                                                                           gCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCa           <  2:440540/150‑1 (MQ=255)
                                                                                                                                                  ggCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACa    <  1:1237116/150‑1 (MQ=255)
                                                                                                                                                   gCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAg   >  1:562288/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGG  >  NC_000913/4633509‑4633806

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: