Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,244,283 G→C 7.1% P83A (CCG→GCG)  hmgA ← homogentisate 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,244,2830GC 7.1% 402.2 / 17.2 126P83A (CCG→GCG) hmgAhomogentisate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (58/59);  new base C (6/3);  total (64/62)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 129 total aligned reads displayed.

GTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGA  >  NC_002947/5244134‑5244430
                                                                                                                                                     |                                                                                                                                                   
gtagatgCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCgg                                                                                                                                                     >  1:1092847/1‑150 (MQ=255)
       cTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGgcgc                                                                                                                                              <  2:326745/150‑1 (MQ=255)
          aCGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGcc                                                                                                                                           >  2:1795366/1‑150 (MQ=255)
          aCGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGcc                                                                                                                                           >  2:79361/1‑150 (MQ=255)
            gCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAg                                                                                                                                         >  2:957626/1‑150 (MQ=255)
             ccAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAgg                                                                                                                                        <  1:1650133/150‑1 (MQ=255)
             ccAGCCGGTTTTTCCCCGCCGGCGGTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAgg                                                                                                                                        <  1:764032/150‑1 (MQ=255)
              cAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGgc                                                                                                                                       <  1:558928/150‑1 (MQ=255)
                  cGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTc                                                                                                                                   >  2:1977057/1‑150 (MQ=255)
                   ggTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGCCTGGCGCGCCAGCCGCt                                                                                                                                    <  2:944218/148‑1 (MQ=255)
                   ggTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGCCTGGCGCGCCAGCCGCt                                                                                                                                    >  1:944218/1‑148 (MQ=255)
                     tttttCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGaa                                                                                                                                <  1:1778595/150‑1 (MQ=255)
                       tttCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGACACCCAGGGGCCCGCCGGGCGGCTGGGGCGCCAGGAGCCAGAAgg                                                                                                                              >  2:1000143/1‑149 (MQ=255)
                         tCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGt                                                                                                                            <  1:1414455/150‑1 (MQ=255)
                         tCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGt                                                                                                                            <  2:1906267/150‑1 (MQ=255)
                           gtccccGGCGTTGGGAGCCATGTGCAGCAAACTTTCGAGCAAATCGATCCGATCGGCAATAATAAGCTGCGGGATTAAGCGCAGGCTGTTGGGGTTGATGCCACGCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTgg                                                                                                                          <  2:421408/146‑1 (MQ=255)
                            gCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGCCTGGCGCGCCAGCCGCt                                                                                                                                    >  1:923566/1‑139 (MQ=255)
                            gCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGCCTGGCGCGCCAGCCGCt                                                                                                                                    <  2:923566/139‑1 (MQ=255)
                                cGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGACACCGGTTGCTTAAACTGAAAGCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGCGCGGCTGGCGCGCCAGGCGCTCGAAGCGTgggagat                                                                                                                     >  1:421408/1‑146 (MQ=255)
                                cGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGACGCGCCAGGCGCTCGAAGCGTGGGTGCa                                                                                                                     <  1:135966/150‑1 (MQ=255)
                                cGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGACGCGCCAGGCGCTCGAAGCGTGCGTGCa                                                                                                                     >  2:135966/1‑150 (MQ=255)
                                   cGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAgg                                                                                                                  <  1:79361/150‑1 (MQ=255)
                                    gTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGc                                                                                                                 >  1:847100/1‑150 (MQ=255)
                                       ggCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGACAAGCAGGAATCGGCTGCGGGCTCCACCGCAGGCGGTTGGGGTTGATGCCACCCAGTTGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGa                                                                                                              >  2:275077/1‑150 (MQ=255)
                                          aGCCAGGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAgg                                                                                                           >  2:898946/1‑150 (MQ=255)
                                           gCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAggg                                                                                                          <  1:670629/150‑1 (MQ=255)
                                             cATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGt                                                                                                                            >  1:469234/1‑130 (MQ=255)
                                             cATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGt                                                                                                                            <  2:469234/130‑1 (MQ=255)
                                             cATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCg                                                                                                        >  1:1082657/1‑150 (MQ=255)
                                              aTGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGa                                                                                                       >  1:2207331/1‑150 (MQ=255)
                                                gggCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAAt                                                                                                     >  1:882014/1‑150 (MQ=255)
                                                   cAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCg                                                                                                  <  1:1506670/150‑1 (MQ=255)
                                                   cAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCg                                                                                                  <  1:2318793/150‑1 (MQ=255)
                                                    aGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGTgg                                                                                                 >  2:752199/1‑150 (MQ=255)
                                                      ccAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTa                                                                                               <  2:1787035/150‑1 (MQ=255)
                                                        aaCCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCTGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACa                                                                                             <  2:1737893/150‑1 (MQ=255)
                                                         aCCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAg                                                                                            <  1:1264921/150‑1 (MQ=255)
                                                          ccTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGc                                                                                           <  2:1180664/150‑1 (MQ=255)
                                                           cTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGcc                                                                                          <  2:709008/150‑1 (MQ=255)
                                                            ttCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGGGGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCGGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCa                                                                                         >  1:792682/1‑150 (MQ=255)
                                                             tCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAg                                                                                        <  1:2344058/150‑1 (MQ=255)
                                                               gATGAAATCGGTCGGTTCAGCAGGAATCGGGTGGGGGCGCCAGCGCAGACGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGt                                                                                      >  2:92685/1‑150 (MQ=255)
                                                               gATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGaa                                                                                                                                >  2:1647051/1‑108 (MQ=255)
                                                               gATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGaa                                                                                                                                <  1:1647051/108‑1 (MQ=255)
                                                                 tGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTAc                                                                                    <  2:218907/150‑1 (MQ=255)
                                                                  gAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAgg                                                                                                                  >  2:2316121/1‑119 (MQ=255)
                                                                  gAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAgg                                                                                                                  <  1:2316121/119‑1 (MQ=255)
                                                                          tcggtTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCt                                                                           >  2:1076808/1‑150 (MQ=255)
                                                                              ttCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACa                                                                                             <  2:1561751/128‑1 (MQ=255)
                                                                              ttCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCt                                                                       <  2:554105/150‑1 (MQ=255)
                                                                                  gcagGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCgg                                                                   >  1:864860/1‑150 (MQ=255)
                                                                                  gcagGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCgg                                                                   >  1:1429445/1‑150 (MQ=255)
                                                                                   cagGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCCCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTTCAAGGCAGAAGGGCTAATGCGGTACAGCCAGTTACGGCGCAGCTCGCTGCggg                                                                  <  1:92685/150‑1 (MQ=255)
                                                                                    agGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGc                                                                 >  2:2067231/1‑150 (MQ=255)
                                                                                    agGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGc                                                                 >  2:1965978/1‑150 (MQ=255)
                                                                                      gAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCAGCTgg                                                                                                                                                 >  1:1458571/1‑68 (MQ=255)
                                                                                        aTCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATg                                                             <  2:1318848/150‑1 (MQ=255)
                                                                                            gCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGa                                                         >  1:828707/1‑150 (MQ=255)
                                                                                                 gggCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGcc                                                    <  1:1839115/150‑1 (MQ=255)
                                                                                                   gCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGt                                                  >  1:1456522/1‑150 (MQ=255)
                                                                                                     tCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCGAGCTGCgc                                                                   >  2:602662/1‑130 (MQ=255)
                                                                                                     tCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCGAGCTGCgc                                                                   <  1:602662/131‑2 (MQ=255)
                                                                                                       cAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGccc                                              >  2:1348481/1‑150 (MQ=255)
                                                                                                       cAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGccc                                              >  2:1511202/1‑150 (MQ=255)
                                                                                                         gcgcAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGGAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCgg                                                                   >  1:404423/1‑127 (MQ=255)
                                                                                                         gcgcAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCgg                                                                   <  2:404423/127‑1 (MQ=255)
                                                                                                           gcAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGAca                                          <  2:2194565/150‑1 (MQ=255)
                                                                                                             aGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGt                                                                                                                            >  2:1157542/1‑66 (MQ=255)
                                                                                                              ggCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGAcagca                                       <  2:279155/150‑1 (MQ=255)
                                                                                                                 ggTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGt                                                  <  1:564677/136‑1 (MQ=255)
                                                                                                                 ggTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGt                                                  >  2:564677/1‑136 (MQ=255)
                                                                                                                   ttGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCg                                  <  1:732930/150‑1 (MQ=255)
                                                                                                                       ggTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCAt                                                              >  2:162514/1‑118 (MQ=255)
                                                                                                                       ggTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCAt                                                              <  1:162514/118‑1 (MQ=255)
                                                                                                                        gTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGGCCGACAGCAGCTCGGAAtt                             >  1:36136/1‑149 (MQ=255)
                                                                                                                        gTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATa                             >  1:1202687/1‑150 (MQ=255)
                                                                                                                        gTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCCAAGGCGGTACAGCCAGGTAAAGAGCCCAACTGGGCGCACCGGGGGTAACGCCGTGCCCGACAGCAGCTCGGCATa                             >  1:541467/1‑150 (MQ=255)
                                                                                                                         ttGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGc                                                                                                                              <  2:1937627/52‑1 (MQ=255)
                                                                                                                         ttGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATAc                            <  2:1202687/150‑1 (MQ=255)
                                                                                                                            aTGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGt                                                                                                                            >  2:349893/1‑51 (MQ=255)
                                                                                                                            aTGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGt                                                                                                                            <  1:349893/51‑1 (MQ=255)
                                                                                                                             tGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGc                        <  1:352241/150‑1 (MQ=255)
                                                                                                                               ccACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCg                                                                                                  >  1:1329334/1‑74 (MQ=255)
                                                                                                                               ccACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCg                                                                                                  <  2:1329334/74‑1 (MQ=255)
                                                                                                                               ccACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCa                      >  1:1314356/1‑150 (MQ=255)
                                                                                                                                 aCCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTAcc                                                                                             <  2:1802047/77‑2 (MQ=255)
                                                                                                                                 aCCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACa                                                                                             >  1:1802047/1‑77 (MQ=255)
                                                                                                                                   ccAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACa                                                                                             >  2:970478/1‑75 (MQ=255)
                                                                                                                                   ccAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACa                                                                                             <  1:970478/75‑1 (MQ=255)
                                                                                                                                   ccAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACg                  >  1:2185459/1‑150 (MQ=255)
                                                                                                                                     aGTGGCCCGCCGAGCGGCTGGCGCGCCAGGAGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACggg                <  1:1795366/150‑1 (MQ=255)
                                                                                                                                       tGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGc              <  1:671536/150‑1 (MQ=255)
                                                                                                                                        ggCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGcc             <  1:1965978/150‑1 (MQ=255)
                                                                                                                                           gcttcgAGCGCCTGGCGCGCCAGCCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGTCCGAc                                           <  2:1070778/113‑1 (MQ=39)
                                                                                                                                           gcttcgAGCGCCTGGCGCGCCAGCCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGTCCGAc                                           >  1:1070778/5‑117 (MQ=39)
                                                                                                                                             ttcgAGCGCCTGGCGCGCCAGCCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTg        >  1:687385/3‑150 (MQ=255)
                                                                                                                                             gccgAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTg        >  2:8524/1‑150 (MQ=255)
                                                                                                                                              ccgAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGc       >  2:1511165/1‑150 (MQ=255)
                                                                                                                                              ccgAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGc       >  1:2124763/1‑150 (MQ=255)
                                                                                                                                                gAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCgg     <  1:1000143/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGA  >  NC_002947/5244134‑5244430

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