Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,046,075 A→C 5.2% intergenic (‑174/+93) 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,0750AC 5.2% 1335.2 / 7.3 431intergenic (‑174/+93)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base A (161/232);  new base C (30/2);  total (192/239)
Fisher's exact test for biased strand distribution p-value = 1.78e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 443 total aligned reads displayed.

GTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTAAAC  >  NC_000913/3045927‑3046224
                                                                                                                                                    |                                                                                                                                                     
gTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaa                                                                                                                                                      >  2:8508290/1‑150 (MQ=255)
  tGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                                                                                                    <  1:260021/150‑1 (MQ=255)
   gACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCAAACCCACACCCCAAAgg                                                                                                                                                   >  2:3877673/1‑150 (MQ=255)
       tATCTTTGAGATTGCTAACAATCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCg                                                                                                                                               <  2:6844784/150‑1 (MQ=255)
       tATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCGAACCCACACCCCACAGGGGCg                                                                                                                                               >  1:8067758/1‑150 (MQ=255)
          cTTTGAGATTGCTAACAAACTGTCTCAATCTTCAATGTTGGCACGCTTTTCCCCCGCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                                            <  2:1003234/150‑1 (MQ=37)
           tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                           >  2:1460484/1‑150 (MQ=255)
           tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTGTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACAACCCAAAAGGGCGAggg                                                                                                                                           >  2:1929157/1‑150 (MQ=255)
              gagaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGAc                                                                                                                                        <  1:1977161/150‑1 (MQ=255)
               agtttGATAACAATCTGGCTCAATCTTAAATGTTGGCACGCTTTTCCGCCTGGCCCTTTCCGGGTGCGGGCCGGGGTGAGGGTAAATATCCGCGCCCGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACt                                                                                                                                       <  1:2257674/147‑1 (MQ=25)
               agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCAAAACCACAACCCAAAGGGGCGAggggccg                                                                                                                                       >  1:2556001/1‑147 (MQ=255)
                gaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCGAACCCACTCCCCCAAGGGGCGAGGGGACTg                                                                                                                                      >  2:6100324/1‑150 (MQ=255)
                   tGCTAACAACCTGGCTCAATCTTCTTTGTTGGCACGCTTTTCCCCCTCGCCCTTTCTGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCt                                                                                                                                   <  1:7829526/150‑1 (MQ=255)
                   tGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACTCCCCACAGGGGCGAGGGGACTGTcg                                                                                                                                   >  2:4255304/1‑149 (MQ=255)
                      tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTGACCCTAACCCTCTCCCCAAAGGGGGGTGGGGAGTGTCTgtg                                                                                                                                >  2:2837073/1‑148 (MQ=255)
                      tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCAAACCCTCTCCCCCAAGGGGCGAGGGGACTGTCTGAg                                                                                                                                >  2:2928565/1‑150 (MQ=255)
                       aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCTGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGc                                                                                                                               >  2:8257635/1‑150 (MQ=255)
                       aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCAAACCCACACCCCAAACGGGCGAGGGGACTGTCTGAGc                                                                                                                               >  1:5887472/1‑150 (MQ=255)
                        aCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCACCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAACCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCCGAgga                                                                                                                              >  2:1315110/1‑148 (MQ=255)
                         cAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCAc                                                                                                                             <  2:1747706/150‑1 (MQ=255)
                              tGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTg                                                                                                                        <  2:621376/150‑1 (MQ=255)
                               ggCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGt                                                                                                                       <  1:579552/150‑1 (MQ=255)
                                    aaTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCCCACCCTAACCCTCGCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACttt                                                                                                                  <  2:5581265/150‑1 (MQ=255)
                                      tCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGt                                                                                                                <  2:2647156/150‑1 (MQ=255)
                                      tCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCCCACCCGAACCCACACCCCAAAGGGGCGCGGGGCCTGTCTGAGCACATTTTTaagtttt                                                                                                                >  2:3013501/1‑144 (MQ=255)
                                       cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTc                                                                                                               <  2:7558705/150‑1 (MQ=255)
                                        ttCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTTAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCa                                                                                                              <  1:1333824/150‑1 (MQ=255)
                                        ttCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCGAACCCTCTCCCCACAGGGGCGAGGGGCCGGTCTGTGCACTTTTGTACTTTGTCa                                                                                                              >  1:8554062/1‑150 (MQ=255)
                                         tCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCGCACCCAAACCCACTCCCCAAAGGGGGGGGGGGACTGGGTGTGCACATTTGTTGTTTTTCAt                                                                                                             >  2:2267599/1‑150 (MQ=255)
                                         tCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCAt                                                                                                             <  1:6724353/150‑1 (MQ=255)
                                         tCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTTTCAt                                                                                                             >  2:6019540/1‑150 (MQ=255)
                                          cAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATc                                                                                                            <  2:3493990/150‑1 (MQ=255)
                                           aaTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCt                                                                                                           <  1:5537440/150‑1 (MQ=255)
                                             tGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGa                                                                                                         <  1:1006960/150‑1 (MQ=255)
                                              gTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCCCACCCAAACCCACACCCCCAAGGGGCGGGGGGCCTGTCGGGGCCCGTTTGTACGTTGTGATCTggc                                                                                                        >  2:1645969/1‑148 (MQ=255)
                                                    aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCCCACCCAAAACCACAACCCCAACGGGCGGGGGGGCTGTGTGAGCACAGTTGTACTTTGTCAACTGACTaaaaa                                                                                                  >  2:3792160/1‑150 (MQ=255)
                                                    aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGAGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCAGAGCACTTTTGTACTTTGTCATCTGACTaaaaa                                                                                                  <  2:2508413/150‑1 (MQ=255)
                                                         tttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcgc                                                                                             <  2:4346897/150‑1 (MQ=255)
                                                          ttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCTCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACGTTTGTACTTTGTCATCTGACTAAAAAGGCGcc                                                                                            >  1:6454513/1‑150 (MQ=255)
                                                            cccccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGCCTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGa                                                                                          >  2:7899554/1‑150 (MQ=255)
                                                               ccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgc                                                                                       <  1:385545/150‑1 (MQ=255)
                                                                cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCAAACCCACACCCCAAAGGGGCGAGGGGACGGTCTGAGCACATTTGTTATTTTTAATCTGTCTAAAAATGCGCCCAAgcc                                                                                      >  1:5228259/1‑149 (MQ=21)
                                                                 tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACAACCCAAAGGGGCGAGGGGACAGTCTGAGCACATTTGTACGTTTTGATCTGACTAAAAAGGCGCCGAAgcgc                                                                                     >  2:1846144/1‑150 (MQ=25)
                                                                    cccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCtt                                                                                  <  2:2100739/150‑1 (MQ=255)
                                                                     ccTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCttt                                                                                 <  2:2797582/150‑1 (MQ=255)
                                                                       tttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCACACCCAAACCCACAAACCCAAAGGGCGCGGGGCCAGTCGGAGCACATTTGTTATTTTTGATAATTAAAAAAATGCGCCCAAGCGCCCTTAg                                                                               >  1:6961030/1‑150 (MQ=255)
                                                                        ttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGa                                                                              <  2:5418290/150‑1 (MQ=255)
                                                                         tCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAAACCAAACCCCCAAGGGGCGAGGGGCCGGTCTGAGCACGTTTGTAATTTGTGATATGACAAAAAATGCGCCCAAGCGCCGTTAGaa                                                                             >  1:3582478/1‑150 (MQ=255)
                                                                           aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                                           <  1:1089633/150‑1 (MQ=255)
                                                                             ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                                         <  1:1778049/150‑1 (MQ=255)
                                                                             ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                                         <  1:5851422/150‑1 (MQ=255)
                                                                             ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTCAGAAAATa                                                                         <  2:1126171/150‑1 (MQ=255)
                                                                             ggAGAGGGCCGGGGTGAGGGTAAATAGTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                                         <  2:575032/150‑1 (MQ=255)
                                                                              gagaGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAg                                                                        >  2:4774324/1‑150 (MQ=255)
                                                                                gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTACCCCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTc                                                                      <  2:7692303/150‑1 (MQ=255)
                                                                                   ggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCAAACCCAAACCCAATCCCCATAGGGGCGAGGGGGCTGTCTGAGCCCATTTGTACTTTGGCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAAGAGTAGaa                                                                   >  2:6982404/1‑150 (MQ=255)
                                                                                    gCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTAACCCAAACCCTCACACCAAAAGGGCCAGGGGCCTGTATGAGCACGTTTTTACTTTGTAATCTGACTAAAAAGGCGCCGCAGCCCCGTTAGAAAATAgaagaaa                                                                  >  2:7491902/1‑144 (MQ=255)
                                                                                        gggTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACACCCCAAAGGGGCGCGGGGACTGTCTGAGCACATTTGTACTTTGTAATATGACAAAAAAGGCGCCGAACCGCCGTTAGAAAAAAGTCGaaaaagt                                                              >  2:4447131/1‑145 (MQ=255)
                                                                                         ggTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGTATCAGTg                                                             >  2:547672/1‑150 (MQ=255)
                                                                                               ggTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTAc                                                       <  2:8651735/150‑1 (MQ=255)
                                                                                               ggTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACTACCCAAAGGGGCGAGGGGACTGTCTGAGCACGTTTGTACTTTTTCATCTGACTAAAAAGGCGCCGAAGCGCCGTTAGAAAATAGTCGAATCAGTGAATTAc                                                       >  2:4660262/1‑150 (MQ=39)
                                                                                                gTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                                                      <  2:5276148/150‑1 (MQ=255)
                                                                                                    atatTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAAACCAATCCCCCAAGGGGCGAGGGGCCAGTCTGAGCACATTTGTACTTTTTCATCTGACAAAAAAGGCGCCCAAGCGCCGTTAGAAAATAGTCGAAACAGTGAAATACTGGTa                                                  >  2:8053989/1‑150 (MQ=255)
                                                                                                    atatTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTa                                                  <  2:5873959/150‑1 (MQ=255)
                                                                                                      atTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTAtt                                                <  1:4165133/150‑1 (MQ=255)
                                                                                                          gcgcCAGTGCCGGCCTGTTCCCCTCACCCTAAACCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTAATGGTATTCGCt                                            >  2:6479949/1‑150 (MQ=255)
                                                                                                              cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                                        <  1:3957430/150‑1 (MQ=255)
                                                                                                              cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                                        <  2:237060/150‑1 (MQ=255)
                                                                                                              cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                                        <  1:4406465/150‑1 (MQ=255)
                                                                                                                gTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCgg                                      <  2:3352954/150‑1 (MQ=255)
                                                                                                                    cGGCCTGTTCCCTTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAATAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACg                                  >  1:5375605/1‑150 (MQ=255)
                                                                                                                        cTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAgg                              >  2:4193786/1‑150 (MQ=255)
                                                                                                                         tGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGa                             <  2:5899524/150‑1 (MQ=255)
                                                                                                                          gTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAg                            >  2:6602043/1‑150 (MQ=255)
                                                                                                                               ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGaa                       >  1:7577589/1‑150 (MQ=255)
                                                                                                                                cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAAc                      >  1:7360699/1‑150 (MQ=255)
                                                                                                                                 tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACa                     >  2:7299536/1‑150 (MQ=255)
                                                                                                                                 tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACa                     <  1:5028589/150‑1 (MQ=255)
                                                                                                                                     ccTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGtt                 >  1:657831/1‑150 (MQ=255)
                                                                                                                                       tAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTAc               <  1:6321303/150‑1 (MQ=255)
                                                                                                                                       tAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTAc               <  1:3056041/150‑1 (MQ=255)
                                                                                                                                       gaaCCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTAc               <  2:8291121/149‑1 (MQ=255)
                                                                                                                                        aaCCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACg              <  2:3994420/150‑1 (MQ=255)
                                                                                                                                         aCCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACgg             <  1:797934/150‑1 (MQ=255)
                                                                                                                                             tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc         >  1:4918132/1‑150 (MQ=255)
                                                                                                                                             tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc         >  2:5107339/1‑150 (MQ=255)
                                                                                                                                             tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc         >  1:2851035/1‑150 (MQ=255)
                                                                                                                                              ctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGcc        <  1:3214033/150‑1 (MQ=255)
                                                                                                                                              ctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGcc        <  1:8130244/150‑1 (MQ=255)
                                                                                                                                              ctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGcc        <  2:2560927/150‑1 (MQ=255)
                                                                                                                                                ccccAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGt      <  1:871313/150‑1 (MQ=255)
                                                                                                                                                 cccAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTa     >  2:7383224/1‑150 (MQ=255)
                                                                                                                                                 cccAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTa     >  1:1245775/1‑150 (MQ=255)
                                                                                                                                                  ccAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTTCAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTaa    <  1:5165985/150‑1 (MQ=255)
                                                                                                                                                  ccAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTaa    >  1:16134/1‑150 (MQ=255)
                                                                                                                                                  ccAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTaa    >  1:4018888/1‑150 (MQ=255)
                                                                                                                                                  ccAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTaa    >  1:6146820/1‑150 (MQ=255)
                                                                                                                                                  ccAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTaa    >  2:3142682/1‑150 (MQ=255)
                                                                                                                                                    aaaGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTAAAc  >  2:3576452/1‑150 (MQ=255)
                                                                                                                                                    aaaGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTCAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTacac  >  2:7155115/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGCCGTAAAC  >  NC_000913/3045927‑3046224

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: