Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,763,687 A→G 5.9% E357G (GAG→GGG) ‡ PP_3325 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,763,6870AG 5.9% 638.5 / 18.0 201E357G (GAG→GGG) ‡PP_3325ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (95/94);  new base G (6/6);  total (101/100)
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.70e-01

Only 100 of 208 total aligned reads displayed.

CCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACC  >  NC_002947/3763538‑3763835
                                                                                                                                                     |                                                                                                                                                    
ccTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGa                                                                                                                                                      <  1:1400518/150‑1 (MQ=255)
ccTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCGTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGa                                                                                                                                                      <  2:2294900/150‑1 (MQ=255)
 cTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAg                                                                                                                                                     >  1:4800413/1‑150 (MQ=255)
 cTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAg                                                                                                                                                     >  1:3311587/1‑150 (MQ=255)
 cTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAg                                                                                                                                                     >  1:2026566/1‑150 (MQ=255)
     ccAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTgg                                                                                                                                                 <  1:2691892/150‑1 (MQ=255)
     ccAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTgg                                                                                                                                                 <  2:4201060/150‑1 (MQ=255)
     ccAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTgg                                                                                                                                                 <  1:493072/150‑1 (MQ=255)
      cAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGCTCGAACTGCCCGTGCAGCTCGGCCTggg                                                                                                                                                <  1:1227905/150‑1 (MQ=38)
       aaCCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGTATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGt                                                                                                                                               <  1:2605319/150‑1 (MQ=255)
       aaCCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGt                                                                                                                                               >  1:2984033/1‑150 (MQ=255)
        cccACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCCGGCCGAGCTGCACGGGCAGTTCGAGCTGGGtg                                                                                                                                              <  1:2633499/149‑1 (MQ=255)
        aCCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACACCACCCAGCTCGAACTGCCCGTGCAGCTCGGCCTGGGtg                                                                                                                                              <  2:351892/150‑1 (MQ=38)
          cACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGtggc                                                                                                                                            >  1:5345136/1‑150 (MQ=255)
            cTACACCCAGGGCACGCCCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGtggctg                                                                                                                                          >  2:2002741/1‑150 (MQ=255)
              acacCCAGGGCCCGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGtggctggc                                                                                                                                        >  2:2096842/1‑150 (MQ=255)
                acCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGgcag                                                                                                                                      >  2:1061688/1‑150 (MQ=255)
                acCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGgcag                                                                                                                                      >  1:3967790/1‑150 (MQ=255)
                  ccAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGa                                                                                                                                                      >  1:4235808/1‑132 (MQ=255)
                  ccAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGa                                                                                                                                                      <  2:4235808/132‑1 (MQ=255)
                  ccAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGgcagca                                                                                                                                    >  2:712086/1‑150 (MQ=255)
                  ccAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGgcagca                                                                                                                                    >  1:1385615/1‑150 (MQ=255)
                   cAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCAt                                                                                                                                   <  1:1225258/150‑1 (MQ=255)
                     gggCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCa                                                                                                                                 >  1:3665364/1‑150 (MQ=255)
                              cAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCtggt                                                                                                                        <  2:624065/150‑1 (MQ=255)
                               aaGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCtggtg                                                                                                                       >  1:506195/1‑150 (MQ=255)
                                aGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCtggtgg                                                                                                                      <  2:3332056/150‑1 (MQ=255)
                                aGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCtggtgg                                                                                                                      <  1:3613295/150‑1 (MQ=255)
                                 gggCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTggg                                                                                                                     >  1:4826863/1‑150 (MQ=255)
                                  ggCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGAAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTgggg                                                                                                                    <  2:4791564/150‑1 (MQ=255)
                                   gCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGc                                                                                                                   <  2:2564278/150‑1 (MQ=255)
                                   gCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGc                                                                                                                   <  2:77184/150‑1 (MQ=255)
                                   gCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGc                                                                                                                   <  2:160179/150‑1 (MQ=255)
                                     gaCAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTg                                                                                                                 >  1:2790649/2‑150 (MQ=255)
                                     aaCAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTg                                                                                                                 <  1:3761545/150‑1 (MQ=255)
                                     aaCAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTg                                                                                                                 <  2:3696296/150‑1 (MQ=255)
                                     aaCAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTg                                                                                                                 <  1:1331100/150‑1 (MQ=255)
                                       cAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGa                                                                                                               <  2:2126066/150‑1 (MQ=255)
                                        aGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGa                                                                                                                               <  2:2478340/133‑1 (MQ=255)
                                        aGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGa                                                                                                                               >  1:2478340/1‑133 (MQ=255)
                                        aGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAg                                                                                                              >  2:539904/1‑150 (MQ=255)
                                         gCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGt                                                                                                             <  2:3781476/150‑1 (MQ=255)
                                          cTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACAGCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTa                                                                                                            <  2:3579158/150‑1 (MQ=255)
                                           tCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTAc                                                                                                           >  2:28182/1‑150 (MQ=255)
                                           tCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTAc                                                                                                           >  2:2959208/1‑150 (MQ=255)
                                             cGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGTACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGa                                                                                                         >  2:5294562/1‑150 (MQ=255)
                                              gAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGaa                                                                                                        <  1:1449402/150‑1 (MQ=255)
                                                 ccccAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAAc                                                                                                     >  2:2705903/1‑150 (MQ=255)
                                                     aGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTAcc                                                                                                 <  2:2872084/150‑1 (MQ=255)
                                                                aGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAg                                                                                      <  1:1360196/150‑1 (MQ=255)
                                                                 gCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGa                                                                                     >  1:5279970/1‑150 (MQ=255)
                                                                   accaccATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAg                                                                                   >  2:2739324/1‑150 (MQ=255)
                                                                     caccaTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCg                                                                                                >  2:5308272/1‑135 (MQ=255)
                                                                     caccaTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCg                                                                                                <  1:5308272/135‑1 (MQ=255)
                                                                     caccaTCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTa                                                                                 >  2:3463835/1‑150 (MQ=255)
                                                                                  ttCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGcc                                                                    >  2:116463/1‑150 (MQ=255)
                                                                                     tATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTgggg                                                                                                                    <  1:2184209/99‑1 (MQ=255)
                                                                                     tATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTgggg                                                                                                                    >  2:2184209/1‑99 (MQ=255)
                                                                                       tCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCACAGCGCCACGTc                                                               >  2:2911848/1‑150 (MQ=255)
                                                                                       tCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTc                                                               >  1:470555/1‑150 (MQ=255)
                                                                                           gAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCa                                                           <  1:841120/150‑1 (MQ=255)
                                                                                           gAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCGGCACGGGCAGGTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCa                                                           <  2:5245300/150‑1 (MQ=255)
                                                                                            aaCGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAg                                                          >  2:5202651/1‑150 (MQ=255)
                                                                                               gcgACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGt                                                       >  2:87556/1‑150 (MQ=255)
                                                                                                 gACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCTCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTAt                                                     <  1:712086/150‑1 (MQ=255)
                                                                                                 gACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTAt                                                     >  2:665552/1‑150 (MQ=255)
                                                                                                  aCTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATg                                                    <  1:1012717/150‑1 (MQ=255)
                                                                                                    ttCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATggg                                                  >  1:2425159/1‑150 (MQ=255)
                                                                                                      cGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGTGTGGCTTGCAGCATCAGACCCTGGTGGGGCTGTAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCt                                                <  1:5247486/150‑1 (MQ=255)
                                                                                                      cGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCt                                                <  1:99153/150‑1 (MQ=255)
                                                                                                       gAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTg                                               >  1:4407854/1‑150 (MQ=255)
                                                                                                       gAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTg                                               <  1:4722869/150‑1 (MQ=255)
                                                                                                       gAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTg                                               >  1:595830/1‑150 (MQ=255)
                                                                                                         aTGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGa                                             <  1:3374797/150‑1 (MQ=255)
                                                                                                          tGGAACGACAAAATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGACCGAAGCGTATGGGCTGGat                                            >  2:467658/1‑150 (MQ=255)
                                                                                                           ggAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGata                                           <  1:915776/150‑1 (MQ=255)
                                                                                                             aaCGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATc                                         <  2:2475289/150‑1 (MQ=255)
                                                                                                                  caacaTCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACaa                                    >  2:183510/1‑150 (MQ=255)
                                                                                                                     agccACCCAGCTCGAACTGCCCGTGCAGCTCGGCCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAAccc                                 >  1:676228/4‑150 (MQ=38)
                                                                                                                      aTCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCg                                <  2:5119240/150‑1 (MQ=255)
                                                                                                                        cACCCAGCTCGAACTGCCCGTGCAGCTCGGCCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGt                              >  2:801346/1‑150 (MQ=38)
                                                                                                                         gccgaagCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGta                             >  2:4181994/7‑150 (MQ=255)
                                                                                                                           ccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGtata                           >  1:3150950/1‑150 (MQ=255)
                                                                                                                           ccAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGtata                           >  1:1777750/1‑150 (MQ=255)
                                                                                                                               gCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGc                       >  1:4152783/1‑150 (MQ=255)
                                                                                                                               gCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGc                       <  2:4418858/150‑1 (MQ=255)
                                                                                                                                 cGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCaa                     <  2:1385615/150‑1 (MQ=255)
                                                                                                                                  gAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTACGCAGCg                                                        <  1:791171/114‑1 (MQ=255)
                                                                                                                                       gCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGaa               <  2:4826863/150‑1 (MQ=255)
                                                                                                                                       gCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGaa               <  2:3665364/150‑1 (MQ=255)
                                                                                                                                       gCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGaa               <  2:3411096/150‑1 (MQ=255)
                                                                                                                                       gCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGaa               <  2:3311587/150‑1 (MQ=255)
                                                                                                                                       gCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGaa               <  1:1829824/150‑1 (MQ=255)
                                                                                                                                         aCGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAagc             <  1:2911848/150‑1 (MQ=255)
                                                                                                                                             gCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAagccagc         <  2:3605716/150‑1 (MQ=255)
                                                                                                                                             gCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAagccagc         <  2:3294180/150‑1 (MQ=255)
                                                                                                                                                gTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCAt      <  2:2425159/150‑1 (MQ=255)
                                                                                                                                                 ttCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCACGCCGAAGCCAGCCATc     >  2:4284727/1‑150 (MQ=255)
                                                                                                                                                 ttCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATc     >  1:1083617/1‑150 (MQ=255)
                                                                                                                                                   cGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCAc   >  2:2179215/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTGGCCAACCACTACACCCAGGGCACGCTCAAGGGCAACAGCTCCGAGCCCCAGGCCCTGGTGAGCACCACCATCACGCGCTTCTATCGCGAACGCGACTTCGAATGGAACGACAACATCACCCAGGCCGAGCTGCACGGGCAGTTCGAGCTGGGTGGCTGGCAGCATCAGACCCTGGTGGGGCTGGAGTACGAAAACTACCGCAACAGCCAGAAGTACCCGCAAAGCGCCACGTCCGCAGCGTATGGGCTGGATATCTACAACCCGGTATACGGCAAGCCGAAGCCAGCCATCACC  >  NC_002947/3763538‑3763835

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: