Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 264,065 G→C 10.4% K272N (AAG→AAC) ‡ gabD‑I → succinate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947264,0650GC10.4% 636.0 / 51.0 203K272N (AAG→AAC) ‡gabD‑Isuccinate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (99/82);  new base C (8/13);  total (108/95)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 204 total aligned reads displayed.

TCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCT  >  NC_002947/263919‑264214
                                                                                                                                                  |                                                                                                                                                     
tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                                                                                                                                    >  1:1090264/1‑150 (MQ=255)
tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                                                                                                                                    >  1:943371/1‑150 (MQ=255)
tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                                                                                                                                    <  1:2850392/150‑1 (MQ=255)
tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGAGGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                                                                                                                                    >  2:2834786/1‑150 (MQ=255)
   cTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTc                                                                                                                                                 >  2:38943/1‑150 (MQ=255)
     ggTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGa                                                                                                                                               <  2:4990115/150‑1 (MQ=255)
       tACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAgg                                                                                                                                             <  1:5343723/150‑1 (MQ=255)
          gCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGATGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGgcg                                                                                                                                          <  2:3176330/150‑1 (MQ=255)
          gCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGgcg                                                                                                                                          >  2:832913/1‑150 (MQ=255)
          gCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGgcg                                                                                                                                          >  2:4742049/1‑150 (MQ=255)
             aGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGa                                                                                                                                       <  1:2118462/150‑1 (MQ=255)
              gCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGat                                                                                                                                      >  1:2261179/1‑150 (MQ=255)
               cTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGatc                                                                                                                                     <  1:4607466/150‑1 (MQ=255)
                ggTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGttca                                                                                                                                    >  1:398575/2‑150 (MQ=255)
                  tCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGatcatc                                                                                                                                  >  1:3074263/1‑150 (MQ=255)
                       ctccGGCTCGACCGAAATCGGTCGCCAGCTGATGGATGTATGCGCCAAGGACTTCAATCAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATGGTGTTCTTCGACGCCGACCTGGACAAGGCGGTAGAGGGCGCGATCTTCTCCaa                                                                                                                             <  1:4904290/148‑1 (MQ=255)
                       cACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCaa                                                                                                                             <  2:3074263/150‑1 (MQ=255)
                       cACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCaa                                                                                                                             <  2:199834/150‑1 (MQ=255)
                          cGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTa                                                                                                                          <  2:684735/150‑1 (MQ=255)
                          cGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTa                                                                                                                          >  1:2914906/1‑150 (MQ=255)
                                gACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaa                                                                                                                    >  1:3921996/1‑150 (MQ=255)
                                gACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaa                                                                                                                    >  2:2265247/1‑150 (MQ=255)
                                 aCCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaac                                                                                                                   <  1:45108/150‑1 (MQ=255)
                                 aCCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaac                                                                                                                   >  1:1978527/1‑150 (MQ=255)
                                   cGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGGCGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaacaa                                                                                                                 >  1:1826817/1‑150 (MQ=255)
                                   cGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaacaa                                                                                                                 >  2:3581730/1‑150 (MQ=255)
                                    gAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaacaac                                                                                                                <  1:1706209/150‑1 (MQ=255)
                                      aaTCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACgg                                                                                                              >  2:3201307/1‑150 (MQ=255)
                                         cGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCa                                                                                                           <  2:1102981/150‑1 (MQ=255)
                                           gTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGa                                                                                                         <  1:3675272/150‑1 (MQ=255)
                                           gTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGa                                                                                                         >  2:5397655/1‑150 (MQ=255)
                                           gTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGa                                                                                                         >  2:5052712/1‑150 (MQ=255)
                                            tCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACACCgg                                                                                                              >  2:4419445/1‑144 (MQ=255)
                                            tCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACgg                                                                                                              <  1:4419445/144‑1 (MQ=255)
                                              gCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCt                                                                                                      <  1:4048896/150‑1 (MQ=255)
                                              gCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCt                                                                                                      <  2:4637124/150‑1 (MQ=255)
                                              gCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCACCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCt                                                                                                      <  1:3401791/150‑1 (MQ=255)
                                                cAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGc                                                                                                    <  1:4122776/150‑1 (MQ=255)
                                                cAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGc                                                                                                    <  2:678214/150‑1 (MQ=255)
                                                cAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGc                                                                                                    <  1:3201307/150‑1 (MQ=255)
                                                cAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGc                                                                                                    <  1:5181560/150‑1 (MQ=255)
                                                cAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGc                                                                                                    >  2:1339455/1‑150 (MQ=255)
                                                   cTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTc                                                                                                 <  2:4478204/150‑1 (MQ=255)
                                                    tGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCt                                                                                                <  2:2330708/150‑1 (MQ=255)
                                                        ggAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTgcgc                                                                                            <  2:5250390/150‑1 (MQ=255)
                                                        ggAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTgcgc                                                                                            >  1:2365814/1‑150 (MQ=255)
                                                        ggAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTgcgc                                                                                            >  2:120222/1‑150 (MQ=255)
                                                        ggAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTgcgc                                                                                            >  1:810993/1‑150 (MQ=255)
                                                            gaaTGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAAc                                                                                        <  2:134920/150‑1 (MQ=255)
                                                                 cgcCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTAt                                                                                   <  1:5061710/150‑1 (MQ=255)
                                                                   cgccgaccATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTctct                                                                                 <  1:2777033/143‑1 (MQ=255)
                                                                   cgccgaccATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTctct                                                                                 <  1:2777218/143‑1 (MQ=255)
                                                                    cAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTa                                                                                >  1:1296172/1‑150 (MQ=255)
                                                                     aaGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTAc                                                                               >  1:2052466/1‑150 (MQ=255)
                                                                        gACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTg                                                                            <  1:432161/150‑1 (MQ=255)
                                                                        accATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGta                                                                            <  2:4836722/148‑2 (MQ=255)
                                                                          cATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACa                                                                          >  1:1110172/1‑150 (MQ=255)
                                                                          cATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACAGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACa                                                                          <  2:713263/150‑1 (MQ=255)
                                                                           aTCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGc                                                                                                             >  1:1538259/1‑114 (MQ=255)
                                                                           aTCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAg                                                                         <  2:3523272/150‑1 (MQ=255)
                                                                              aagaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGAc                                                                      >  2:1429239/1‑150 (MQ=255)
                                                                              aagaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGAc                                                                      >  1:835213/1‑150 (MQ=255)
                                                                                gaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACgg                                                                    >  2:2024866/1‑150 (MQ=255)
                                                                                         ccTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTacga                                                           <  1:4495837/150‑1 (MQ=255)
                                                                                          cTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTacgac                                                          >  1:4218684/1‑150 (MQ=255)
                                                                                          cTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTacgac                                                          <  1:3794638/150‑1 (MQ=255)
                                                                                               gCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGtt                                                     >  1:755668/1‑150 (MQ=255)
                                                                                                cttggcggCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCttt                                                    <  1:5409097/144‑2 (MQ=255)
                                                                                                   ggTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTc                                                    <  1:1256096/147‑1 (MQ=255)
                                                                                                    gTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCg                                                <  2:1534841/150‑1 (MQ=255)
                                                                                                     tGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCga                                               <  1:2666721/150‑1 (MQ=255)
                                                                                                       gCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGTGTTAGCCgaga                                             >  1:3481564/1‑150 (MQ=255)
                                                                                                       gCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCgaga                                             >  1:2286302/1‑150 (MQ=255)
                                                                                                        cAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGaa                                            >  1:1152410/1‑150 (MQ=255)
                                                                                                        cAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGaa                                            <  2:2914906/150‑1 (MQ=255)
                                                                                                            gcgcctttCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCtc                                        <  1:5152214/144‑2 (MQ=255)
                                                                                                            gCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTg                                        <  2:4590625/150‑1 (MQ=255)
                                                                                                            gCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTg                                        <  2:4913864/150‑1 (MQ=255)
                                                                                                              cccGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTa                                                                                    >  1:1068363/1‑104 (MQ=255)
                                                                                                                    cATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGc                                >  2:185900/1‑150 (MQ=255)
                                                                                                                     aTCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACg                                                                     <  2:3557660/112‑1 (MQ=255)
                                                                                                                     aTCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACg                                                                     >  1:3557660/1‑112 (MQ=255)
                                                                                                                     aTCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCa                               <  2:5159840/150‑1 (MQ=255)
                                                                                                                          gtTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGc                          >  2:3427562/1‑150 (MQ=255)
                                                                                                                          gtTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGc                          <  1:5398953/150‑1 (MQ=255)
                                                                                                                           ttCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGc                                          >  2:4216746/1‑133 (MQ=255)
                                                                                                                           ttCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGcc                         >  1:2389793/1‑150 (MQ=255)
                                                                                                                             cgacgaCGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCaa                       <  2:2261091/150‑1 (MQ=255)
                                                                                                                              gacgacGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCt                                                                                                      >  1:1284288/1‑70 (MQ=255)
                                                                                                                              gacgacGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCt                                                                                                      <  2:1284288/70‑1 (MQ=255)
                                                                                                                               acgacgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGc                     >  1:1731880/1‑150 (MQ=255)
                                                                                                                               acgacgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGc                     >  1:4475376/1‑150 (MQ=255)
                                                                                                                                 gacgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTg                   >  2:1024/1‑150 (MQ=255)
                                                                                                                                   cgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGaa                 >  2:618293/1‑150 (MQ=255)
                                                                                                                                    gCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAg                >  2:2351320/1‑150 (MQ=255)
                                                                                                                                     ccGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGa               >  2:4714243/1‑150 (MQ=255)
                                                                                                                                      cGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGAt              <  1:4871576/150‑1 (MQ=255)
                                                                                                                                       gACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATc             <  1:4742049/150‑1 (MQ=255)
                                                                                                                                                aaGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGc    >  2:3700627/1‑150 (MQ=255)
                                                                                                                                                aaGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGcc   >  2:1638463/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCT  >  NC_002947/263919‑264214

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