Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,046,068 T→A 15.7% 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,0680TA15.7% 1129.7 / 55.8 486intergenic (‑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 (101/285);  new base A (81/2);  total (199/287)
Fisher's exact test for biased strand distribution p-value = 8.77e-37
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 519 total aligned reads displayed.

TAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGC  >  NC_000913/3045921‑3046217
                                                                                                                                                   |                                                                                                                                                      
 aattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCtctc                                                                                                                                                     <  2:7319792/150‑1 (MQ=255)
   ttCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAAACCACAccc                                                                                                                                                   >  2:2882420/1‑150 (MQ=255)
       ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                                               <  2:778231/150‑1 (MQ=255)
          aCGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGcccctttggg                                                                                                                                            <  1:4058802/150‑7 (MQ=255)
            gTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGc                                                                                                                                          >  2:8726722/1‑150 (MQ=255)
               tCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGACACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                                       <  2:3395265/150‑1 (MQ=255)
                cTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                                      <  1:8789979/150‑1 (MQ=255)
                cTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                                      <  1:4279293/150‑1 (MQ=255)
                cTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                                      <  1:1030421/150‑1 (MQ=37)
                cTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                                      <  1:2908880/150‑1 (MQ=37)
                cTTTGAGATTGCTAACAAACTGGCACAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                                      <  1:5469439/150‑1 (MQ=255)
                 tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTCACCCTCTCCCCCAAGGGGCGAggg                                                                                                                                     >  2:8496199/1‑150 (MQ=255)
                 tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                     <  2:9468696/150‑1 (MQ=255)
                     agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACt                                                                                                                                 >  2:2361011/1‑150 (MQ=255)
                       agtGCTAACACACTGGCTCAATCTTCACTGTTGGCACGCTTGTCCCCCTCGCCCTGGCAGGGAGCGGGCCGGGGTGAGGGTAAATCGTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGt                                                                                                                               <  1:3550233/148‑1 (MQ=25)
                       aTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGt                                                                                                                               <  2:4511965/150‑1 (MQ=255)
                        ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACGGTc                                                                                                                              >  2:5202120/1‑150 (MQ=255)
                          gCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCTCTACCCAAAGGGGCGAGGGGACTGTCTg                                                                                                                            >  2:2230803/1‑150 (MQ=255)
                           cTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGa                                                                                                                           <  2:8299769/150‑1 (MQ=255)
                            tAACAAACTGGCTCAATCTTCCATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                                                                                                                          >  2:4077251/1‑150 (MQ=255)
                             aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGc                                                                                                                         >  2:8774638/1‑150 (MQ=255)
                             aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCAAAACCACTCCCCCAAGGGGCGCGGGGGCTGTCTGAGc                                                                                                                         >  2:7443611/1‑150 (MQ=255)
                              aCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCa                                                                                                                        >  2:8085406/1‑150 (MQ=255)
                              aCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCa                                                                                                                        >  2:9568352/1‑150 (MQ=255)
                                   cTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACtttt                                                                                                                  <  1:6488708/150‑1 (MQ=255)
                                       cTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTAc                                                                                                              <  2:1732500/150‑1 (MQ=255)
                                       cTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAAACCACTACCCCAAGGGGCGAGGGGAATGTCTGAGCACATTTTTaa                                                                                                              >  2:6805975/1‑149 (MQ=255)
                                         cAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCCCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACtt                                                                                                            <  2:1602909/150‑1 (MQ=255)
                                         cAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACtt                                                                                                            <  1:7206988/150‑1 (MQ=255)
                                         cAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACtt                                                                                                            <  1:2149460/150‑1 (MQ=255)
                                         cAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACtt                                                                                                            <  2:4122151/150‑1 (MQ=255)
                                         cAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCACTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACtt                                                                                                            <  1:8641591/150‑1 (MQ=255)
                                          aaTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCGAACCCACACCCCAAAGGGGCGCGGGGGCGGGCGGAGCCCTTTTGTTCttt                                                                                                           >  2:466495/1‑150 (MQ=255)
                                             cTTCAATGTTGGCACGCTTTTCCCCCCCGCCCTGTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTc                                                                                                        <  1:14/150‑1 (MQ=255)
                                                 aaTGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCt                                                                                                    <  2:3919994/150‑1 (MQ=255)
                                                   tGTTGGCACGCTTGTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGa                                                                                                  <  1:2526284/150‑1 (MQ=255)
                                                     ttGGCTCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACt                                                                                                <  1:1852355/150‑1 (MQ=255)
                                                       ggCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaa                                                                                              >  2:5896038/1‑150 (MQ=255)
                                                        gCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaa                                                                                             >  2:560236/1‑150 (MQ=255)
                                                        gCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaa                                                                                             <  1:1684159/150‑1 (MQ=255)
                                                        gCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGATGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaa                                                                                             <  2:7992592/150‑1 (MQ=255)
                                                         cACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGATCACTTTTGTACTTTGTCATCTGACTaaaa                                                                                            >  2:3189643/1‑150 (MQ=255)
                                                         cACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaaa                                                                                            <  1:8211358/150‑1 (MQ=255)
                                                          aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaaaa                                                                                           <  1:1203374/150‑1 (MQ=255)
                                                          aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaaaa                                                                                           <  1:7667104/150‑1 (MQ=255)
                                                              ttttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcg                                                                                       >  2:5093594/1‑150 (MQ=255)
                                                              ttttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACAACCCAAAAGGGCGATGGGACTGTCTGTGCACAGTTGTACTATGTAATATGAATaaaaattcg                                                                                       >  2:271757/1‑146 (MQ=255)
                                                              ttttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCGAACCCACTCCCCAAACGGGCGAGGGGACTGGCGGAGCACTTTTGTACTTTGTCCTCTTTCAAAAAATgcg                                                                                       >  2:6615004/1‑150 (MQ=255)
                                                              ttttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACCCCCTCACCCTCTCCCCCAAGGGGCGGGGGGCCTGTCTGAGCACTTGTGTACTTTGTCATATGAATAAAAAGgcg                                                                                       >  2:4255243/1‑150 (MQ=255)
                                                              ccccggcccTCTCCCTGAAAGGGCGAGGGGCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACTGTCTGAGCACGTTTGTACTTTGTCATCTGACTAAAAAGgcg                                                                                       >  2:7907806/7‑150 (MQ=25)
                                                              ccccggcccTCTCCCTGAAAGGGCGAGGGGCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCGAACCCTCTCCCCAAAGGGGCGCGGGGGCGGTCTGTGCACGGTTGTACGTTGTCATCTTTCTAAAAATgcg                                                                                       >  2:2059050/7‑150 (MQ=12)
                                                               tttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcgc                                                                                      <  1:9007131/150‑1 (MQ=255)
                                                                   ccccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCTAACCCTCTCCCCCAAGGGGCGCGGGGACTGTCTGAGCACATTTGTACGTTGTCATCGGACTAAAAAGGCGCCGaa                                                                                  >  2:150675/1‑150 (MQ=255)
                                                                    cccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAg                                                                                 <  1:3490861/150‑1 (MQ=255)
                                                                     ccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTGCCCCTCACCCTAACCCTCTCCCCATAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgc                                                                                <  1:7451588/150‑1 (MQ=255)
                                                                      cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCCCTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcg                                                                               >  2:5995819/1‑150 (MQ=255)
                                                                      cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTAACCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTCCGTTGTCATCTGACTAAAAAGGCGCCGAAgcg                                                                               >  2:4894394/1‑150 (MQ=255)
                                                                         gCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCGCACCCAAACCCACAACCCCAAGGGGCGAGGGGGCTGGGGGAGCACATTTGTTCGTTTTGATATGAATAAAAAAGCGCCGAAGCGccg                                                                            >  1:8493585/1‑149 (MQ=255)
                                                                                 aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                                    <  2:9223894/150‑1 (MQ=255)
                                                                                   ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                                  <  2:6137617/150‑1 (MQ=255)
                                                                                   ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCACTACCCAAAGGGGCGAGGGGACAGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCCTAGCGCCATTAGaaaaaa                                                                  >  2:3816786/1‑148 (MQ=37)
                                                                                   ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCACACCCCAAACCACAACCCAAAGGGGCGCGGGGCCAGTCAGAGCACAGTTGTACATTTTGATAAGAATAAAAATGCGCCGAAGCGCCGTTAGAAAATa                                                                  >  1:9658880/1‑150 (MQ=255)
                                                                                   ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCACACCCAAAACCACGCCCCAAAGGGGCGCGGGGACTGTCTGAGGACGTTTGTTATTTGTAAACTGACAAAAAAAGCGCCCAAGCGCCGGTAGaaaaaa                                                                  >  1:5347762/1‑148 (MQ=255)
                                                                                   ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGCGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                                  <  2:5802728/150‑1 (MQ=255)
                                                                                    gagaGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGCCTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTGTAGAAAATAg                                                                 >  2:2931888/1‑150 (MQ=255)
                                                                                     agagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCAAACCCTCACCCCAAAGGGGCGCGGGGACGGTCGGAGCCCGGTTGTACGTTGTCATCTGACGAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGt                                                                >  2:5421813/1‑150 (MQ=255)
                                                                                      gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTc                                                               <  2:722171/150‑1 (MQ=255)
                                                                                        gggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGa                                                             <  2:5500694/150‑1 (MQ=255)
                                                                                           ccGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACAACCAAAACCACAACCCCAAGGGGCGCGGGGGAAGTGTGAGCACATTTTTTAATTTTCAAATTAAGAAAAAATCGCCGAAGCGCCGTTAGAAAATATTCGaaaa                                                          >  1:8053443/1‑148 (MQ=255)
                                                                                                gTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCGTTAGAAAATAGTCGAATCAGTGa                                                     >  2:8870820/1‑150 (MQ=255)
                                                                                                     ggTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTAc                                                <  2:2577864/150‑1 (MQ=255)
                                                                                                     ggTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCTCAACCCCAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTAc                                                >  2:7529469/1‑150 (MQ=255)
                                                                                                      gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                                               <  2:1910270/150‑1 (MQ=255)
                                                                                                      gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                                               <  1:7324704/150‑1 (MQ=255)
                                                                                                      gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                                               <  2:1109029/150‑1 (MQ=255)
                                                                                                       tAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAAACCAAACCCCAAAGGGGCGAGGGGACGGTCAGAGCACTTTTGTACTTTGTCATATGACAAAAAATGCGCCGAAGCGCCGTTAGAAAATAGTCGAATCAGTGAATTAATg                                              >  2:3743898/1‑150 (MQ=255)
                                                                                                              tCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCg                                       <  1:6805975/150‑1 (MQ=255)
                                                                                                                 cgcCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTa                                    <  2:9000105/150‑1 (MQ=255)
                                                                                                                 cgcCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTa                                    <  2:2414616/150‑1 (MQ=255)
                                                                                                                   ccAGTGCAGGCACGTTCCCCTCACCCTAACCCTCTCCCGAAAGGGTAGATGTGTCTTTCTTAGCTCTGTTGTACTTTGTTATGTGAATTATATGGCGCCGAAGCGCCTTTATAAATTAGTCGTTTCAGTGATTTACTGGTATTCGCTTAt                                  <  1:6999227/150‑1 (MQ=255)
                                                                                                                    cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                                 <  1:7797501/150‑1 (MQ=255)
                                                                                                                    cAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACCTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                                 <  1:4495722/150‑1 (MQ=255)
                                                                                                                      gTGCCGGCCTGTGCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCgg                               <  1:582982/150‑1 (MQ=255)
                                                                                                                         ccGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTAc                            <  2:9657956/150‑1 (MQ=255)
                                                                                                                                   ccccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAg                  <  1:560236/150‑1 (MQ=255)
                                                                                                                                   ccccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAg                  <  1:492575/150‑1 (MQ=255)
                                                                                                                                     ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGaa                <  2:3526709/150‑1 (MQ=255)
                                                                                                                                       tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACa              <  2:6436432/150‑1 (MQ=255)
                                                                                                                                            cTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGttc         >  1:9518401/1‑149 (MQ=255)
                                                                                                                                            cTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTa         >  1:6053321/1‑150 (MQ=255)
                                                                                                                                             tAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTAc        >  2:3679492/1‑150 (MQ=255)
                                                                                                                                              aaCCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACg       <  1:9600502/150‑1 (MQ=255)
                                                                                                                                               aCCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACgg      <  1:7976824/150‑1 (MQ=255)
                                                                                                                                                 ccTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTc    <  1:9307624/150‑1 (MQ=255)
                                                                                                                                                  ctctCCCTGAAAGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCg   <  1:3307878/150‑1 (MQ=255)
                                                                                                                                                  ctctCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCg   >  1:2319772/1‑150 (MQ=255)
                                                                                                                                                   tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGc                                      >  2:1339293/1‑114 (MQ=255)
                                                                                                                                                   tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc  >  2:8381330/1‑150 (MQ=255)
                                                                                                                                                   tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc  >  2:6876392/1‑150 (MQ=255)
                                                                                                                                                   tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc  >  1:3808557/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGC  >  NC_000913/3045921‑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: