Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,046,056 T→A 8.5% intergenic (‑155/+112) ygfF ← / ← gcvP putative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,046,0560TA 8.5% 2086.6 / 20.6 816intergenic (‑155/+112)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base T (221/440);  new base A (81/3);  total (372/444)
Fisher's exact test for biased strand distribution p-value = 3.28e-31
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 829 total aligned reads displayed.

TTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACA  >  NC_000913/3045907‑3046205
                                                                                                                                                     |                                                                                                                                                     
   cctacaaaaaTAAATTCGTTGACGTAACTTTGAGATTGCAAAAAAACTGGCTAAATCTTCAATGTTGGCACGCTTTTCCCCCACGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAc                                                                                                                                                    <  1:7648209/145‑1 (MQ=255)
    cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGGGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACAcc                                                                                                                                                   >  1:4058080/1‑150 (MQ=255)
    cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                                                                                   >  1:7806772/1‑150 (MQ=255)
      agaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCCGTGCCGGCCTGTTCCCCGCAccca                                                                                                                                                 >  1:13535973/1‑149 (MQ=255)
      agaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCt                                                                                                                                                 >  2:11444597/1‑150 (MQ=255)
      agaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCAccca                                                                                                                                                 >  2:14429804/1‑149 (MQ=255)
      agaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCACAccca                                                                                                                                                 >  1:3253816/1‑149 (MQ=255)
        aaaaaTTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTaa                                                                                                                                               >  2:14065333/1‑150 (MQ=255)
        aaaaaTTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTaa                                                                                                                                               >  2:7534048/1‑150 (MQ=255)
          aaaTTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAAcc                                                                                                                                             >  1:9193753/1‑150 (MQ=255)
           aattaattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAccc                                                                                                                                            <  2:7051448/150‑1 (MQ=255)
            attaattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAcccg                                                                                                                                           >  1:14760601/1‑149 (MQ=255)
            attaattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCct                                                                                                                                           >  1:8989601/1‑150 (MQ=255)
               aattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCtctc                                                                                                                                        <  2:15349430/150‑1 (MQ=255)
                attCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTcc                                                                                                                                       <  2:13427266/150‑1 (MQ=255)
                  tCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCGCTcccc                                                                                                                                     >  1:9916058/1‑150 (MQ=255)
                    gTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAAACCACAccccca                                                                                                                                   >  1:14402846/1‑148 (MQ=255)
                    gTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCGCACCCCAACCCACACCCCaa                                                                                                                                   >  1:6892235/1‑150 (MQ=255)
                    gTTGACGTATCTTTGAGATTGCTAACAAACTGACTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACCCCCAAACCCGCAccccca                                                                                                                                   >  1:251696/1‑148 (MQ=255)
                     ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                                  <  1:15454063/150‑1 (MQ=255)
                       gACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAAgg                                                                                                                                >  1:11709163/1‑150 (MQ=255)
                        aCGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAggg                                                                                                                               <  2:12251171/150‑1 (MQ=255)
                               tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGTCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCAAAACCACAACCCAAAGGGGCGAggg                                                                                                                        >  1:4285809/1‑150 (MQ=37)
                               tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGCggg                                                                                                                        >  2:1014830/1‑150 (MQ=255)
                               tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCACCCTCGCCCTTTCAGGGAGAGGGCCGCGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTACCCGCCCCCAAAACCACACCCCCCAGGGGCGAggg                                                                                                                        >  1:6046929/1‑150 (MQ=255)
                                ttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCACTCCCCCAAGGGGCGAgggg                                                                                                                       >  2:12480792/1‑150 (MQ=255)
                                 tGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                                      <  2:8861498/150‑1 (MQ=255)
                                       tGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAAGGGGCGAGGGGGCTGTCt                                                                                                                >  2:6518978/1‑150 (MQ=255)
                                         cTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTGTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGAGGGGACTGTCTGa                                                                                                              <  2:9133738/150‑1 (MQ=255)
                                          tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                                                                                                             <  1:5206760/150‑1 (MQ=255)
                                          tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCTCACCCAAAACCACAACCCAAAGGGGCGAGGGGACAGTCTGAg                                                                                                             >  1:8436283/1‑150 (MQ=255)
                                          tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCGAACCCTCTCCCCCAAGGGGCGAGGGGACTGTCTGAg                                                                                                             >  1:9721203/1‑150 (MQ=255)
                                          tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCAAACCCACTCCCCAAACGGGCGCGGGGCCTGTCTGAg                                                                                                             >  2:5476552/1‑150 (MQ=255)
                                           aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTGTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGc                                                                                                            <  2:14776203/150‑1 (MQ=255)
                                              aaaCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACt                                                                                                         <  2:421532/150‑1 (MQ=255)
                                              aaaCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAACCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACt                                                                                                         >  2:13602434/1‑150 (MQ=255)
                                                    gCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTa                                                                                                   <  1:6477416/150‑1 (MQ=255)
                                                       cAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCCCACCCAAACCCACAACCCAAACGGGCCGGGGGGCAGTCGGAGCACAGTTGttcat                                                                                                >  2:6739783/1‑146 (MQ=255)
                                                        aaTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACttt                                                                                               <  2:10855109/150‑1 (MQ=255)
                                                        aaTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgggttagg                                                                                                                   <  2:3707447/130‑6 (MQ=25)
                                                         aTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTg                                                                                              <  1:15187545/150‑1 (MQ=255)
                                                           cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTc                                                                                            <  1:8076696/150‑1 (MQ=255)
                                                            ttCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCa                                                                                           >  2:2471483/1‑150 (MQ=255)
                                                             tCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGTGCACGTTTGTACGTTGTCAt                                                                                          >  2:2150591/1‑150 (MQ=255)
                                                             tCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAAACCACTACCCCAAGGGGCGGGGGGCCGGTCTGAGCACTGTTGTACATTTTCAt                                                                                          >  1:12368679/1‑150 (MQ=255)
                                                              cAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACAACCCCAAGGGGCGAGGGGACTGTGTGAGCACGTTTTTAAtttttaaaa                                                                                         >  1:4784223/1‑144 (MQ=255)
                                                               aaTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCt                                                                                        <  2:3089594/150‑1 (MQ=255)
                                                               aaTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCTAACCCGCACCCCACAGGGGGGGGGGGCCCGTCTGGGCCCGTTTGTGCGATGTCATCt                                                                                        >  2:11774025/1‑150 (MQ=255)
                                                                      gCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCACCCCAAACGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaa                                                                                 >  2:4909711/1‑150 (MQ=255)
                                                                       cACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACGTTGTCCTCTGACAaaaa                                                                                >  2:15030143/1‑150 (MQ=255)
                                                                          gCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCCCCCCAACCCGCTCCCCCCAGGGGCGAGGGGACGGGCTGCGCACGGTTGTCCTTGGTCATCAGACTAAAAAgg                                                                             >  1:1671113/1‑150 (MQ=255)
                                                                              ttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGcc                                                                         <  2:9821850/150‑1 (MQ=255)
                                                                              ttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGcc                                                                         <  1:3666433/150‑1 (MQ=255)
                                                                                  cccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAg                                                                     <  1:2865703/150‑1 (MQ=255)
                                                                                    cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcg                                                                   >  2:1040382/1‑150 (MQ=255)
                                                                                     tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCGAACCCACTCCCCCAAGGGGCGCGGGGCCAGTCGGCGCACCTTTGTACGTTGTCATCGGACTAAAAAGGCGCCCCAgcgc                                                                  >  2:4723176/1‑150 (MQ=255)
                                                                                     ccGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcgc                                                                  <  1:5679010/149‑1 (MQ=255)
                                                                                       gCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCAAACCCACACCCCCAAGGGGCGGGGGGACTGTCGGAGCACATTTGTACTTTTTGATATGACAAAAAAAGCGCCCAAGCGCCt                                                                >  1:10809737/1‑150 (MQ=255)
                                                                                        cccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCtt                                                               <  1:9505243/150‑1 (MQ=255)
                                                                                         ccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCttt                                                              <  2:6132250/150‑1 (MQ=255)
                                                                                              cAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaa                                                         >  2:9643203/1‑150 (MQ=255)
                                                                                              aaGGGCGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaa                                                         <  1:15045744/149‑1 (MQ=255)
                                                                                               aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                        <  2:6307999/150‑1 (MQ=255)
                                                                                               aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                        <  2:10624497/150‑1 (MQ=255)
                                                                                               aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAAACCTCACCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTAATTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                        >  2:1316942/1‑150 (MQ=255)
                                                                                               aGGGAGAGGGCCGGGGAGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAAACCTCTCCCCAAACGGGCGAGGGGACTGTCGGAGCACTTTTGTACGTTGTCATCAGACGAAAAAGGCGCCGAAGCGCCGTTAGcaaa                                                        >  2:14541916/1‑150 (MQ=255)
                                                                                                    gagGGCCGGGGTGAGGGTACATATTCGCGCCAGTGCCGGCCTGTTCCCCTCCCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTc                                                   <  2:12083753/150‑1 (MQ=255)
                                                                                                    gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTc                                                   <  1:823893/150‑1 (MQ=255)
                                                                                                       ggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGaa                                                >  2:8073539/1‑150 (MQ=255)
                                                                                                        gCCGGGGTGAGCGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCTCTCCCCAAAGGCGCGAGGCGACTGTCTCAGCACTTTTGTACTTTGTCAACTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAAt                                               <  1:1396188/150‑1 (MQ=38)
                                                                                                            gggTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGt                                           <  1:1054778/150‑1 (MQ=255)
                                                                                                                gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAAt                                       <  1:10021674/150‑1 (MQ=255)
                                                                                                                    gTAACTAGTCGCGCCAGTGCCGGCCTGCTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGTCTAATAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGACTCAGTGAATTACt                                   <  2:8257763/150‑1 (MQ=38)
                                                                                                                    gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                                   <  1:857530/150‑1 (MQ=255)
                                                                                                                    gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCGCTCCCCAAAGGGGCGTGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                                   >  2:1283065/1‑150 (MQ=255)
                                                                                                                    gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTAACCCTCAACCCAAAGGGGCGAGGGGCCTGTCTGAGCACTTTTGTACTTTTTAATAATTCTAAAAAGGCGCCGAAGCGCCGTTAGAAAATAGTCGAATAAGTGAAATACt                                   >  2:3328111/1‑150 (MQ=255)
                                                                                                                     tAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                         >  1:12909382/1‑63 (MQ=18)
                                                                                                                       aaTATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCGAACCCACACCCCAAACGGGCGAGGGGCCAGTCGGAGCACTTTTGTACTTTGTCATATGACTAAAAAGGCGCCGAAGCGCCGTTAGAAAATAGTCGAATCAGTGAATTaattgt                                >  2:11068744/1‑145 (MQ=255)
                                                                                                                        atatTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTa                               <  1:1913469/150‑1 (MQ=255)
                                                                                                                         tatTCGCGCCAGTGCCGGCCTGTTCCCCGCAACCTAACCCACACCCCAAACGGGCGAGGGGACAGTCTTAGCACATTTTTACTTTTTAATCTGACTAAAAAGGCGCCGAAGAGCCGTTAGAAAAAATTAGAATAAGTTAATTACTTGTAt                              >  2:6903954/1‑150 (MQ=255)
                                                                                                                            tCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCg                           <  2:2020928/150‑1 (MQ=255)
                                                                                                                             cgcgCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGc                          <  1:6383994/150‑1 (MQ=255)
                                                                                                                             cgcgCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGc                          <  2:5632367/150‑1 (MQ=255)
                                                                                                                              gcgcCAGTGCCGGCCTGTTCCCCTCACCCAAAACCACTACCCAAAGGGGCGGGGGGGCAGTAAGAGCACATTTGTAATTTTTGAAATGTCTAAAAATGCGCCCAAACGCCATTATAAAAAAATCTAAAAACTGAATAAATGTTAttagat                         >  2:1550714/1‑146 (MQ=255)
                                                                                                                               cgcCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTa                        <  2:4064449/150‑1 (MQ=255)
                                                                                                                                  cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                     <  2:11310600/150‑1 (MQ=255)
                                                                                                                                  cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                     <  2:5164256/150‑1 (MQ=255)
                                                                                                                                  cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                     <  2:2322982/150‑1 (MQ=255)
                                                                                                                                  cAGTGCCGGCCTGTTCCCCTCACCCAAACCCACTACCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTTTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATAAGTGAAATACAGGTATTCGCTaaac                     >  2:1730981/1‑148 (MQ=255)
                                                                                                                                      gCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTa                 <  2:7789048/150‑1 (MQ=255)
                                                                                                                                        cGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACg               <  2:7298497/150‑1 (MQ=255)
                                                                                                                                         ggCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGc              <  1:1317114/150‑1 (MQ=255)
                                                                                                                                         ggCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGc              <  1:11112097/150‑1 (MQ=255)
                                                                                                                                         ggCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGc              <  1:2223047/150‑1 (MQ=255)
                                                                                                                                          gCCTGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCa             <  2:5827803/150‑1 (MQ=255)
                                                                                                                                             tggtCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGa          <  2:4594247/147‑1 (MQ=255)
                                                                                                                                               ttCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGc        <  2:95873/150‑1 (MQ=255)
                                                                                                                                                 ccccTCCCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAg      <  2:9329760/150‑1 (MQ=255)
                                                                                                                                                   ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGaa    >  1:7384287/1‑150 (MQ=255)
                                                                                                                                                    cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTGTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAAc   >  1:11392812/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACA  >  NC_000913/3045907‑3046205

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