Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 243,969 G→A 5.5% Q68* (CAA→TAA) ‡ fkl ← FKBP‑type peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947243,9690GA 5.5% 1052.8 / 31.0 329Q68* (CAA→TAA) ‡fklFKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base G (155/156);  new base A (11/7);  total (166/163)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 330 total aligned reads displayed.

GTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAGGGGTGCCATCAGGC  >  NC_002947/243821‑244118
                                                                                                                                                    |                                                                                                                                                     
gTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGc                                                                                                                                                      >  2:3230000/1‑150 (MQ=255)
gTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGc                                                                                                                                                      >  1:4169924/1‑150 (MQ=255)
gTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGc                                                                                                                                                      >  2:4754782/1‑150 (MQ=255)
gTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGc                                                                                                                                                      >  2:1282900/1‑150 (MQ=255)
 ttAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGcc                                                                                                                                                     <  2:2222567/150‑1 (MQ=255)
 ttAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGcc                                                                                                                                                     >  1:964865/1‑150 (MQ=255)
  tAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGccc                                                                                                                                                    <  2:4423240/150‑1 (MQ=255)
   aGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGcccc                                                                                                                                                   <  1:3040802/150‑1 (MQ=255)
    gTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCt                                                                                                                                                  >  1:3762160/1‑150 (MQ=255)
    gTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCt                                                                                                                                                  >  2:4030656/1‑150 (MQ=255)
    gTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCt                                                                                                                                                  <  2:2339737/150‑1 (MQ=255)
    gTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCt                                                                                                                                                  <  1:1647315/150‑1 (MQ=255)
    gTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCt                                                                                                                                                  <  2:1513852/150‑1 (MQ=255)
    gTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCt                                                                                                                                                  <  2:4612117/150‑1 (MQ=255)
     ttCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTg                                                                                                                                                 <  1:1121141/150‑1 (MQ=255)
      tCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCGTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGa                                                                                                                                                <  2:2361306/150‑1 (MQ=255)
        gTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATa                                                                                                                                              <  2:5949164/150‑1 (MQ=255)
        gTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATa                                                                                                                                              <  2:161876/150‑1 (MQ=255)
        gTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATa                                                                                                                                              <  1:3560344/150‑1 (MQ=255)
        gTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATa                                                                                                                                              <  2:5476062/150‑1 (MQ=255)
        gTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGCAACCCCTGaaa                                                                                                                                              <  1:5833858/150‑3 (MQ=255)
         tGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAg                                                                                                                                             <  1:6658758/150‑1 (MQ=255)
         tGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAg                                                                                                                                             <  1:1939803/150‑1 (MQ=255)
          gCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACAGCCTGCATGCGTGCCTTGTCGAGTATCAGGGGCAACCCCTGaaact                                                                                                                                            >  1:5375423/1‑146 (MQ=255)
           cACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGa                                                                                                                                           <  2:2702329/150‑1 (MQ=255)
            aCCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGccc                                                                                                                                                    <  2:6392497/140‑1 (MQ=255)
            aCCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGAc                                                                                                                                          >  2:4524999/1‑150 (MQ=255)
            aCCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGAc                                                                                                                                          <  1:4778955/150‑1 (MQ=255)
              cccATAACGCCCCCGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTGTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGGCGAGTTTCAGGGGTTGCCCCTGCTAGGACTg                                                                                                                                        >  2:4188163/1‑150 (MQ=255)
              cccATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTg                                                                                                                                        >  1:6467351/1‑150 (MQ=255)
              cccATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTg                                                                                                                                        >  2:3104433/1‑150 (MQ=255)
                cATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCt                                                                                                                                      >  1:6331672/1‑150 (MQ=255)
                 aTAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCtt                                                                                                                                     <  2:255319/150‑1 (MQ=255)
                   accGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCa                                                                                                                                   <  2:6751635/148‑1 (MQ=255)
                     cGCCCACGTTCATTAGCCCTGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAgg                                                                                                                                 <  1:1362639/150‑1 (MQ=255)
                        ccACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCt                                                                                                                              >  1:5000106/1‑150 (MQ=255)
                        ccACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCt                                                                                                                              >  2:4763811/1‑150 (MQ=255)
                         cACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCtt                                                                                                                             <  1:1845284/150‑1 (MQ=255)
                            gTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGa                                                                                                                          >  2:35171/1‑150 (MQ=255)
                                aTTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAg                                                                                                                      >  2:1166820/1‑150 (MQ=255)
                                 ttAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGTCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGt                                                                                                                     <  2:3580053/150‑1 (MQ=255)
                                 ttAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCGGGCCTTCGACCAGt                                                                                                                     <  2:3868109/150‑1 (MQ=255)
                                  tAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTg                                                                                                                    <  2:6224065/150‑1 (MQ=255)
                                     ccATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCAt                                                                                                                 >  2:1982384/1‑150 (MQ=255)
                                      cATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATc                                                                                                                >  1:2041260/1‑150 (MQ=255)
                                      cATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATc                                                                                                                >  2:4213523/1‑150 (MQ=255)
                                        tGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCa                                                                                                              >  1:5677966/1‑150 (MQ=255)
                                           aaCGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCt                                                                                                           <  1:1227261/150‑1 (MQ=255)
                                                 tttCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGc                                                                                                     >  2:5876894/1‑150 (MQ=255)
                                                    cggcggCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTg                                                                                                  <  2:2416727/150‑1 (MQ=255)
                                                     ggcggcCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTgg                                                                                                 >  1:2539880/1‑150 (MQ=255)
                                                     ggcggcCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTgg                                                                                                 <  1:383992/150‑1 (MQ=255)
                                                      gcggcCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGc                                                                                                <  1:3230000/150‑1 (MQ=255)
                                                          ccTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCAtct                                                                                            <  2:4998540/150‑1 (MQ=255)
                                                            tGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTcc                                                                                          >  1:6362309/1‑150 (MQ=255)
                                                             gCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCt                                                                                         >  1:2204225/1‑150 (MQ=255)
                                                             gCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCt                                                                                         >  2:5365769/1‑150 (MQ=255)
                                                             gCAGCTTTTGCCCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCt                                                                                         >  2:2579802/1‑150 (MQ=255)
                                                                  ttttGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGAc                                                                                                                         <  2:5172690/113‑1 (MQ=255)
                                                                      gcgcAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGc                                                                                >  1:1432823/1‑150 (MQ=255)
                                                                      gcgcAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGc                                                                                >  1:5998888/1‑150 (MQ=255)
                                                                               cgTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCa                                                                       >  1:2779036/1‑150 (MQ=255)
                                                                                      ttCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGc                                                                <  2:5126398/150‑1 (MQ=255)
                                                                                       tCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCa                                                               <  1:4309514/150‑1 (MQ=255)
                                                                                         cTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCAcc                                                             >  1:2851722/1‑150 (MQ=255)
                                                                                         cTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCAcc                                                             >  2:5381094/1‑150 (MQ=255)
                                                                                          tGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCAccc                                                            >  1:33421/1‑150 (MQ=255)
                                                                                           gggCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCa                                                           <  2:2468408/150‑1 (MQ=255)
                                                                                               ttCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCt                                                       <  1:1290299/150‑1 (MQ=255)
                                                                                               ttCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCt                                                       >  2:3361965/1‑150 (MQ=255)
                                                                                                   tgctgctgcAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAg                                                   <  1:4465681/150‑1 (MQ=255)
                                                                                                    gctgctgcAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAgg                                                  >  2:1752900/1‑150 (MQ=255)
                                                                                                        ctgcAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTACAAGCCTGGCATCTCCTGGCGCAGGCGTTCa                                                                           >  1:6134361/1‑121 (MQ=255)
                                                                                                         tgcAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCaaa                                             >  1:1702326/1‑150 (MQ=255)
                                                                                                              caCTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTg                                        <  2:881353/150‑1 (MQ=255)
                                                                                                              caCTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTg                                        <  2:375638/150‑1 (MQ=255)
                                                                                                                    ctgcATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGgtcgtc                                  <  1:2321821/150‑1 (MQ=255)
                                                                                                                    ctgcATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGgtcgtc                                  >  2:4776609/1‑150 (MQ=255)
                                                                                                                    ctgcATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGgtcgtc                                  <  1:1523837/150‑1 (MQ=255)
                                                                                                                     tgcATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGtcgtcg                                 >  1:5139189/1‑150 (MQ=255)
                                                                                                                       cATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGAGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGtt                               >  2:2952260/1‑150 (MQ=255)
                                                                                                                          gCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTc                            <  1:3144248/150‑1 (MQ=255)
                                                                                                                             tGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGc                         >  1:4386247/1‑150 (MQ=255)
                                                                                                                             tGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGc                         <  1:3974512/150‑1 (MQ=255)
                                                                                                                               ccTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCa                       >  2:6560330/1‑150 (MQ=255)
                                                                                                                               ccTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCa                       >  1:3368849/1‑150 (MQ=255)
                                                                                                                                  tGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAgg                    <  2:5721064/150‑1 (MQ=255)
                                                                                                                                   gTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGc                   >  1:1919152/1‑150 (MQ=255)
                                                                                                                                      gAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAg                >  1:799529/1‑150 (MQ=255)
                                                                                                                                      gAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAg                >  2:4299276/1‑150 (MQ=255)
                                                                                                                                       aGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAgg               >  2:1910432/1‑150 (MQ=255)
                                                                                                                                        gTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAggg              <  1:3587688/150‑1 (MQ=255)
                                                                                                                                         tatCAGGGGCAACCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAgggg             >  1:1196552/3‑150 (MQ=255)
                                                                                                                                         tatCAGGGGCAACCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAgggg             >  1:5404909/3‑150 (MQ=255)
                                                                                                                                          ttCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAGGGGt            <  1:3721480/150‑1 (MQ=255)
                                                                                                                                           tCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAGGGGTg           >  1:2546161/1‑150 (MQ=255)
                                                                                                                                           tCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAGGGGTg           >  2:6839479/1‑150 (MQ=255)
                                                                                                                                               gggTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAGGGGTGCCAt       <  1:5912551/150‑1 (MQ=255)
                                                                                                                                                 ggtGCCCCTGAGAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCGGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAGGGGTGCCATCa     <  2:5814939/148‑1 (MQ=255)
                                                                                                                                                 gTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAGGGGTGCCATCa     >  2:4880492/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTAGTTCGTGCACCCCATAACGCCCACGTTCATTAGCCATGAAACGCGTTTCGGCGGCCTGCAGCTTTTGCGCAGCGGCGTCACCTTCCTGGGCTTCATGCTGCTGCAGCACTGCCTGCATGCGTGCCTTGTCGAGTTTCAGGGGTTGCCCCTGATAGGACTGCTTCAGGCCTTCGACCAGTGCATCCAGCTGCAAGCCTGGCATCTCCTGGCGCAGGCGTTCACCCAGGCTGGCACCCAGGCTGTAGGCCAAATCGTGGTCGTCGTTGTCGGCCAAGGCGAGGGGTGCCATCAGGC  >  NC_002947/243821‑244118

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: