Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,919,901 G→C 10.5% G282G (GGC→GGG) ‡ gudX ← glucarate dehydratase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,919,9010GC10.5% 263.7 / 19.0 86G282G (GGC→GGG) ‡gudXglucarate dehydratase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base C (5/4);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

CCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAA  >  NC_000913/2919752‑2920049
                                                                                                                                                     |                                                                                                                                                    
ccccAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAg                                                                                                                                                      <  1:574969/150‑1 (MQ=255)
ccccAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAg                                                                                                                                                      >  2:589097/1‑150 (MQ=255)
      gcgtcgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGGCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTc                                                                                                                                                <  1:820612/149‑1 (MQ=255)
       cgtcgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCg                                                                                                                                               <  2:195810/150‑1 (MQ=255)
         tcgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGcc                                                                                                                                             <  1:1125038/150‑1 (MQ=255)
           gCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCg                                                                                                                                           <  1:1417009/150‑1 (MQ=255)
           gCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCg                                                                                                                                           <  2:487312/150‑1 (MQ=255)
            cAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGt                                                                                                                                          <  2:642011/150‑1 (MQ=255)
              aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCg                                                                                                                                        <  2:1278581/150‑1 (MQ=255)
              aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCg                                                                                                                                        <  2:144656/150‑1 (MQ=255)
              aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCg                                                                                                                                        <  1:636935/150‑1 (MQ=255)
              aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCg                                                                                                                                        <  2:27191/150‑1 (MQ=255)
              aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCAACCGCATGGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCGACCGGCTTGCCCGTCg                                                                                                                                        <  2:1417514/150‑1 (MQ=255)
               aGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCGACCGGCTTGCCCGTCgc                                                                                                                                       <  1:838071/150‑2 (MQ=255)
               aGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGa                                                                                                                                       <  1:629019/150‑1 (MQ=255)
                gCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCGACCGGCTTGCCCGTCgcg                                                                                                                                      <  1:95948/150‑3 (MQ=255)
                      ccACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCg                                                                                                                                <  1:981149/150‑1 (MQ=255)
                      ccACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCg                                                                                                                                <  1:622644/150‑1 (MQ=255)
                       cacaCGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGGGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGc                                                                                                                               <  2:1598935/150‑1 (MQ=255)
                                 gCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCa                                                                                                                     >  1:743700/1‑150 (MQ=255)
                                          aGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGgag                                                                                                            <  2:25133/150‑1 (MQ=255)
                                            cGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGaa                                                                                                          <  2:1232018/150‑1 (MQ=255)
                                             gTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAg                                                                                                         >  1:113973/1‑150 (MQ=255)
                                               ccAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGcc                                                                                                       >  1:697943/1‑150 (MQ=255)
                                               ccAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGcc                                                                                                       >  1:871674/1‑150 (MQ=255)
                                                   aaGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGt                                                                                                   <  2:1391959/150‑1 (MQ=255)
                                                       gCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTGGGCCGTCGAAATTCCGCCATCAGTTCACGTCCGGAGAAGCCCTGTTCTg                                                                                               <  2:698591/150‑1 (MQ=255)
                                                        cGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTgc                                                                                              <  1:1121468/150‑1 (MQ=255)
                                                             cGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGcgccgc                                                                                         <  1:910285/150‑1 (MQ=255)
                                                              ggCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCa                                                                                        >  1:153/1‑150 (MQ=255)
                                                               gCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCAt                                                                                       >  1:1535975/1‑150 (MQ=255)
                                                                 aaGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATgg                                                                                     <  2:1167882/150‑1 (MQ=255)
                                                                  aGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGa                                                                                    >  2:805478/1‑150 (MQ=255)
                                                                     ggAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCt                                                                                 >  1:1492337/1‑150 (MQ=255)
                                                                      gAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTACGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCtt                                                                                >  2:508112/1‑150 (MQ=255)
                                                                        atatCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCg                                                                              <  1:659153/150‑1 (MQ=255)
                                                                          atCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGc                                                                            <  1:169052/150‑1 (MQ=255)
                                                                             tACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATa                                                                         <  2:1476640/150‑1 (MQ=255)
                                                                             tACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATa                                                                         >  1:596126/1‑150 (MQ=255)
                                                                               ccGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAgg                                                                       <  2:1206613/150‑1 (MQ=255)
                                                                                 gCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTa                                                                     >  1:994389/1‑150 (MQ=255)
                                                                                 gCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTa                                                                     <  1:1514426/150‑1 (MQ=255)
                                                                                      gAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTc                                                                                                                 >  2:684307/1‑101 (MQ=255)
                                                                                      gAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTc                                                                                                                 <  1:684307/101‑1 (MQ=255)
                                                                                      gAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAAc                                                                >  1:381343/1‑150 (MQ=255)
                                                                                        gCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAAcat                                                              >  2:704367/1‑150 (MQ=255)
                                                                                        gCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCAGCATAGGTAAGAAcat                                                              >  2:602009/1‑150 (MQ=255)
                                                                                             aCCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCTCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTc                                                         >  2:1143083/1‑150 (MQ=255)
                                                                                             aCCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTc                                                         <  2:1492337/150‑1 (MQ=255)
                                                                                                gCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGc                                                      <  1:526265/150‑1 (MQ=255)
                                                                                                  aTGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGccc                                                    >  1:1276337/1‑150 (MQ=255)
                                                                                                      cccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGGAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTg                                                >  1:1069982/1‑150 (MQ=255)
                                                                                                      cccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTg                                                >  2:35120/1‑150 (MQ=255)
                                                                                                      cccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTg                                                >  1:1457311/1‑150 (MQ=255)
                                                                                                      cccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTg                                                >  2:1144946/1‑150 (MQ=255)
                                                                                                       ccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGc                                               <  2:601632/150‑1 (MQ=255)
                                                                                                         aTTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCAt                                             <  2:113973/150‑1 (MQ=255)
                                                                                                                gcCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGaaa                                      >  2:53413/1‑150 (MQ=255)
                                                                                                                 ccAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAAt                                     >  1:915064/1‑150 (MQ=255)
                                                                                                                  cAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATg                                    <  2:697943/150‑1 (MQ=255)
                                                                                                                    gTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGc                                  <  1:1419733/150‑1 (MQ=255)
                                                                                                                            cGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCaa                          >  1:1328948/1‑150 (MQ=255)
                                                                                                                            cGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCaa                          <  2:1535975/150‑1 (MQ=255)
                                                                                                                              aTCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGTCGCATTTATATTATTCTTAGGTAAGATAATCATTCAGGCTTTTTCATTAAGAAATGTCTtgatgatgc                        <  2:1298945/150‑6 (MQ=255)
                                                                                                                                 aTGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAagcagc                     >  1:686010/1‑150 (MQ=255)
                                                                                                                                  tGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCttt                                                 >  1:663888/1‑121 (MQ=255)
                                                                                                                                  tGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCttt                                                 <  2:663888/121‑1 (MQ=255)
                                                                                                                                  tGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGcc                    >  1:1458469/1‑150 (MQ=255)
                                                                                                                                       gTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCa               >  1:156202/1‑150 (MQ=255)
                                                                                                                                        tttcgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCAc              >  2:1450376/4‑150 (MQ=255)
                                                                                                                                         ttcgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCAcc             >  1:5454/3‑150 (MQ=255)
                                                                                                                                         ttcgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCAcc             >  1:810939/3‑150 (MQ=255)
                                                                                                                                         ttcgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCAcc             >  2:1408161/3‑150 (MQ=255)
                                                                                                                                         cgcgACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCAcc             >  2:1155101/1‑150 (MQ=255)
                                                                                                                                          gcgACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCg            >  1:535801/1‑150 (MQ=255)
                                                                                                                                           cgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGt           >  1:527897/1‑150 (MQ=255)
                                                                                                                                           cgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGt           <  2:527897/150‑1 (MQ=255)
                                                                                                                                           cgACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGt           >  2:647778/1‑150 (MQ=255)
                                                                                                                                              cGGGCAAGCCGGTCGCCCGTCGGAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTgg        <  2:1069982/150‑1 (MQ=255)
                                                                                                                                              cGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTg         >  2:364347/1‑149 (MQ=255)
                                                                                                                                              cGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTg         <  1:364347/149‑1 (MQ=255)
                                                                                                                                               gggCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTggg       <  2:1427560/150‑1 (MQ=255)
                                                                                                                                                 gCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCTCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTgatcc     <  1:1143083/150‑5 (MQ=255)
                                                                                                                                                 gCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGAt     >  2:1467566/1‑150 (MQ=255)
                                                                                                                                                  cAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGc                                                                            >  2:405717/1‑78 (MQ=255)
                                                                                                                                                  cAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGc                                                                            <  1:405717/78‑1 (MQ=255)
                                                                                                                                                    aGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCaa  >  1:349210/1‑150 (MQ=255)
                                                                                                                                                    aGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCaa  <  2:743700/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAA  >  NC_000913/2919752‑2920049

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