Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 451,273 A→T 6.5% R201R (CGT→CGA PP_0371 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947451,2730AT 6.5% 271.2 / 6.2 89R201R (CGT→CGAPP_0371LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (35/48);  new base T (3/3);  total (38/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTG  >  NC_002947/451124‑451422
                                                                                                                                                      |                                                                                                                                                     
tcatcTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGa                                                                                                                                                       <  1:1054006/150‑1 (MQ=255)
tcatcTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGa                                                                                                                                                       <  2:2120135/150‑1 (MQ=255)
tcatcTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGa                                                                                                                                                       <  1:1461232/150‑1 (MQ=255)
tcatcTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGa                                                                                                                                                       <  1:1305938/150‑1 (MQ=255)
 catcTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGAc                                                                                                                                                      <  1:183355/150‑1 (MQ=255)
  atcTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACg                                                                                                                                                     >  2:2334252/1‑150 (MQ=255)
       cggcggCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCt                                                                                                                                                <  2:1028084/150‑1 (MQ=255)
         gcggcCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCAACAACAGCCTTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTgg                                                                                                                                              <  2:593404/150‑1 (MQ=255)
               tcctTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAActcc                                                                                                                                        >  2:1322466/1‑150 (MQ=255)
               tcctTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAActcc                                                                                                                                        <  1:1322466/150‑1 (MQ=255)
                 ctTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAActcctc                                                                                                                                      >  1:124413/1‑150 (MQ=255)
                      cgcgcCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCca                                                                                                                                 <  1:486546/150‑1 (MQ=255)
                       gcgcCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAActcctc                                                                                                                                      >  1:1034163/1‑144 (MQ=255)
                       gcgcCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAActcctc                                                                                                                                      <  2:1034163/144‑1 (MQ=255)
                         gcCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCcagca                                                                                                                              <  2:2383700/150‑1 (MQ=255)
                            ccGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCCGTCAGATCTggca                                                                                                                                            <  1:1008083/133‑3 (MQ=255)
                            ccGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCCGTCAGATCTggca                                                                                                                                            >  2:1008083/1‑131 (MQ=255)
                             cGGCAGAATGGTCGCCAGCGCTGTATTGATGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCcatgctgat                                                                                                                          >  1:1271346/1‑143 (MQ=255)
                               gCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTgg                                                                                                                                              <  2:1623430/128‑1 (MQ=255)
                               gCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTgg                                                                                                                                              >  1:1623430/1‑128 (MQ=255)
                                cAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCcagcatcagcat                                                                                                                       <  2:1147148/150‑1 (MQ=255)
                                 aGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATGCCTGCCAGATTGTGCATCTATGCCTGCAGCTGCGGCGGTACGGGCCTGGTTGCCTGTGGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCGGCATCAGCATg                                                                                                                      >  1:550341/1‑150 (MQ=255)
                                  gAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCCGTCAGATCTggcag                                                                                                                                           >  1:387213/1‑125 (MQ=255)
                                  gAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCCGTCAGATCTggcag                                                                                                                                           >  1:386377/1‑125 (MQ=255)
                                  gAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCCGTCAGATCTggcag                                                                                                                                           <  2:386377/128‑4 (MQ=255)
                                  gAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCCGTCAGATCTggcag                                                                                                                                           <  2:387213/128‑4 (MQ=255)
                                   aaTGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATggg                                                                                                                    <  1:586863/150‑1 (MQ=255)
                                    aTGGTCGCCAGCGTTTTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGc                                                                                                                   <  2:2244346/150‑1 (MQ=255)
                                      ggTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCca                                                                                                                                 >  2:1101649/1‑134 (MQ=255)
                                      ggTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCca                                                                                                                                 <  1:1101649/134‑1 (MQ=255)
                                            cAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGACAACTctgc                                                                                                           >  1:714655/1‑150 (MQ=255)
                                                    tgtgCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGAc                                                                                                   >  2:116560/1‑150 (MQ=255)
                                                       gcgcGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGa                                                                                                >  1:2095789/1‑150 (MQ=255)
                                                        cgcgAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGAt                                                                                               <  1:1117480/150‑1 (MQ=255)
                                                         gcgAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATg                                                                                              <  2:733589/150‑1 (MQ=255)
                                                           gAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCcagcatc                                                                                                                            <  2:401672/118‑1 (MQ=255)
                                                           gAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCcagcatc                                                                                                                            >  1:401672/1‑118 (MQ=255)
                                                            aGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGac                                                                                           <  2:925362/150‑1 (MQ=255)
                                                             ggCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGaca                                                                                          <  2:2095789/150‑1 (MQ=255)
                                                                 aTCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACctt                                                                                      >  2:1337416/1‑150 (MQ=255)
                                                                  tCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGa                                                                                            >  2:610396/1‑143 (MQ=255)
                                                                  tCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGa                                                                                            <  1:610396/143‑1 (MQ=255)
                                                                    gagaATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACcttcta                                                                                   <  2:1656950/150‑2 (MQ=255)
                                                                     agaATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTg                                                                                  >  2:2128398/1‑150 (MQ=255)
                                                                       aaTCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGaa                                                                                <  1:961178/150‑1 (MQ=255)
                                                                       aaTCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGaa                                                                                <  1:37461/150‑1 (MQ=255)
                                                                        aTCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAg                                                                               <  2:1886637/150‑1 (MQ=255)
                                                                         tCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACcttc                                                                                     <  2:1203779/143‑1 (MQ=255)
                                                                         tCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACcttc                                                                                     >  1:1203779/1‑143 (MQ=255)
                                                                         tCCCTGCCATATTGCTCCTCTATGCCTGCACCTGCGGCCGTCAGCCCTGACTGGTCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGt                                                                              <  2:127175/150‑1 (MQ=255)
                                                                               ccATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCa                                                                        >  1:918267/1‑150 (MQ=255)
                                                                                cATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAg                                                                       >  1:346621/1‑150 (MQ=255)
                                                                                 atatTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGc                                                                      <  1:31357/150‑1 (MQ=255)
                                                                                     tGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCTGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCgggt                                                                  <  2:714655/150‑1 (MQ=255)
                                                                                     tGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCgggt                                                                  <  1:2313390/150‑1 (MQ=255)
                                                                                     tGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGACTTCTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCgggt                                                                  <  1:1869695/150‑1 (MQ=255)
                                                                                       ttCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCgggtgg                                                                <  2:124413/150‑1 (MQ=255)
                                                                                            ctctGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATg                                                                                              >  1:1541182/1‑115 (MQ=255)
                                                                                            ctctGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATg                                                                                              <  2:1541182/115‑1 (MQ=255)
                                                                                             tctGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCg                                                          >  1:777596/1‑150 (MQ=255)
                                                                                              ctgcctgcACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCgg                                                         <  1:2183219/150‑1 (MQ=255)
                                                                                                 cctgcACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTa                                                      <  1:903897/150‑1 (MQ=255)
                                                                                                  ctgcACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGATTGCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTAt                                                     <  1:287129/150‑1 (MQ=255)
                                                                                                   tgcACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGacac                                                                                         <  2:1067744/113‑1 (MQ=255)
                                                                                                   tgcACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGacac                                                                                         >  1:1067744/1‑113 (MQ=255)
                                                                                                   tgcACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAg                                                                               <  2:1139645/123‑1 (MQ=255)
                                                                                                   tgcACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAg                                                                               >  1:1139645/1‑123 (MQ=255)
                                                                                                     cACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGAc                                                  >  1:711321/1‑150 (MQ=255)
                                                                                                     cACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGAc                                                  <  2:918267/150‑1 (MQ=255)
                                                                                                      aCCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGAcc                                                 >  1:499771/1‑150 (MQ=255)
                                                                                                        cTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGcc                                                                                                       >  2:1241343/1‑94 (MQ=255)
                                                                                                        cTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGcc                                                                                                       <  1:1241343/94‑1 (MQ=255)
                                                                                                         tGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGTACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGt                                              >  2:711087/1‑150 (MQ=255)
                                                                                                              ccGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAg                                         <  1:2334252/150‑1 (MQ=255)
                                                                                                                   cGGCCCTGCCTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCa                                     <  2:630469/150‑1 (MQ=255)
                                                                                                                               gCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCcagcatca                                                                                                                           <  2:941557/52‑1 (MQ=255)
                                                                                                                               gCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCcagcatca                                                                                                                           >  1:941557/1‑52 (MQ=255)
                                                                                                                               gCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCggcggc                                                                                                        >  1:1472132/1‑71 (MQ=255)
                                                                                                                               gCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCggcggc                                                                                                        <  2:1472132/71‑1 (MQ=255)
                                                                                                                                      gcgcCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGc                  <  1:996120/150‑1 (MQ=255)
                                                                                                                                        gcCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGAt                                                                                               >  1:1100696/1‑71 (MQ=255)
                                                                                                                                        gcCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGAt                                                                                               <  2:1100696/71‑1 (MQ=255)
                                                                                                                                        gcCTGCCAGATCTGACGGAACTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCAGCTTCGCACAACAGCAAGCCa                >  2:897805/1‑150 (MQ=255)
                                                                                                                                            gCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCa                                     <  1:323947/125‑1 (MQ=255)
                                                                                                                                            gCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCa                                     >  2:323947/1‑125 (MQ=255)
                                                                                                                                             ccAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAg           >  1:1265894/1‑150 (MQ=255)
                                                                                                                                             ccAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAg           >  1:839659/1‑150 (MQ=255)
                                                                                                                                               aGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCgggtgg                                                                <  2:1605155/95‑1 (MQ=255)
                                                                                                                                               aGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCgggtgg                                                                >  1:1605155/1‑95 (MQ=255)
                                                                                                                                               aGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGcc         >  1:823429/1‑150 (MQ=255)
                                                                                                                                                 aTCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCct       <  2:346621/150‑1 (MQ=255)
                                                                                                                                                  tCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGCGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCctg      <  1:1337416/150‑1 (MQ=255)
                                                                                                                                                   cTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCctgg     <  2:651147/150‑1 (MQ=255)
                                                                                                                                                      aCGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGcc         >  2:1554997/1‑143 (MQ=255)
                                                                                                                                                      aCGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGcc         <  1:1554997/143‑1 (MQ=255)
                                                                                                                                                      aCGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCctggctg  <  1:711087/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCATCTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTG  >  NC_002947/451124‑451422

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: