Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,050,713 A→T 5.4% intergenic (‑46/+402) gcvT ← / ← ubiI aminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex/2‑octaprenylphenol hydroxylase, FAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,050,7130AT 5.4% 386.1 / 8.5 130intergenic (‑46/+402)gcvT/ubiIaminomethyltransferase, tetrahydrofolate‑dependent, subunit (T protein) of glycine cleavage complex/2‑octaprenylphenol hydroxylase, FAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (63/60);  new base T (5/2);  total (68/62)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 132 total aligned reads displayed.

CGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATC  >  NC_000913/3050565‑3050860
                                                                                                                                                    |                                                                                                                                                    
cgATTTGCGAACCGTAATGCAGCGGCATCATCCAGACGTGGAAATCCACCATGCGAGCGCCGCAAAGCGAGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCCTCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGAGaaa                                                                                                                                                     >  1:292114/1‑150 (MQ=255)
 gATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGAAAACTTTTTCATGGTGAAAt                                                                                                                                                    >  1:1252304/1‑150 (MQ=255)
  aTTTGCGATCCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCCTCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAAtt                                                                                                                                                   >  2:1185426/1‑150 (MQ=255)
        gAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCAATCTTGTCCTCATTGAATAAGCTGCGCTGACAACTTTTTCATGGTGAAATTATcacc                                                                                                                                             >  1:1966396/1‑150 (MQ=255)
            cGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGa                                                                                                                                         >  1:1641065/1‑150 (MQ=255)
               aaTGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGATAATTTCACCATGaaaa                                                                                                                                      <  1:1732868/150‑2 (MQ=255)
               aaTGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAAc                                                                                                                                      <  1:1388076/150‑1 (MQ=255)
                  gCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCACAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTACTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATGACAACGATAcccc                                                                                                                                   >  1:94163/1‑149 (MQ=255)
                  gCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAACCCCa                                                                                                                                   >  2:1952880/1‑150 (MQ=255)
                   cAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAg                                                                                                                                  <  2:2005561/150‑1 (MQ=255)
                    aGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCCCGCAACCCAGc                                                                                                                                 >  2:276516/1‑150 (MQ=255)
                        gCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCTGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCg                                                                                                                             <  1:1096919/150‑1 (MQ=255)
                        gCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCCGCATCg                                                                                                                             >  2:1695991/1‑150 (MQ=255)
                         catcatCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCgg                                                                                                                            >  1:1473679/1‑150 (MQ=255)
                           caagccATCCGTGGATATACACCATGCGAGCGCCGCAAATCGGGGGTTGTGCGTACAAAGGAGTCTGTTGTGACATATTGTCCTCATTGAATAAGCGTGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAAACCAGCATCGGAg                                                                                                                          <  2:825472/146‑1 (MQ=255)
                             atcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGcc                                                                                                                        <  1:910223/150‑1 (MQ=255)
                               ccAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACTAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCAc                                                                                                                      <  2:701137/150‑1 (MQ=255)
                                   ccGTGGAAATCCACCATGCGAGCGCCGCAAAGCGAGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTCTCCTCATAGAATAAGCGGGGCTGACAACTTTTTCATGGTCAAATTATCACCCCGAATCCCAGCATCGGAGCCACTccc                                                                                                                  >  1:33975/1‑150 (MQ=255)
                                    cGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCCGCATCGGAGCCACTCCCg                                                                                                                 >  2:639050/1‑150 (MQ=255)
                                     gTGGAACTCCACCATGCGAGCGCCGCAAAGCGGGTGTTGGTCGTACAAAGGAGTCTGTTGGGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCgg                                                                                                                <  2:1231228/150‑1 (MQ=255)
                                     gTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCgg                                                                                                                >  2:906334/1‑150 (MQ=255)
                                      tGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGCAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATAATCACCACGAAACCCAGCATCGGAGCCACACCCGGt                                                                                                               >  1:1551991/1‑150 (MQ=255)
                                       ggAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTc                                                                                                              <  1:1931510/150‑1 (MQ=255)
                                        gAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTcc                                                                                                             <  1:1388190/150‑1 (MQ=255)
                                         aaaTCCACCATGCGAGCGCCGCAAAGCGGGGGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTccc                                                                                                            <  2:977602/150‑1 (MQ=255)
                                         aaaTCCACCATGCGAGCGCCGCAAAACGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTccc                                                                                                            <  2:139099/150‑1 (MQ=255)
                                          aaTCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCCCTCCCGGTcccc                                                                                                           >  1:633489/1‑150 (MQ=255)
                                           aTCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCa                                                                                                          <  1:1642005/150‑1 (MQ=255)
                                           aTCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCa                                                                                                          >  1:529432/1‑150 (MQ=255)
                                             ccaccaTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGTTGCTGACATCTTTTTCATGGTGAAATTATCACCACGAAACGCAGCATCGGAGCCACTCCCGGTCCCCAAc                                                                                                        <  2:1162172/150‑1 (MQ=255)
                                             ccaccaTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAAc                                                                                                        <  2:1185103/150‑1 (MQ=255)
                                              caccaTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACg                                                                                                       >  1:747406/1‑150 (MQ=255)
                                                  aTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTAGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAAt                                                                                                   <  2:1102954/150‑1 (MQ=255)
                                                   tGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATc                                                                                                  >  2:1153374/1‑150 (MQ=255)
                                                    gCGACCGCCGCAAAGCGGGTGTTGTTCGTACAAAGGAGGCTGTTGTGCCATCTTGGCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCg                                                                                                 <  1:1142692/150‑1 (MQ=255)
                                                     cGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGt                                                                                                >  2:1631438/1‑150 (MQ=255)
                                                         ctccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTTTTGTGCCATCTTGTCCTCATTGAATAAGCGGGTCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTtctc                                                                                            <  1:1239168/148‑1 (MQ=255)
                                                         cgccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTtctc                                                                                            <  1:1232858/150‑1 (MQ=255)
                                                          gccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTctct                                                                                           >  2:103251/1‑150 (MQ=255)
                                                           ccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCtt                                                                                          >  1:1814598/1‑150 (MQ=255)
                                                            cgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCACCGCAATCGTTCTCttt                                                                                         >  2:878446/1‑150 (MQ=255)
                                                             gcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGAGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCtttt                                                                                        >  1:597038/1‑150 (MQ=255)
                                                              ctaattgtgtgtTTTTCGTACAATGGAGTCTGTTGTGCCATGTTTTCCTCATTGAATAAGAGGGGCTGTCTACTTTTTCATGGTGAAATTATCACCAGGAATCCCAGCATCGGATCCACTCCAGGTCCCCAACGCAATCGTTCTCTTTTg                                                                                       <  1:2009130/144‑1 (MQ=255)
                                                                aaGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGcc                                                                                     <  2:1146763/150‑1 (MQ=255)
                                                                aaGCGTGTGTTGGTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGcc                                                                                     <  1:978394/150‑1 (MQ=255)
                                                                    gtgtgtTGTTCGTACTTAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAATTTATCACCACGAAACCCAGGATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGaa                                                                                 <  2:1619499/150‑1 (MQ=255)
                                                                     tgtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGTGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGTGCCTCTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAAc                                                                                <  1:1040964/150‑1 (MQ=255)
                                                                     tgtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAAc                                                                                <  1:85193/150‑1 (MQ=255)
                                                                     tgtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAAc                                                                                <  1:1972165/150‑1 (MQ=255)
                                                                      gtgtTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACt                                                                               <  2:322443/150‑1 (MQ=255)
                                                                        gttgttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTa                                                                             >  2:1586562/1‑150 (MQ=255)
                                                                           gttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCCCCACGAAACACAGCATCGGCGCCACTCACGGTCCCCAACTCACTCGTTCTCTTTTGCATGAACTTcccc                                                                          >  1:1223161/1‑149 (MQ=255)
                                                                            ttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTaccac                                                                         >  1:1149890/1‑150 (MQ=255)
                                                                            ttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTaccac                                                                         <  2:758613/150‑1 (MQ=255)
                                                                             tCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCTTCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTATGCCTGAACTTaccacc                                                                        >  2:488916/1‑150 (MQ=255)
                                                                              cGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCg                                                                       >  1:1849901/1‑150 (MQ=255)
                                                                               gTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGa                                                                      >  2:794543/1‑150 (MQ=255)
                                                                                tACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGaa                                                                     <  1:1141313/150‑1 (MQ=255)
                                                                                 aCAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCATGAACTTACCACCGaaa                                                                    >  2:1481540/1‑150 (MQ=255)
                                                                                  cAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAAc                                                                   >  1:667855/1‑150 (MQ=255)
                                                                                             gttgtGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAc                                                        <  1:775626/150‑1 (MQ=255)
                                                                                              ttgtGCCATCTAGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAAATATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAcc                                                       <  2:1850488/150‑1 (MQ=255)
                                                                                               tgtgCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCa                                                      <  1:1631438/150‑1 (MQ=255)
                                                                                                    cATCTTGTCCTCATTGAATAAGCTGGGCTGACAACTTTATCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGGCCCCAACGCACTCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAgg                                                 <  2:154545/150‑1 (MQ=255)
                                                                                                      tCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTa                                               >  2:1287598/1‑150 (MQ=255)
                                                                                                           tCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAAtt                                          >  2:570099/1‑150 (MQ=255)
                                                                                                           tCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAAtt                                          <  2:1814598/150‑1 (MQ=255)
                                                                                                             cTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGa                                        >  2:1246682/1‑150 (MQ=255)
                                                                                                             cTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTATCTTTTTCCTGAACTTACCACCGAAACAGACTGTTAACCAAAATGTAAAATTGa                                        >  2:1848597/1‑150 (MQ=255)
                                                                                                             cTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATAACCACGAAACCCAGCATCGGAGCCACTCACGGTCCCCAACGCAATCGTTCTATTTTGCCTGAACTTACCACCGAAACAGACAGTTAACCATAAGGTAAAATTGa                                        >  1:1861408/1‑150 (MQ=255)
                                                                                                               cATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatc                                      <  2:1473679/150‑1 (MQ=255)
                                                                                                                aTTGAATAAGCGGGGCTGACAACTTTTTCATTGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatca                                     <  2:493204/150‑1 (MQ=255)
                                                                                                                aTTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATAACCACGAAACCCAGCATCGGAGCCACTCACGGTCCCCAACGCAAACGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTATACCATAAGGTAAAATTGatca                                     >  1:1226868/1‑150 (MQ=255)
                                                                                                                   gAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACCGACTGTTAACCATAAGGTAAAATTGAtcatca                                 <  2:1037509/150‑1 (MQ=255)
                                                                                                                    aaTAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATcac                                <  2:1474383/150‑1 (MQ=255)
                                                                                                                     aTAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATcaca                               >  1:957245/1‑150 (MQ=255)
                                                                                                                      tAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACAt                              >  2:1056242/1‑150 (MQ=255)
                                                                                                                      tAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGGTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACAt                              >  1:964726/1‑150 (MQ=255)
                                                                                                                        aGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTa                            >  2:1544481/1‑150 (MQ=255)
                                                                                                                         gCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAg                           <  1:909198/150‑1 (MQ=255)
                                                                                                                          cGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGc                          >  1:458728/1‑150 (MQ=255)
                                                                                                                           ggggCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCt                         >  1:1721847/1‑150 (MQ=255)
                                                                                                                               cTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATAGTTCTCTTTTGCCTGAACTAACCACCCAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTAtt                     >  2:1554178/1‑149 (MQ=255)
                                                                                                                                tGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATgg                    <  2:1727493/150‑1 (MQ=255)
                                                                                                                                tGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATgg                    <  2:1136895/150‑1 (MQ=255)
                                                                                                                                tGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCAAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATgg                    <  2:1746335/150‑1 (MQ=255)
                                                                                                                                tGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGACACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATgg                    <  1:386999/150‑1 (MQ=255)
                                                                                                                                  tgctgggtttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGtt                  >  2:491025/8‑150 (MQ=255)
                                                                                                                                  aCAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTGATGGtt                  >  2:1380582/1‑150 (MQ=255)
                                                                                                                                   cAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCTGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTa                 <  1:169920/150‑1 (MQ=255)
                                                                                                                                      cTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCTTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGTTTaaaa              >  2:280340/1‑150 (MQ=255)
                                                                                                                                      cTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaa              >  1:1367180/1‑150 (MQ=255)
                                                                                                                                       tttttCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGCTAATTTGATGATCAATTTTACATTATGGTTaacag             <  1:1449273/150‑4 (MQ=255)
                                                                                                                                        gtttCGTGGTGATAATTTCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACGACCGAAACAGACGGTTAACAATAAGGTAAAATTGATCATCACTTTAGCTTATGGTTaaaaaa            >  1:1024830/2‑150 (MQ=255)
                                                                                                                                          ttCATTGTGAAAATATCAAAACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTTAACTTACCACCGAAACCGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGGTAAAAAATg          <  2:1718548/150‑1 (MQ=255)
                                                                                                                                           tCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTTCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGc         <  2:1223161/150‑1 (MQ=255)
                                                                                                                                               ggTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaaa     >  2:1912831/1‑150 (MQ=255)
                                                                                                                                               ggTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaaa     >  1:554294/1‑150 (MQ=255)
                                                                                                                                                gTGAAATTATCACCACGAAACTCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaaaa    <  1:1188582/150‑1 (MQ=255)
                                                                                                                                                  gAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATc  <  1:1056242/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATC  >  NC_000913/3050565‑3050860

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: