Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,406,104 C→G 7.4% A184A (GCG→GCC PP_4750 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,406,1040CG 7.4% 300.9 / 9.8 94A184A (GCG→GCCPP_4750amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (38/49);  new base G (4/3);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATA  >  NC_002947/5405959‑5406248
                                                                                                                                                 |                                                                                                                                                 
acaAAGTCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCgg                                                                                                                                               <  1:723972/150‑1 (MQ=255)
 caAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                              >  2:840903/1‑150 (MQ=255)
 caAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                              <  1:1411459/150‑1 (MQ=255)
  aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCCCTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgtg                                                                                                                                             >  2:313884/1‑150 (MQ=255)
  aaaGCACTTCCACCACCCCGACAACCATTGCATCCACCTCGCCAGCCGTTCGCTCATCGGGCAATCCTTTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGCCCACTGGCACGCTGAACCCCAGGTACACCACCTCCGCCGCGatg                                                                                                                                             <  2:151524/150‑1 (MQ=255)
       ccTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGAc                                                                                                                                        <  2:314180/150‑1 (MQ=255)
             acaACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAg                                                                                                                                 <  1:2086293/150‑1 (MQ=255)
             acaACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGGTTCAACGTGGCCGCCAGCCACCGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGGGTAGACGCTCAg                                                                                                                                 >  1:549071/1‑150 (MQ=255)
                 cTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGccc                                                                                                                             >  2:760891/1‑150 (MQ=255)
                 cTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGccc                                                                                                                             <  2:513842/150‑1 (MQ=255)
                       cAACCATTGCATCCACATCGTCAGCAGTTCGCTCCTCGGGCAATCCTCTGTTTCCACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCATGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAg                                                                                                                       <  1:209980/150‑1 (MQ=255)
                       cAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAg                                                                                                                       <  2:1059685/150‑1 (MQ=255)
                        aaCCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAgg                                                                                                                      <  2:851104/150‑1 (MQ=255)
                         aCCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGt                                                                                                                     <  2:1499643/150‑1 (MQ=255)
                         aCCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGt                                                                                                                     <  1:840903/150‑1 (MQ=255)
                          ccATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGtt                                                                                                                    <  1:2257604/150‑1 (MQ=255)
                           cATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGttt                                                                                                                   <  1:1013488/150‑1 (MQ=255)
                             ttGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACTCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGc                                                                                                                 >  1:2079538/1‑150 (MQ=255)
                             ttGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTg                                                                                                                  >  2:1632450/1‑149 (MQ=255)
                             ttGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTg                                                                                                                  <  1:1632450/149‑1 (MQ=255)
                             ttGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGc                                                                                                                 <  1:1050193/150‑1 (MQ=255)
                                        tCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGATCGTCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTg                                                                                                      <  1:150054/150‑1 (MQ=255)
                                        tCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTg                                                                                                      <  2:585984/150‑1 (MQ=255)
                                        tCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTg                                                                                                      <  2:1455778/150‑1 (MQ=255)
                                          gTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGt                                                                                                    >  1:1418217/1‑150 (MQ=255)
                                           tCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACAACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTa                                                                                                   >  1:1196690/1‑150 (MQ=255)
                                             agcagTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggt                                                                                                                                           <  2:1060415/109‑3 (MQ=255)
                                             agcagTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggt                                                                                                                                           >  1:1060415/1‑107 (MQ=255)
                                              gcagTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCcgcg                                                                                                                                                <  1:991618/103‑1 (MQ=255)
                                              gcagTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCcgcg                                                                                                                                                >  2:991618/1‑103 (MQ=255)
                                                    cGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGaaa                                                                                          <  2:1320132/150‑1 (MQ=255)
                                                     gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAg                                                                                         >  1:2170304/1‑150 (MQ=255)
                                                      cTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGc                                                                                        >  1:1416118/1‑150 (MQ=255)
                                                          tCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCAccgccgccg                                                                                                                                                  >  2:390219/1‑89 (MQ=255)
                                                          tCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCAccgccgccg                                                                                                                                                  <  1:390219/89‑1 (MQ=255)
                                                                aaTCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCAcc                                                                              >  1:708731/1‑150 (MQ=255)
                                                                  tCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGc                                                                                                                                      >  2:786294/1‑93 (MQ=255)
                                                                  tCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGc                                                                                                                                      <  1:786294/93‑1 (MQ=255)
                                                                     tctGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGc                                                                         >  1:1003826/1‑150 (MQ=255)
                                                                         tttCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGAt                                                                                    <  1:478088/135‑1 (MQ=255)
                                                                         tttCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGAt                                                                                    >  2:478088/1‑135 (MQ=255)
                                                                               cGTGTCCGACAGCCAACGAGAGCACAAGCCTACTGGCACGCTGAGCACCAGGTACACCACCGCCGTCGCGGTGTAGCCGCTCAGGCCCTTAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGa                                                               <  2:1149135/150‑1 (MQ=255)
                                                                               cGTGGCCGACAGCCAACGACATAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGa                                                               <  1:126091/150‑1 (MQ=255)
                                                                                      gACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGaaa                                                                                                                        >  1:1084391/1‑86 (MQ=255)
                                                                                      gACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGaaa                                                                                                                        <  2:1084391/86‑1 (MQ=255)
                                                                                      gACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACg                                                        <  1:1001783/150‑1 (MQ=255)
                                                                                      gACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACg                                                        >  1:1126887/1‑150 (MQ=255)
                                                                                      gACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACg                                                        >  1:1137353/1‑150 (MQ=255)
                                                                                           ccAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTATACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGa                                                   >  2:2149567/1‑150 (MQ=255)
                                                                                                gagagaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCt                                              <  2:708731/150‑1 (MQ=255)
                                                                                                  gagaACAAGGCCACTTTCACGCGGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGCCTGCCGAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTg                                            <  1:760891/150‑1 (MQ=255)
                                                                                                  gagaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGc                                                                         >  2:368019/1‑121 (MQ=255)
                                                                                                  gagaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGc                                                                         <  1:368019/121‑1 (MQ=255)
                                                                                                      aCAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCa                                        <  2:1003826/150‑1 (MQ=255)
                                                                                                       cAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTgg                                                                                                     >  1:2076063/1‑88 (MQ=255)
                                                                                                       cAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTgg                                                                                                     <  2:2076063/88‑1 (MQ=255)
                                                                                                       cAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTgg                                                                                                     <  2:2075489/88‑1 (MQ=255)
                                                                                                       cAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTgg                                                                                                     >  1:2075489/1‑88 (MQ=255)
                                                                                                       cAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGACGCCACACTATATCGTCTGATGTCCTTAAAGGTTTGCTGGCAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAg                                       <  1:543202/150‑1 (MQ=255)
                                                                                                        aaGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAgg                                      >  2:973951/1‑150 (MQ=255)
                                                                                                        aaGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAgg                                      <  2:1137353/150‑1 (MQ=255)
                                                                                                        aaGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAgg                                      >  2:1724587/1‑150 (MQ=255)
                                                                                                          ggCCACTGGCACGCTGAGCACCGGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAAtcgt                                                                  <  2:1114134/120‑1 (MQ=255)
                                                                                                          ggCCACTGGCACGCTGAGCACCGGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAAtcgt                                                                  >  1:1114134/1‑120 (MQ=255)
                                                                                                          ggCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGcc                                    >  2:589769/1‑150 (MQ=255)
                                                                                                           gCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCt                                              >  2:1043492/1‑139 (MQ=255)
                                                                                                           gCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCt                                              <  1:1043492/139‑1 (MQ=255)
                                                                                                           gCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCa                                   >  1:1570392/1‑150 (MQ=255)
                                                                                                              ccTGGCACGATGAGCACCAGGTACACCGCTGCCGCCGCGGTGTAGACGGTCAGGCCCTGAAAGGTTTGCTGTGAACCCTGGTAACCCTGAAAGCGGATTTCACCAACGCCACTCGGCGAGGTCACGGCCTAGTCCTTGAGAAGGCCAATg                                <  2:549071/149‑1 (MQ=255)
                                                                                                               cTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGa                                                   >  2:259571/1‑130 (MQ=255)
                                                                                                               cTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGa                                                   <  1:259571/130‑1 (MQ=255)
                                                                                                                  gCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGcc                                                     >  1:414402/1‑125 (MQ=255)
                                                                                                                  gCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGcc                                                     <  2:414402/125‑1 (MQ=255)
                                                                                                                  gCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGt                            >  1:2242035/1‑150 (MQ=255)
                                                                                                                      gCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGt                        >  1:438965/1‑150 (MQ=255)
                                                                                                                        tGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTgg                      >  2:19141/1‑150 (MQ=255)
                                                                                                                          aGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCt                    <  1:589769/150‑1 (MQ=255)
                                                                                                                          aGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCt                    <  2:1418217/150‑1 (MQ=255)
                                                                                                                             aCCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAg                 <  2:1416118/150‑1 (MQ=255)
                                                                                                                               cAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGt               <  2:2242035/150‑1 (MQ=255)
                                                                                                                               cAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCAACGGCCGAGTCTTTTAGCAGGCCAATGGCGTAGGTGGCTGAGGt               >  2:2200461/1‑150 (MQ=255)
                                                                                                                                aGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTg              <  1:313884/150‑1 (MQ=255)
                                                                                                                                 ggTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCCACGCCAAGGTGCGAGTCGGCGTTCGAGTCCTTTAGCAGGCCAATGGCGTAGGTGGCTGAGGTgg             <  1:658034/150‑1 (MQ=255)
                                                                                                                                 ggTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTgg             <  1:2200461/150‑1 (MQ=255)
                                                                                                                                    acacCACCGCCGCCGCGGTGTAGCACTCTAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGaa          <  1:1209839/150‑1 (MQ=255)
                                                                                                                                     ctacaccGCGGCGGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAAt         >  1:1154947/4‑150 (MQ=255)
                                                                                                                                      tacaccGCGGCGGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCa                                        >  2:1492131/3‑118 (MQ=255)
                                                                                                                                      tacaccGCGGCGGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCa                                        <  1:1492131/116‑1 (MQ=255)
                                                                                                                                      tacaccGCGGCGGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAg                 >  2:1710944/3‑141 (MQ=255)
                                                                                                                                      tacaccGCGGCGGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAg                 <  1:1710944/139‑1 (MQ=255)
                                                                                                                                       ccaccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg       <  2:1894443/150‑1 (MQ=255)
                                                                                                                                       ccaccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg       >  1:1894443/1‑150 (MQ=255)
                                                                                                                                        caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAgg                                      >  2:3819/1‑118 (MQ=255)
                                                                                                                                        caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAgg                                      <  1:3819/118‑1 (MQ=255)
                                                                                                                                         accGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCtt                                             >  2:1718161/1‑110 (MQ=255)
                                                                                                                                         accGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCtt                                             <  1:1718161/110‑1 (MQ=255)
                                                                                                                                         accGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCg     >  2:291494/1‑150 (MQ=255)
                                                                                                                                         accGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCg     >  2:1152414/1‑150 (MQ=255)
                                                                                                                                          ccgccgccgCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGa    <  1:1152414/150‑1 (MQ=255)
                                                                                                                                            gccgccgcGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGATGGGAATGGCGATa  >  2:472310/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATA  >  NC_002947/5405959‑5406248

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