Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,046,068 T→A 14.9% intergenic (‑167/+100) 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,0680TA14.9% 527.1 / 15.8 218intergenic (‑167/+100)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base T (48/125);  new base A (34/2);  total (91/127)
Fisher's exact test for biased strand distribution p-value = 2.89e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

AATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGC  >  NC_000913/3045922‑3046217
                                                                                                                                                  |                                                                                                                                                     
aattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCtctc                                                                                                                                                    <  1:2485759/150‑1 (MQ=255)
 attCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACTcc                                                                                                                                                   >  2:969701/1‑150 (MQ=255)
   tCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCGAACCCTCTcccc                                                                                                                                                 >  2:411976/1‑150 (MQ=255)
   tCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCGAACCCACTcccc                                                                                                                                                 >  1:2501962/1‑150 (MQ=255)
    cGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCa                                                                                                                                                <  2:2791463/150‑1 (MQ=255)
     gTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCAATACCCaa                                                                                                                                               >  1:1534855/1‑150 (MQ=255)
             aTCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                                                                                                                       <  2:470299/150‑1 (MQ=255)
             aTCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                                                                                                                       <  2:2620403/150‑1 (MQ=255)
              tCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                                      <  1:4245395/150‑1 (MQ=255)
              acTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                                      <  1:2978150/149‑1 (MQ=255)
               cTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGTGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                                     <  1:1122701/150‑1 (MQ=37)
               cTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                                     <  1:3156975/150‑1 (MQ=37)
                tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAggg                                                                                                                                    <  2:3132913/150‑1 (MQ=37)
                tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACACCCCAAAGGGGCGAggg                                                                                                                                    >  1:3542100/1‑150 (MQ=255)
                tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCCCCCAAACCCACACCCCCAAGGGGCGAggg                                                                                                                                    >  1:1916922/1‑150 (MQ=255)
                tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCACTCCCCAAAGGGGCGAggg                                                                                                                                    >  1:228405/1‑150 (MQ=255)
                 ttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCAAACCCACTCCCCAAAGGGGCGAgggg                                                                                                                                   >  1:3504271/1‑150 (MQ=255)
                  tGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCGGGGAGAGGGCCGGGGTGAGGGTAAATCTTCGCGCCAGTGCCGGCCTGTTCCCCTCGCCCTAACCCTCTCCCCAAAGGGGCGCggggc                                                                                                                                  >  2:2667283/1‑149 (MQ=255)
                    agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACt                                                                                                                                <  2:2955256/150‑1 (MQ=255)
                      aTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGt                                                                                                                              <  1:60717/150‑1 (MQ=255)
                       ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCTAAACCAAAAAAAAAACACAAATGGGCGAGGGGACTGTc                                                                                                                             >  2:949033/1‑150 (MQ=255)
                        tGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCt                                                                                                                            <  2:4074341/150‑1 (MQ=255)
                          cTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGa                                                                                                                          <  2:1371855/150‑1 (MQ=255)
                           tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCTCACCCGAAACCACTCCCCAAAGGGGCGAGGGGAGtgtgtggg                                                                                                                         >  1:553008/1‑148 (MQ=37)
                           tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCTAACCCTCTCCCCCAAGGGGCGGGGGGGCTGTCTGAg                                                                                                                         >  2:1171026/1‑150 (MQ=255)
                            aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGc                                                                                                                        >  2:3542339/1‑150 (MQ=255)
                              cAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCGTTCAGGGAGAGGGCCGGGGTGAGGGTAAATAGTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCAc                                                                                                                      <  1:2164982/150‑1 (MQ=255)
                                aaCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCGCTCCCCAAAGGGGCGAGGGGACTGTCTGTGCCCtt                                                                                                                    >  1:3719873/1‑150 (MQ=37)
                                  cTGGCTCAATCTCCAATGTTGGCACGCTGTGCCCCCTCGCCCTTTCCGGGAGAGGGCCGGGGTGAGGGTAAAGATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTACCCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACtttt                                                                                                                  <  2:3171296/150‑1 (MQ=25)
                                     gCTCAATCTTCAAAGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAACCAAACCCGCTCCCCCAAGGGGCGAGGGGACTGTATGAGCACGTTTGTa                                                                                                               >  2:3260389/1‑150 (MQ=255)
                                       tCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACAGTCTGAGCACGTTTGTACt                                                                                                             >  2:2378396/1‑150 (MQ=255)
                                             ttCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCAAACCCACTCCCCCAAGGGGCGAGGGGACAGTCTGAGCACTTTTGTACTTTGTCa                                                                                                       >  1:3898532/1‑150 (MQ=255)
                                              tCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGTACTGTCTGAGCACTGTTGTACTTTGTCATc                                                                                                     >  2:1581782/1‑150 (MQ=255)
                                              tCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCAAACCCTCACCCCCAAGGGGCGAGGGGACAGTCTGAGCACGTTTGTACTTTGTAAt                                                                                                      >  1:2711434/1‑150 (MQ=255)
                                              tCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCAAACCCACACCCCAAAGGGGCGAGGGGCCGGTCTGAGCACTTTTGTACTTTGTCaa                                                                                                      >  2:4097991/1‑149 (MQ=255)
                                               cAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTACATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGGGGGGCCTGTCGGAGCACGGTTGTACGTTGTCATc                                                                                                     >  1:1944965/1‑150 (MQ=25)
                                               cAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATc                                                                                                     <  1:943909/150‑1 (MQ=255)
                                                  tGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGTGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGTCTGTCTGAGCACTTTTGTACTTTGTCATCTGa                                                                                                  <  2:1257010/150‑1 (MQ=255)
                                                  tGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGa                                                                                                  <  1:2457598/150‑1 (MQ=255)
                                                   gTTGGCACGCTTTTCCCCCTCGCCCTGTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCCGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGAc                                                                                                 <  2:3941097/150‑1 (MQ=255)
                                                       gCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGTCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGAGGGGACAGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaa                                                                                             >  2:2116404/1‑150 (MQ=255)
                                                        cACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaaa                                                                                            >  2:1780277/1‑150 (MQ=255)
                                                        cACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCAATTTTGTACTTTGTGTTGTGTCTaaat                                                                                            <  1:3583759/150‑2 (MQ=25)
                                                         aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaaaa                                                                                           <  2:3669886/150‑1 (MQ=255)
                                                         aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaaaa                                                                                           <  2:2289631/150‑1 (MQ=255)
                                                         aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCACAACCAAACCCACAACCCCAAAGGGCGGGGGGCCGGTCGGGGCACATTTGTTATTTTTGATAAGACaaaaaa                                                                                           >  1:18303/1‑150 (MQ=255)
                                                          cGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAg                                                                                          <  1:1665320/150‑1 (MQ=255)
                                                          cGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAg                                                                                          <  2:1909060/150‑1 (MQ=255)
                                                             ttttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCCGTGCCGGCCGGGTCCCCGCCCCCCAACCCACACCCCCCACGGGCGCGGGGCCGGGCTGAGCCCGGTTGGTCGGTGGCAACTGCCAAAAAAGgcg                                                                                       >  1:2931050/1‑150 (MQ=255)
                                                             ttttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcg                                                                                       <  2:170279/150‑1 (MQ=255)
                                                             ttttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcg                                                                                       <  2:170408/150‑1 (MQ=255)
                                                                   cccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACACCCCAAAGGGGCGGGGGGACTGTCTGAGCACGTTTGTACTTTTTCATCTGACTAAAAAGGCGCCGAAg                                                                                 >  2:114136/1‑150 (MQ=37)
                                                                     cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTAACCCACTACCCAAAGGGGCGAGGGGCCAGTCAGAGCACATTTGTACAGTGTCATCTGACTAAAAAGGCGCCGAAgcg                                                                               >  2:2321882/1‑150 (MQ=255)
                                                                      tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcgc                                                                              <  1:3422827/150‑1 (MQ=255)
                                                                      tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcgc                                                                              <  1:368060/150‑1 (MQ=255)
                                                                      tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcgc                                                                              >  2:304489/1‑150 (MQ=255)
                                                                      tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCAAACCCACTCCCCAAAGGGGCGCGGGGCAGGTAGGAGCACATTTGTTCTTTGTCAAATGACTAAAAAAGCGCCGAAacgc                                                                              >  1:2385427/1‑150 (MQ=255)
                                                                      tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCACACCCTAACCCCCAACCCAAAGGGGCGAGGGGCCTGTCTGAGCACATTTTTAATTTTTAATCTTAAGAAAAAGGCGCCCAAgccc                                                                              >  2:819140/1‑148 (MQ=255)
                                                                      tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCACTCCCCCAAGGGGCGAGGGGACAGTCTGAGCACTTTTGTACTTTGTCATATGACTAAAAAGGCGCCGAAgcgc                                                                              >  2:3767486/1‑150 (MQ=255)
                                                                      tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCAAACCCTCTCCCCAAAGGGGCGCGGGGCCAGTCGGAGCACATTTGTACGTTGTCAACTGACAAAAAAGGCGccccagccc                                                                              >  2:1661284/1‑143 (MQ=21)
                                                                        gCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCGCACCCAAACCCACACCCCCAAGGGGCGAGGGGACCGTCTGAGCACATTTGTACTTTGTCATCTGACTAAAAATGCGCCCAAGCGccg                                                                            >  2:1520636/1‑149 (MQ=255)
                                                                          ccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCttt                                                                          <  1:2654666/150‑1 (MQ=255)
                                                                          ccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCttt                                                                          <  1:1734972/150‑1 (MQ=255)
                                                                          ccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGCTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCTTCGGCGCCttt                                                                          <  1:2544908/150‑1 (MQ=37)
                                                                            tttCAGGGAGAGGGCCGGGGTGAGGGTAAATCTTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAg                                                                        <  1:317737/150‑1 (MQ=255)
                                                                              tccgggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaa                                                                      <  2:314952/147‑1 (MQ=255)
                                                                              tCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaa                                                                      <  2:2698960/150‑1 (MQ=255)
                                                                               cAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaa                                                                     <  2:2464818/150‑1 (MQ=255)
                                                                               cAGGGAGAGGGCCGGGGTGAGGGAAAATATACGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaa                                                                     <  2:1562640/150‑1 (MQ=255)
                                                                                aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCGAACCCAATACCCAAAGGGGGGAGGGGACAGTCTGAGCACTTTTTTACTTTGTCAACTGACAAAAAAGGCGCCGACGCGCCGTTAGacaa                                                                    >  2:2873072/1‑150 (MQ=255)
                                                                                  ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                                  <  1:969701/150‑1 (MQ=255)
                                                                                           cGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCa                                                         <  1:4080449/150‑1 (MQ=255)
                                                                                                tGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGaa                                                    <  2:1321278/150‑1 (MQ=255)
                                                                                                     gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                                               <  2:374079/150‑1 (MQ=255)
                                                                                                     gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCTCTACCCAAAGGGGCGAGGGGACTGTCTGAGCACATTTTTACTTTGTCATCTTAATAAAAAGGCGCCGAAGCGCCTTTATAAAATAGTAGAATCATTGAATTACt                                               >  2:94177/1‑150 (MQ=255)
                                                                                                       accgcttCGCGCCAGTGCCGGCCGGTGCCCCTCCCCCGACCCCTCGCCCCCCCGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGGCATCTGACTAAAACGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTgg                                             <  2:3381340/145‑1 (MQ=255)
                                                                                                       aaaTATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTgg                                             <  1:325740/150‑1 (MQ=255)
                                                                                                         atatTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTa                                           <  2:325554/150‑1 (MQ=255)
                                                                                                             tCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCg                                       <  2:148899/150‑1 (MQ=255)
                                                                                                                cgcCAGTGCCGGCCTGTTACCCTCACCCAAAACCTCTCCCCAAAGGGGCGAGGGGACTGTCAGATCACTTTTGTACTTTGTAAACAGACTAAAAAGGCGCCGAAGCGCCGTTAGAAAAAAGTCGAATCAGTGAAATACAGGTATTAGCTa                                    >  2:897498/1‑150 (MQ=37)
                                                                                                                  ccAGTGCCGGCCTGTTCCCCACACCCTAAACCTCTACCCACAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACAATGTCATCTGACTAAAAACGCGCCGAAGCGCCATTAGAAAATAGTCGAATCAGTGAAATACAGGAATTCGCTaaa                                  >  2:3782094/1‑149 (MQ=255)
                                                                                                                   cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                                 <  1:1083138/150‑1 (MQ=255)
                                                                                                                         cGGCCTGTTACCCTCACCCTAACCCTCAACCCAAAGGGGCGAGGGGACAGTCTGAGCACTTTTGTACTTTGTCATATGACTAAAAAGGCGCCCAAGCGCCTTTAGAAAAAAGTAGAATAAGTGAAATAATGGTATTAGCTAATCGGTACg                           >  2:2466567/1‑150 (MQ=255)
                                                                                                                            cctctaccccACACCCTAACCCACTCACCAAAGGGGCGAGGGGACTGACAGAGCACTATTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAg                        <  2:744423/144‑1 (MQ=255)
                                                                                                                            ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAg                        <  1:3098372/150‑1 (MQ=255)
                                                                                                                              tGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGa                      <  2:175551/150‑1 (MQ=255)
                                                                                                                                ttCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGGGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGc                    >  2:3089768/1‑150 (MQ=255)
                                                                                                                                    ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGaa                <  1:3014970/150‑1 (MQ=255)
                                                                                                                                    ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGaa                <  2:255198/150‑1 (MQ=255)
                                                                                                                                     cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAAc               >  1:3070467/1‑150 (MQ=255)
                                                                                                                                     cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAAc               >  2:1247249/1‑150 (MQ=255)
                                                                                                                                     cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAAc               >  1:1976396/1‑150 (MQ=255)
                                                                                                                                      tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACa              <  2:854393/150‑1 (MQ=255)
                                                                                                                                      tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACa              <  2:311732/150‑1 (MQ=255)
                                                                                                                                          ccTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGt           <  2:2775404/149‑1 (MQ=255)
                                                                                                                                             aaCCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACg       <  2:3092785/150‑1 (MQ=255)
                                                                                                                                             aaCCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACg       <  2:3423204/150‑1 (MQ=255)
                                                                                                                                                ccTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTc    >  1:724877/1‑150 (MQ=255)
                                                                                                                                                  tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc  >  2:698370/1‑150 (MQ=255)
                                                                                                                                                  tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc  >  2:2932365/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGC  >  NC_000913/3045922‑3046217

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: