Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,017,882 A→T 5.3% intergenic (‑175/+59) PP_5679 ← / ← PP_4421 hypothetical protein/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,017,8820AT 5.3% 441.6 / 9.3 131intergenic (‑175/+59)PP_5679/PP_4421hypothetical protein/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (58/66);  new base T (3/4);  total (61/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 135 total aligned reads displayed.

GAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAGAGGA  >  NC_002947/5017735‑5018031
                                                                                                                                                    |                                                                                                                                                     
gagatgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTaaa                                                                                                                                                     >  1:3967729/1‑150 (MQ=255)
    tgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTTGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCgg                                                                                                                                                <  2:524530/150‑1 (MQ=255)
     gaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCgg                                                                                                                                                >  1:5034480/1‑150 (MQ=255)
     gaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCgg                                                                                                                                                <  1:1635118/150‑1 (MQ=255)
     gaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCgg                                                                                                                                                >  1:3219462/1‑150 (MQ=255)
      aaTACGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTTATGACAGCGTCTTGGCTTTGATTTTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCggg                                                                                                                                               <  2:3048379/150‑1 (MQ=255)
      aaTAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCggg                                                                                                                                               >  2:5348104/1‑150 (MQ=255)
       aaaaGCCCGCACCGGTATCCGGCTCGCGATATTACTTTCCGCCCATTCGACATTATCTCACCCGAAACAACTACCTTGTGAACTCATCACACCTACTTGCAATTGATTGTGAAACAACAACAAAAAAAAAAAGCCCGCCTAAAGGCGGGc                                                                                                                                              <  2:2863425/148‑1 (MQ=255)
        tacgcccgcccCGGTATCCGGCGGGGGCTTTTACTTTCCGGGCCTTGTACTGTATCTTACCCGAAAGTAGGTCCTTGTGAACTGTTTACAGCTTCTGGGCTTTTATTGTTTTAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCt                                                                                                                                             <  1:1351790/147‑1 (MQ=255)
         aaGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCtt                                                                                                                                            <  1:436866/150‑1 (MQ=255)
           gCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTa                                                                                                                                           >  2:4310412/1‑150 (MQ=255)
           gCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTTTAGGCGGGCtttt                                                                                                                                          <  1:2051957/150‑3 (MQ=255)
             ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCATTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTaaaa                                                                                                                                        <  2:1987696/150‑1 (MQ=255)
                 aCCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAAAGaa                                                                                                                                    >  2:10780/1‑150 (MQ=255)
                 aCCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGAGGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGaa                                                                                                                                    <  2:1816858/150‑1 (MQ=255)
                   cGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAaga                                                                                                                                  <  1:1019057/150‑1 (MQ=255)
                     gTCTCCGGGGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAATTTACGTCCTTGTTAACTTATGACAGCGTCTTGGCTTTTATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAgg                                                                                                                                <  2:3869948/150‑1 (MQ=255)
                     gTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAgg                                                                                                                                >  2:4551177/1‑150 (MQ=255)
                     gTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAgg                                                                                                                                >  1:4630478/1‑150 (MQ=255)
                     gTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAgg                                                                                                                                <  1:4447117/150‑1 (MQ=255)
                     gTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAgg                                                                                                                                >  2:4551200/1‑150 (MQ=255)
                     gTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAgg                                                                                                                                <  2:3869958/150‑1 (MQ=255)
                           ggCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAAtt                                                                                                                          <  2:234634/150‑1 (MQ=255)
                           ggCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAAtt                                                                                                                          <  2:3019688/150‑1 (MQ=255)
                            gcgcGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGACGAGGGTAATTc                                                                                                                         >  2:4844751/1‑150 (MQ=255)
                             cgcgGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCa                                                                                                                        <  1:193708/150‑1 (MQ=255)
                                    ttttACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGca                                                                                                                 >  1:2927866/1‑150 (MQ=255)
                                      ttACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGcatc                                                                                                               >  2:5124224/1‑150 (MQ=255)
                                        aCTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGcatcat                                                                                                             >  2:3660244/1‑150 (MQ=255)
                                        aCTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGcatcat                                                                                                             >  1:1885563/1‑150 (MQ=255)
                                         cTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGTCAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCAtt                                                                                                            <  1:2717361/150‑1 (MQ=255)
                                          tttCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAATCCCGCCAAAAGGCGGGCTTAAAAAGAAGAGGGTAATTCAGCTGGCATCATTg                                                                                                           >  1:4123181/1‑150 (MQ=255)
                                          tttCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTg                                                                                                           <  2:1780330/150‑1 (MQ=255)
                                          tttCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTg                                                                                                           <  1:698049/150‑1 (MQ=255)
                                          tttCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTCGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGc                                                                                                          >  2:2165713/1‑150 (MQ=255)
                                           ttCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGc                                                                                                          <  2:2548168/150‑1 (MQ=255)
                                                  ccTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTGAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGtg                                                                                                   >  2:4101844/1‑149 (MQ=255)
                                                    ttGGACAGTATCTGACCCGAAAGTACTTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGa                                                                                                 <  1:3349242/150‑1 (MQ=255)
                                                    ttGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTTTAGGCGGGCtttttt                                                                                                                                        >  1:3781485/1‑107 (MQ=255)
                                                    ttGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTTTAGGCGGGCtttttt                                                                                                                                        <  2:3781485/111‑5 (MQ=255)
                                                    ttGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt                                                                                                <  2:2203663/150‑1 (MQ=255)
                                                     tGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt                                                                                                >  1:3066383/1‑150 (MQ=255)
                                                      tgACAGTATCTTACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt                                                                                                <  2:1521926/149‑1 (MQ=255)
                                                         cAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCt                                                                                            <  1:1062364/150‑1 (MQ=255)
                                                           gTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCa                                                                                          <  1:2633740/150‑1 (MQ=255)
                                                              tCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTTATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGAc                                                                                       <  1:2403800/150‑1 (MQ=255)
                                                               cTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGAcc                                                                                      >  1:1197752/1‑150 (MQ=255)
                                                               cTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGAcc                                                                                      <  1:1175560/150‑1 (MQ=255)
                                                                tGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACTAAGAGGGTAATACAGCTGGCATCATTGCAAAATTTGATGGCTCAGAcat                                                                                     >  2:2732502/1‑148 (MQ=255)
                                                                 gACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTTTAGGCGGGCtttt                                                                                                                                          >  2:2305328/1‑95 (MQ=255)
                                                                 gACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTTTAGGCGGGCtttt                                                                                                                                          <  1:2305328/97‑3 (MQ=255)
                                                                  aCCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGTTTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCt                                                                                                                     <  1:3027215/116‑1 (MQ=255)
                                                                  aCCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCt                                                                                                                     >  2:3027215/1‑116 (MQ=255)
                                                                  aCCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTc                                                                                   <  2:1111261/150‑1 (MQ=255)
                                                                   cccGAAAGTACGTCCTTGTGAACTGATGACAGCTTATTGGCTTTTATTGTTATAGAACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGTCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCa                                                                                  <  2:318573/150‑1 (MQ=255)
                                                                    ccGAAAGTACGTCCTTGTGAACTGATGTCAGCGTCTTGGCTTTGATTGTGCTAGCACTACAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCa                                                                                  <  2:1334493/150‑1 (MQ=255)
                                                                    ccGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAg                                                                                 <  1:848957/150‑1 (MQ=255)
                                                                    ccGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAg                                                                                 >  2:306889/1‑150 (MQ=255)
                                                                     cGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGc                                                                                <  1:2947706/150‑1 (MQ=255)
                                                                      gAAATTACGTCCTTTTGACCTGATGACAGCGTCTTGGCTTTGATTGTGATAGAACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCt                                                                               <  1:352952/150‑1 (MQ=255)
                                                                      gAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtc                                                                              >  1:4075795/1‑150 (MQ=255)
                                                                         aGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtcat                                                                            <  2:730633/150‑1 (MQ=255)
                                                                            aCGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAATACCCCCAAAAGGCGGGCTGAAAAAGAAGAGGGGGATTAAGCTGGCAGCATATACAACGATGATGTCTCAGACCTACTGCTCATCaa                                                                         >  1:644817/1‑150 (MQ=255)
                                                                              gTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAgg                                                                       >  1:4892187/1‑150 (MQ=255)
                                                                                ccTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGAc                                                                     <  2:2999108/150‑1 (MQ=255)
                                                                                 cTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCt                                                                  <  2:3647207/150‑1 (MQ=255)
                                                                                    gtgAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTc                                                                 >  1:1370428/1‑150 (MQ=255)
                                                                                     tgAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCa                                                                <  2:3967729/150‑1 (MQ=255)
                                                                                     tgAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAg                                                               >  1:3927067/1‑150 (MQ=255)
                                                                                          tgatgcCAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGcaca                                                           <  1:5394291/150‑1 (MQ=255)
                                                                                            aggaCAGCGTCTTGGCTTTGATTGTGATAGCACTCCAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGc                                                         <  1:5380017/148‑1 (MQ=255)
                                                                                                cAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt                                                                                                >  2:4958602/1‑107 (MQ=255)
                                                                                                cAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt                                                                                                <  1:4958602/107‑1 (MQ=255)
                                                                                                 aGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgt                                                    >  1:201048/1‑150 (MQ=255)
                                                                                                 aGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgt                                                    >  1:1416812/1‑150 (MQ=255)
                                                                                                 aGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgt                                                    >  2:3387748/1‑150 (MQ=255)
                                                                                                  gCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgtgt                                                  <  1:1178964/150‑1 (MQ=255)
                                                                                                    gTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGtgtg                                                 <  1:4149474/150‑1 (MQ=255)
                                                                                                     tCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGc                                                >  2:719575/1‑150 (MQ=255)
                                                                                                           cTTGGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCt                                          <  2:4648717/150‑1 (MQ=255)
                                                                                                            tttGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtg                                         >  2:4379759/1‑150 (MQ=255)
                                                                                                             ttGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtga                                        <  1:4551200/150‑1 (MQ=255)
                                                                                                              tGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtgat                                       <  1:4844751/150‑1 (MQ=255)
                                                                                                               gATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtcatc                                                                           >  2:1210085/1‑113 (MQ=255)
                                                                                                                              aCAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTcc                       >  1:4706015/1‑150 (MQ=255)
                                                                                                                                 aaaaaaaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGGAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGAACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTg                     >  1:3336836/1‑150 (MQ=255)
                                                                                                                                   aaaaaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGtt                   <  1:1033262/150‑1 (MQ=255)
                                                                                                                                    aaaaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTc                  >  1:427724/1‑150 (MQ=255)
                                                                                                                                     aaaaaaaGCCCGCCTAAAGGCGGGCATAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCg                 <  1:5348104/150‑1 (MQ=255)
                                                                                                                                      aaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCgg                >  2:3681151/1‑150 (MQ=255)
                                                                                                                                       aaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGt               >  2:2142638/1‑150 (MQ=255)
                                                                                                                                        aaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTc              >  2:1516745/1‑150 (MQ=255)
                                                                                                                                            gCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACGTCCTGTTCGGTCGTGa          <  2:1671212/150‑1 (MQ=255)
                                                                                                                                             cccGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGAc         <  1:3981627/150‑1 (MQ=255)
                                                                                                                                              ccGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTgatgat                                    <  2:1284176/122‑1 (MQ=255)
                                                                                                                                              ccGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTgatgat                                    >  1:1284176/1‑122 (MQ=255)
                                                                                                                                              ccGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACa        >  1:3207898/1‑150 (MQ=255)
                                                                                                                                                gCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAag      >  1:4100548/1‑150 (MQ=255)
                                                                                                                                                   tAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAgagg   >  1:1459346/1‑150 (MQ=255)
                                                                                                                                                    aaaGGCGGGCATAACACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAgagga  <  2:3422290/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCTGTTCGGTCGTGACAAGAGGA  >  NC_002947/5017735‑5018031

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