Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,013,335 G→C 11.1% L36L (CTG→CTC PP_5678 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,013,3350GC11.1% 261.7 / 21.0 91L36L (CTG→CTCPP_5678hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/37);  new base C (6/4);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

ACCCGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTCT  >  NC_002947/5013187‑5013483
                                                                                                                                                    |                                                                                                                                                    
aCCCGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCTCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGa                                                                                                                                                     >  1:1586750/1‑150 (MQ=255)
aCCCGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCTCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGa                                                                                                                                                     >  1:1586721/1‑150 (MQ=255)
     aGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGcgcg                                                                                                                                                >  1:860248/1‑150 (MQ=255)
       cAGTGCAAAAAGCCTCAGACTTAGGGGAGAGACAAGGATGATTAAAGATTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCacggc                                                                                                                                              <  2:327588/150‑1 (MQ=255)
         gTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGTCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGt                                                                                                                                            <  2:1286866/150‑1 (MQ=255)
         gTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCACGGCGt                                                                                                                                            <  2:1342510/150‑1 (MQ=255)
          tGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCTCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTg                                                                                                                                           <  2:1586721/150‑1 (MQ=255)
          tGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCTCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTg                                                                                                                                           <  2:1586750/150‑1 (MQ=255)
          tGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTg                                                                                                                                           >  1:252350/1‑150 (MQ=255)
            caacaaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGa                                                                                                                                         <  1:580268/150‑1 (MQ=255)
             aacaaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGaa                                                                                                                                        >  1:1561725/1‑150 (MQ=255)
              acaaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGaaa                                                                                                                                       >  2:1392900/1‑150 (MQ=255)
                aaGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATc                                                                                                                                     >  2:2221523/1‑150 (MQ=255)
                  gCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTCTAATTCAg                                                                                                                                   <  2:692460/150‑1 (MQ=255)
                                 gagaGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTATTGTCTGGTATGACGGCCATCCACAGCACGCCGTGAAAGCCGCGGCGTGCAAATCAGACCCGCAGCAAAAcc                                                                                                                    >  1:414186/1‑150 (MQ=255)
                                   gagCCTATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGa                                                                                                                                                     <  2:1482258/115‑1 (MQ=255)
                                   gagCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGa                                                                                                                                                     >  1:1482258/1‑115 (MQ=255)
                                       cAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAg                                                                                                              >  1:1254424/1‑150 (MQ=255)
                                           gatgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCTGGCTCGCAGCAAAACCCCGAAGCCTa                                                                                                          >  1:2341859/1‑150 (MQ=255)
                                           gatgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACATGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTa                                                                                                          >  2:748843/1‑150 (MQ=255)
                                            atgatCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTAt                                                                                                         <  2:1198672/150‑1 (MQ=255)
                                                tCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCt                                                                                                           <  2:2422745/144‑1 (MQ=255)
                                                tCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCt                                                                                                           >  1:2422745/1‑144 (MQ=255)
                                                tCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAAc                                                                                                     >  1:1326926/1‑150 (MQ=255)
                                                   aaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTcc                                                                                                  >  1:1939802/1‑150 (MQ=255)
                                                   aaGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTcc                                                                                                  >  2:1191437/1‑150 (MQ=255)
                                                           aCCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGc                                                                                          <  1:1392900/150‑1 (MQ=255)
                                                             ccTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCga                                                                                        <  2:252350/150‑1 (MQ=255)
                                                               tCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCgaga                                                                                      <  2:672327/150‑1 (MQ=255)
                                                                  cAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAAcc                                                                                                                    >  2:1313206/1‑117 (MQ=255)
                                                                  cAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAAcc                                                                                                                    <  1:1313206/117‑1 (MQ=255)
                                                                       gTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCacggcg                                                                                                                                             <  1:1365992/87‑1 (MQ=255)
                                                                       gTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGCGCTCAGCacggcg                                                                                                                                             >  2:1365992/1‑87 (MQ=255)
                                                                         tGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGcttctt                                                                            <  2:860248/150‑1 (MQ=255)
                                                                          gCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCttt                                                                           <  2:890692/150‑1 (MQ=255)
                                                                          gCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCttt                                                                           <  1:130216/150‑1 (MQ=255)
                                                                            ggTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGcc                                                                        >  1:2085063/1‑150 (MQ=255)
                                                                             gTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGcc                                                                        <  2:2390597/150‑1 (MQ=255)
                                                                              tCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCa                                                                       <  1:443691/150‑1 (MQ=255)
                                                                                  cGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGTGAGAGCTGCTTCTTTGCCACGCt                                                                   >  1:1534932/1‑150 (MQ=255)
                                                                                   gTAGTTTGCGACCGTTGCGGCCAGACTTGTACTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAgg                                                                                                                                  >  2:329299/1‑86 (MQ=255)
                                                                                   gTAGTTTGCGACCGTTGCGGCCAGACTTGTACTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAgg                                                                                                                                  <  1:329299/86‑1 (MQ=255)
                                                                                     aGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTc                                                                >  1:2394535/1‑150 (MQ=255)
                                                                                       tttGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCaa                                                              >  1:2321561/1‑150 (MQ=255)
                                                                                        ttGCGACCGTTGCGGCCAGACTTGTCCTCCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAAc                                                             <  1:1407078/150‑1 (MQ=255)
                                                                                            gACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCGCGCCGTGCTGAGCGCGGCTTGGGAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGGACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATa                                                         >  2:672479/1‑150 (MQ=255)
                                                                                                gTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGc                                                     >  2:450057/1‑150 (MQ=255)
                                                                                                  tGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGt                                                   <  2:2321561/150‑1 (MQ=255)
                                                                                                  tGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGt                                                   <  2:1326926/150‑1 (MQ=255)
                                                                                                      gCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAAc                                               >  2:141078/1‑150 (MQ=255)
                                                                                                        cAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGc                                             <  1:334138/150‑1 (MQ=255)
                                                                                                          gACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTg                                           <  2:1561725/150‑1 (MQ=255)
                                                                                                           aCTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTgg                                          >  2:1992916/1‑150 (MQ=255)
                                                                                                            cTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTggg                                         >  2:1476115/1‑150 (MQ=255)
                                                                                                            cTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTggg                                         <  1:450057/150‑1 (MQ=255)
                                                                                                             ttGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGc                                        >  1:2442777/1‑150 (MQ=255)
                                                                                                             ttGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTACATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGc                                        >  2:1225894/1‑150 (MQ=255)
                                                                                                                  cTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCt                                                                                              >  2:2188228/1‑91 (MQ=255)
                                                                                                                  cTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCt                                                                                              <  1:2188228/91‑1 (MQ=255)
                                                                                                                  cTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAAt                                   <  1:1191437/150‑1 (MQ=255)
                                                                                                                   ttCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAAtt                                  >  2:28462/1‑150 (MQ=255)
                                                                                                                      tgatgaCGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGAt                               <  1:2221523/150‑1 (MQ=255)
                                                                                                                        atgaCGGCCATCCACAGCACGCCCCATTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGt                             >  1:1519986/1‑150 (MQ=255)
                                                                                                                          gaCGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAACCGAACCCCGAAGACTATAAACCCTGTCTTTTCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTtc                           >  1:1485360/1‑150 (MQ=255)
                                                                                                                           aCGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTAGTAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTtct                          <  1:748843/150‑1 (MQ=255)
                                                                                                                             ggCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTctct                        >  1:782672/1‑150 (MQ=255)
                                                                                                                               ccATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGa                      <  2:1896082/150‑1 (MQ=255)
                                                                                                                               ccATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGa                      >  1:2359379/1‑150 (MQ=255)
                                                                                                                                cATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCaa                                    <  2:473233/135‑1 (MQ=255)
                                                                                                                                cATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCaa                                    >  1:473233/1‑135 (MQ=255)
                                                                                                                                  tCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCttt                                                                                            >  1:635530/1‑77 (MQ=255)
                                                                                                                                  tCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCttt                                                                                            <  2:635530/77‑1 (MQ=255)
                                                                                                                                  tCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGa                   >  2:2383921/1‑150 (MQ=255)
                                                                                                                                  tCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGa                   >  1:734664/1‑150 (MQ=255)
                                                                                                                                  cccACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAgctgc                                                                                 <  1:489003/87‑1 (MQ=255)
                                                                                                                                  cccACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAgctgc                                                                                 >  2:489003/2‑88 (MQ=255)
                                                                                                                                   ccACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGaa                  <  1:2214292/150‑1 (MQ=255)
                                                                                                                                     acaGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGc                                                                                          <  1:2954/76‑1 (MQ=255)
                                                                                                                                     acaGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGc                                                                                          >  2:2954/1‑76 (MQ=255)
                                                                                                                                     acaGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGc                <  1:469020/150‑1 (MQ=255)
                                                                                                                                     acaGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGc                >  2:1780753/1‑150 (MQ=255)
                                                                                                                                     acaGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTATTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGc                <  2:1939802/150‑1 (MQ=255)
                                                                                                                                       tccACGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTg              >  1:1563399/3‑150 (MQ=255)
                                                                                                                                         cACGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCgaga                                                                                      <  2:2115086/76‑1 (MQ=38)
                                                                                                                                         cACGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCgaga                                                                                      >  1:2115086/1‑76 (MQ=38)
                                                                                                                                         cACGCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGAc            >  1:55404/1‑150 (MQ=255)
                                                                                                                                            gCCGCGCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTcc         >  1:226027/1‑150 (MQ=255)
                                                                                                                                               gTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCAGCTTTGCGAGAGCTGCTTCTAAGCGACGCTTTCAACGATAAAGCGTGAACGCTGGGCACAATTGATGTTCTGTGATGAAGCTTACTCCTTa      >  1:748394/1‑150 (MQ=255)
                                                                                                                                               gTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTACGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCACCGATAAAGCGTGAACGCTGGGCCCAATCGATGTTCTCTGAGGAAGCTGACTCCTTa      <  2:1254424/150‑1 (MQ=255)
                                                                                                                                                 gCTCAGCACGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGagag                                                                                     >  1:1452457/1‑69 (MQ=255)
                                                                                                                                                   tGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCNGAAGCCTATNAACTCCATCTTTGCAAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTActct  >  1:53967/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCCGAGCAGTGCAACAAGCCTCAGACTTAGGGGAGAGCCAATGATGATCAAAGAGTCGACCCTCGCAGCCGTTGCGGTCAGCGTAGTTTGCGACCGTTGCGGCCAGACTTGTCCTTCTGATGACGGCCATCCACAGCACGCCGTGCTGAGCGCGGCGTGGAAATCAGGCTCGCAGCAAAACCCCGAAGCCTATGAACTCCATCTTTGCGAGAGCTGCTTCTTTGCCACGCTTTCAACGATAAAGCGTGAACGCTGGGCCCAATTGATGTTCTCTGAGGAAGCTGACTCCTTACTCT  >  NC_002947/5013187‑5013483

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