Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,919,905 G→C 10.4% T281S (ACC→AGC)  gudX ← glucarate dehydratase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,919,9050GC10.4% 266.2 / 18.9 87T281S (ACC→AGC) gudXglucarate dehydratase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (35/43);  new base C (5/4);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGG  >  NC_000913/2919758‑2920052
                                                                                                                                                   |                                                                                                                                                   
gcgtcgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGGCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTc                                                                                                                                                   <  1:820612/149‑1 (MQ=255)
 cgtcgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCg                                                                                                                                                  <  2:195810/150‑1 (MQ=255)
   tcgCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGcc                                                                                                                                                <  1:1125038/150‑1 (MQ=255)
     gCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCg                                                                                                                                              <  2:487312/150‑1 (MQ=255)
     gCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCg                                                                                                                                              <  1:1417009/150‑1 (MQ=255)
      cAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGt                                                                                                                                             <  2:642011/150‑1 (MQ=255)
        aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCg                                                                                                                                           <  2:1278581/150‑1 (MQ=255)
        aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCg                                                                                                                                           <  1:636935/150‑1 (MQ=255)
        aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCg                                                                                                                                           <  2:27191/150‑1 (MQ=255)
        aaGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCg                                                                                                                                           <  2:144656/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:622644/150‑1 (MQ=255)
                ccACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCg                                                                                                                                   <  1:981149/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                                                                            >  1:596126/1‑150 (MQ=255)
                                                                       tACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATa                                                                            <  2:1476640/150‑1 (MQ=255)
                                                                         ccGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAgg                                                                          <  2:1206613/150‑1 (MQ=255)
                                                                           gCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTa                                                                        <  1:1514426/150‑1 (MQ=255)
                                                                           gCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTa                                                                        >  1:994389/1‑150 (MQ=255)
                                                                                gAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTc                                                                                                                    <  1:684307/101‑1 (MQ=255)
                                                                                gAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTc                                                                                                                    >  2:684307/1‑101 (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:1144946/1‑150 (MQ=255)
                                                                                                cccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTg                                                   >  2:35120/1‑150 (MQ=255)
                                                                                                cccATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTg                                                   >  1:1457311/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                             <  2:1535975/150‑1 (MQ=255)
                                                                                                                      cGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCaa                             >  1:1328948/1‑150 (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                >  2:1408161/3‑150 (MQ=255)
                                                                                                                                   ttcgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCAcc                >  1:810939/3‑150 (MQ=255)
                                                                                                                                   ttcgACGGGCGACCGGCTTGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCAcc                >  1:5454/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                                                                               <  1:405717/78‑1 (MQ=255)
                                                                                                                                            cAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGc                                                                               >  2:405717/1‑78 (MQ=255)
                                                                                                                                              aGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCaa     >  1:349210/1‑150 (MQ=255)
                                                                                                                                              aGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCaa     <  2:743700/150‑1 (MQ=255)
                                                                                                                                                 cGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTTAGAACATCATTGAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACgg  <  2:1100452/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAAGTGCGGATCGGCAAGTGGAATATCTACCGCATTGAGCATCACCGCATGACCCATTTCGCGCCAGTTGGTGGCGATCATGTTAGTCGCGACGGGCAAGCCGGTCGCCCGTCGAAATTCCGCCATCACTTCACGTCCGGAGAAGCCCTGTTCTGCGCCGCATGGATCTTCGGCATAGGTAAGAACATCATTCAGCCCTTTGCATAAAGAAATGGCTTCATCAAGCAGCCATGCACCGTTGGGATCAACGG  >  NC_000913/2919758‑2920052

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: