Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 836,592 T→G 16.5% V81G (GTC→GGC)  ybiC → putative dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913836,5920TG16.5% 160.6 / 2.9 65V81G (GTC→GGC) ybiCputative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (12/39);  new base G (12/2);  total (24/41)
Fisher's exact test for biased strand distribution p-value = 3.81e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCATCGCTTTGATGCTCAGACGCTGCACAGTTTTATTCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGC  >  NC_000913/836360‑836824
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           
ggTCATCGCTTTGATGCTCAGACGCTGCACAGTTTTATTCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGCGGCGGGGGGGGCGGTCAcag                                                                                                                                                                                                                                        >  2:956125/1‑236 (MQ=255)
 gTCATCGCTTTGATGCTCAGACGCTGCACAGTTTTATTCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGcc                                                                                                                                                                                                                                       >  1:396679/1‑237 (MQ=255)
   cATCGCTTTGATGCTCAGACGCTGCACAGTTTTATTCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGa                                                                                                                                                                                                                                    <  1:1225639/239‑1 (MQ=255)
         tttGATGCTCAGACGCTGCACAGTTTTATTCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACGCCCGCTGGCGa                                                                                                                                                                                                                             >  1:851714/1‑239 (MQ=255)
                                 tttATTCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTcgcg                                                                                                                                                                                                     <  2:142233/237‑1 (MQ=255)
                                    agtCAGTCTGTATTTCGTCAGATGGGTAGCGAGTAACAAGAAGCTAAATTAGTTGCCGATCATTTAATCGCGGCAAGCCAGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCGCGATGGCGATCGCGCATTTGGTCAGGTCGCGGcaca                                                                                                                                                                                                <  2:49950/237‑1 (MQ=255)
                                    aTTCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGcaca                                                                                                                                                                                                <  2:738949/239‑1 (MQ=255)
                                    aTTCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGcaca                                                                                                                                                                                                <  2:505696/239‑1 (MQ=255)
                                     ttCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGGCGCGGCACAt                                                                                                                                                                                               >  1:1149796/1‑239 (MQ=255)
                                      tCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATg                                                                                                                                                                                              <  2:1149796/239‑1 (MQ=255)
                                          gCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCCCCCCCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGc                                                                                                                                                                                          >  1:1199808/1‑239 (MQ=255)
                                          gCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCCCGGCACATGAAGc                                                                                                                                                                                          >  1:1038776/1‑239 (MQ=255)
                                          gCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCACTTGCGATCGCGCTTTTGGTACAGTCGCGGGACATGaatc                                                                                                                                                                                          >  1:1058052/1‑237 (MQ=255)
                                             gTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAt                                                                                                                                                                                       <  2:463047/239‑1 (MQ=255)
                                              tATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCCCATGAAGCGATg                                                                                                                                                                                      >  1:6982/1‑239 (MQ=255)
                                               aTTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATg                                                                                                                                                                                      <  2:90531/238‑1 (MQ=255)
                                                   cGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGGGGTCACGCTCGCTGGCGAACCCGCCTTTTGGCAGGGCGCGGCACATGACGCGGTGgcgcc                                                                                                                                                                                 >  2:484789/1‑238 (MQ=255)
                                                        gATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACGCTCGATTGGGATCTCGCATTTGGGCAGGTCGCGGCCCACGAAACGATGGGGCCGGGGa                                                                                                                                                                            >  1:1093936/1‑239 (MQ=255)
                                                        gATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGGGGGCCCGCCTTTGGTGAGGGGGCGGGCCATAAAGCGAGCGCGCTggggg                                                                                                                                                                            >  2:499975/1‑238 (MQ=255)
                                                        gATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTggggg                                                                                                                                                                            >  1:643172/1‑238 (MQ=255)
                                                        gATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCCCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCCCCCCCCCTGGCGAGCGCGCACTTGTGCACGGCGGGGCCCCTGGAGCGGTGGGGCCggggg                                                                                                                                                                            >  2:74212/1‑238 (MQ=255)
                                                             gTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCCCGCTCGGCGGCGGCACCCGCCTTGGGCCGGTCGCGGGGCCATGAGGCGCTGGCGCGGGGGATTgg                                                                                                                                                                        >  2:877542/1‑238 (MQ=255)
                                                                  gAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAg                                                                                                                                                                   <  2:551010/238‑1 (MQ=255)
                                                                   aGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCTCGATGGCGAACGCGCATTTTGGTAGGTCGCGGCACATGAAACGGTGGGGCTGGGGGTTggggaaag                                                                                                                                                                  >  1:850897/1‑231 (MQ=255)
                                                                   aGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCGATCGCGCCTTTTGGCAGGGCGCGGGACATGAAGCGAGGGGGCCGGGGGTTGAGaaaagg                                                                                                                                                                 >  1:699259/1‑236 (MQ=255)
                                                                       actggagtcGAAAATAGTTGCCGATCTTTTAATCGCGGCAATCCTGTCAGCTCATGATTCACATGGTATTCGCATGATCCCAATCTATGTACAGTCCTCGAGTCAGGGGCACCTGGAAATCAACCTTCATGCGAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGGTCGATGGCGATCGCGCATTTGGTCAGGGCGCGGCATATGAAGCAATGGCGCTGGGGATTGAGAAAGCGCAt                                                                                                                                                              <  2:1058052/230‑1 (MQ=255)
                                                                             agCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAc                                                                                                                                                        <  2:553288/238‑1 (MQ=255)
                                                                             agCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg                                                                                                                                                       <  1:418498/239‑1 (MQ=255)
                                                                                                       tCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAAc                                                                                                                             <  1:606073/239‑1 (MQ=255)
                                                                                                        cgcgGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACt                                                                                                                            <  2:1038776/239‑1 (MQ=255)
                                                                                                                          gCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGcc                                                                                                             <  1:856714/236‑1 (MQ=255)
                                                                                                                                 tCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTa                                                                                                   <  2:305889/239‑1 (MQ=255)
                                                                                                                                   acaTGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACt                                                                                                 <  1:499975/239‑1 (MQ=255)
                                                                                                                                   acaTGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACt                                                                                                 <  1:415496/239‑1 (MQ=255)
                                                                                                                                   acaTGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGAAGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACt                                                                                                 <  1:775059/239‑1 (MQ=255)
                                                                                                                                        gTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGt                                                                                                                                                                                                               <  1:958157/125‑1 (MQ=255)
                                                                                                                                         tATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAAc                                                                                                                             <  1:206197/205‑1 (MQ=255)
                                                                                                                                         tATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGc                                                                                                                         <  1:982297/209‑1 (MQ=255)
                                                                                                                                         tATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCg                                                                                            <  2:643172/238‑1 (MQ=255)
                                                                                                                                         tATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCg                                                                                            <  1:854607/238‑1 (MQ=255)
                                                                                                                                           ttGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAg                                                                                         <  1:561357/239‑1 (MQ=255)
                                                                                                                                                   aTCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCCCGATGGCGATCGCGCATTTGGGCCGGGGGCGGCACCTGAAGCGAGGGGGCCGGGGGTTTAGAAAGGGCCTCAGCACGGGATTTCCGCCGGGGCGCGACATAACTCCCACCACAACGGCCGTATCGGGTGCTTGGGGGGGCAGtgtgg                                                                                 >  2:58087/1‑238 (MQ=255)
                                                                                                                                                   aTCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCCAGCTCGATGGCGATCGCGCCTTTGGTCAGGGCGCGGGCCATGAAGAGGTGGCGCTGGGGATTGGGAAAGCGCATCAGCACGGTATTTGCGCCGTGGGGCTCCATTACACCCATCCTACAGGCCGTATCGGTTACTGGGGGGAGGAGggtgt                                                                                 >  2:764108/1‑238 (MQ=255)
                                                                                                                                                           cTAGGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGGCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCGACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGa                                                                                          <  1:534212/222‑1 (MQ=255)
                                                                                                                                                                     tCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACCCCCCATGGCGATCCCGCCTTTTGTCCGGTCGCGGCACCTGAAACCATGGGGGCGGGGGTTTAGGAAGGGCAACCACCCGGGGGTTGCGCCCGGGCGCCCCCTAACTCGCATCCTCTCGGCCGGATCGGGTTCTGGGGGGGGCCgtggggcgctgtcttttttacac                                                                >  2:534212/1‑218 (MQ=255)
                                                                                                                                                                     tCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGGGGGCACGCCCGATGGCGATCGCGCCTTTGGGCCGGGCGCGGGACATGAAGCGATGGCGCTGGGGGGTGAGAAAGCGCATCAGCACGGTATTGGCGCCGGGGCGCGACAATACTCGCATCCTATCGGGCGGGTCCGGTACTGGGGGGGGCCGGGGGCAGGGGCggggggtgtg                                                                 >  2:559497/1‑231 (MQ=255)
                                                                                                                                                                     tCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGtctc                                                               <  1:722284/239‑1 (MQ=255)
                                                                                                                                                                     tCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCCCGCCCCGAGGCGAACGCGCCTTTTGGTCGGGCGCGGGCCATGAAGCGGAGGGGCCGGGGGGTGGGAAAACGGCTCAGCACGGGGTTGGCGCCCGGGGGGTTCATAAACCGCCTCCTATCGGCCGGATCGGTTTCTGGGGGGGGCCGCGGGGAGCGggggggggtggg                                                                 >  2:669967/1‑231 (MQ=255)
                                                                                                                                                                            gTCAGGTGCACCTGCAAATTAACCATCATGCCAAACCCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCa                                                         <  2:371621/238‑1 (MQ=255)
                                                                                                                                                                            gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCa                                                         <  1:764108/238‑1 (MQ=255)
                                                                                                                                                                            gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCa                                                         <  1:484789/238‑1 (MQ=255)
                                                                                                                                                                               aGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGACGGCGAGCGCGCATTTTGCCAGGGCGCGGCACATGAAGCGAGGGCGCTGGGGATTGGGAAAGCGCCTCCGCACGGTAGTGACGCCGGGGCGCTGCCAAACTCCCATCATATCCGCCGCATCGGTTACTTGGGGGAGCCGGGGGGAGCGGCGGGGGTGGTCTCTATcacccat                                                     >  2:37490/1‑233 (MQ=255)
                                                                                                                                                                               aGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCAACCCCGATGGCGATCCCGCCTTTGGGCAGGGCCGGGCACATGAAACGATGGCGCGGGGGGGTGAGAAAACGCCGCAGCACCGGATTTCCGCCGGGGCGCTACATAAACCGCATCCAATCGGGCGGATAGGTTACGTGGCGGGGCAGAGTGGAGCGGCGGGGGTTTtctctctctcacct                                                     >  2:1115487/1‑233 (MQ=255)
                                                                                                                                                                                 gggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCCAACCCGCCTTTTGTCTGGGCCCGGCCCCTTAAACGGTGGGGCCGGGGGGTGAGAAAGGGCACTAGCCCCGGTTTGCCCCCGTGGCGCCACATAAATCGCAACCTCTTCGCCCGATCCGTTTCTGGGCGGGGGCGTGGGGCGCGGGGGGGGGTGtctttttttcctttct                                                   >  2:110097/1‑226 (MQ=255)
                                                                                                                                                                                           aaaTTAACCATCATGCCAAAACCGTTAAAGGGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTGGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGt                                          <  2:81468/238‑1 (MQ=255)
                                                                                                                                                                                           aaaTTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCTCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCTGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGt                                          <  2:100025/238‑1 (MQ=255)
                                                                                                                                                                                           aaaTTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGt                                          <  1:767429/238‑1 (MQ=255)
                                                                                                                                                                                           aaaTTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGt                                          <  1:747258/238‑1 (MQ=255)
                                                                                                                                                                                           aaaTTAACCATCAGGCAAAAACGGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGt                                          <  2:396679/238‑1 (MQ=255)
                                                                                                                                                                                           aaaTGAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGt                                          <  2:850897/238‑1 (MQ=255)
                                                                                                                                                                                              ttAACCATCATGCGAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGt                                      <  2:851714/239‑1 (MQ=255)
                                                                                                                                                                                                                   atAAAGAGGCGGGGGCGGCGGGCACGGTCGAGGGCGAGCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCACCAGTACGGTATTGCCGCCGTGGCGCTACATAACCCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGGGGCGGGGTTTGTCTCTATCCACTTTTTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTc                 <  2:545498/238‑1 (MQ=255)
                                                                                                                                                                                                                     aaaGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCa               <  1:662848/239‑1 (MQ=255)
                                                                                                                                                                                                                                   ggcggTCCCGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCAGCAGCCCGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGc  <  2:699259/238‑1 (MQ=255)
                                                                                                                                                                                                                                   ggcggTCACTCTGGAGGGCGATCGCGCATTTGGGCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGc  <  2:466746/238‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           
GGTCATCGCTTTGATGCTCAGACGCTGCACAGTTTTATTCAGGCTGTATTTCGTCAGATGGGTAGCGAGGAACAAGAAGCGAAATTAGTTGCCGATCATTTAATCGCGGCAAACCTGGCAGGGCATGATTCACATGGTATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGC  >  NC_000913/836360‑836824

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