Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,187,836 A→C 5.1% D110A (GAC→GCC)  PP_4562 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,187,8360AC 5.1% 1007.3 / 30.2 293D110A (GAC→GCC) PP_4562hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (130/148);  new base C (10/5);  total (140/153)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 302 total aligned reads displayed.

CACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACG  >  NC_002947/5187687‑5187984
                                                                                                                                                     |                                                                                                                                                       
cacGCCGGAAAGCGACGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGa                                                                                                                                                         <  2:1960909/150‑1 (MQ=255)
 acGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCACGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGac                                                                                                                                                        >  2:8565281/1‑150 (MQ=255)
   gCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGacac                                                                                                                                                      <  1:3951945/150‑1 (MQ=255)
     cGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACAccc                                                                                                                                                    <  1:2123150/150‑1 (MQ=255)
        aaaGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTc                                                                                                                                                 <  1:6467733/150‑1 (MQ=255)
        aaaGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCAGCTGCTCGCCACCCgga                                                                                                                                                 <  2:4538876/150‑3 (MQ=255)
           gCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCgg                                                                                                                                              <  2:6101206/150‑1 (MQ=255)
           gCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCgg                                                                                                                                              <  1:77451/150‑1 (MQ=255)
             gATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGt                                                                                                                                            <  1:6765342/150‑1 (MQ=255)
                gAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑Gc                                                                                                                                       <  1:3325734/150‑1 (MQ=255)
                gAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGGACGGGTgtc                                                                                                                                         <  1:297174/150‑3 (MQ=255)
                 aGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGCCACCCGGACGGGTG‑‑TCg                                                                                                                                      <  2:5746384/150‑1 (MQ=255)
                  gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑a                                                                                                                                    >  1:2499615/1‑150 (MQ=255)
                  gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑a                                                                                                                                    >  1:7243362/1‑150 (MQ=255)
                  gCTGCAGCACCTGAAGGTGCCCACCAATTTCCACAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCAGGACCTCAGCGCCGACCGGGACTCGCTGCTCCGCGACACCCGTACGGGTG‑‑GCG‑a                                                                                                                                    >  1:8648839/1‑150 (MQ=255)
                   cTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑Gc                                                                                                                                       <  1:1874357/147‑1 (MQ=255)
                     gcagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑agca                                                                                                                                 <  1:6043527/150‑1 (MQ=255)
                       agcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑agcagc                                                                                                                               <  1:22034/150‑1 (MQ=255)
                        gcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGAAGTC‑CTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTg                                                                                                                             <  2:4160170/150‑1 (MQ=255)
                            cTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑ATCAGCTGGAc                                                                                                                          <  2:7869153/150‑1 (MQ=255)
                            cTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGAc                                                                                                                          <  1:1566165/150‑1 (MQ=255)
                              gAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACtt                                                                                                                        >  1:1109835/1‑150 (MQ=255)
                                aGGTGCCCACCAAGTTCCTCACGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGCCCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTcc                                                                                                                      <  1:5060286/150‑1 (MQ=255)
                                  gTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCACGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGc                                                                                                                    >  2:5385166/1‑150 (MQ=255)
                                  gTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGc                                                                                                                    >  2:4065670/1‑150 (MQ=255)
                                     cccACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTTGCAAGCGCCGCTCCTCCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAAc                                                                                                                 <  1:4491871/150‑1 (MQ=255)
                                       caccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACtt                                                                                                               >  2:4259588/1‑150 (MQ=255)
                                        accaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTc                                                                                                              >  2:846014/1‑150 (MQ=255)
                                        accaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTc                                                                                                              >  2:8036629/1‑150 (MQ=255)
                                        accaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTc                                                                                                              <  2:1283751/150‑1 (MQ=255)
                                        accaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTc                                                                                                              <  1:1359058/150‑1 (MQ=255)
                                         ccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCg                                                                                                             <  2:7196387/150‑1 (MQ=255)
                                          ccaTTTGCTCAAGGTGAGCCTTGAGGAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGt                                                                                                            >  1:2642257/3‑150 (MQ=255)
                                          caATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGt                                                                                                            >  1:101550/1‑150 (MQ=255)
                                                ccTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGc                                                                                                      >  2:6953523/1‑150 (MQ=255)
                                                 cTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCg                                                                                                     >  1:5246956/1‑150 (MQ=255)
                                                  tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCACGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                    >  1:7050830/1‑150 (MQ=255)
                                                  tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                    >  2:2412403/1‑150 (MQ=255)
                                                  tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                    >  2:7046807/1‑150 (MQ=255)
                                                  tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                    >  1:8398665/1‑150 (MQ=255)
                                                  tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                    >  1:399215/1‑150 (MQ=255)
                                                  tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                    >  1:5369879/1‑150 (MQ=255)
                                                     aggtgagccTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAg                                                                                                 >  2:4364316/8‑150 (MQ=255)
                                                     aGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAg                                                                                                 <  2:3084580/150‑1 (MQ=255)
                                                      ggtgagccTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTTCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAgg                                                                                                >  2:1724608/7‑150 (MQ=255)
                                                           aCCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTc                                                                                           <  1:5918754/150‑1 (MQ=255)
                                                              ttGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAgg                                                                                        >  2:5505663/1‑150 (MQ=255)
                                                              ttGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAgg                                                                                        <  2:1058156/150‑1 (MQ=255)
                                                              gtGAGCACATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAgg                                                                                        <  2:510312/149‑1 (MQ=255)
                                                                  gCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCttt                                                                                    <  1:2576761/150‑1 (MQ=255)
                                                                  gCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCttt                                                                                    <  2:4216016/150‑1 (MQ=255)
                                                                   cAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtc                                                                                   <  1:4217428/150‑1 (MQ=255)
                                                                     aaTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtctc                                                                                 <  1:5472053/150‑1 (MQ=255)
                                                                      aTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCg                                                                                >  2:3917257/1‑150 (MQ=255)
                                                                       tGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCAGCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACGCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACTGGCGTAAGGCCGTCTGGCTTTCTCGa                                                                               <  2:5279653/150‑1 (MQ=255)
                                                                        ggAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                              <  1:7878040/150‑1 (MQ=255)
                                                                          aaGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGc                                                                            >  1:5694298/1‑150 (MQ=255)
                                                                            gTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGa                                                                          <  1:985266/150‑1 (MQ=255)
                                                                             tGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGAt                                                                         <  1:7591319/150‑1 (MQ=255)
                                                                                     ccTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATgg                                                                 <  1:1801970/150‑1 (MQ=255)
                                                                                            aGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGa                                                          <  1:2345667/150‑1 (MQ=255)
                                                                                             gcgcCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGAt                                                         >  2:2798578/1‑150 (MQ=255)
                                                                                             gcgcCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGAt                                                         <  1:5732045/150‑1 (MQ=255)
                                                                                             gcgcCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAACGCCGTCAGGCTTTCTCGATGCGATCGTCATCCATGACGAt                                                         <  2:6209226/150‑1 (MQ=255)
                                                                                               gccgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCa                                                       >  2:646381/1‑150 (MQ=255)
                                                                                                     ccaccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCt                                                 <  2:1112074/150‑1 (MQ=255)
                                                                                                      caccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTg                                                >  2:4857400/1‑150 (MQ=255)
                                                                                                       accaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGc                                               >  2:2591928/1‑150 (MQ=255)
                                                                                                       accaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACGTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGc                                               <  2:7381616/150‑1 (MQ=255)
                                                                                                        ccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCt                                              <  1:8703173/150‑1 (MQ=255)
                                                                                                        ccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCt                                              <  1:3129847/150‑1 (MQ=255)
                                                                                                        ccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAAGcagggcc                                              <  2:4657567/150‑7 (MQ=255)
                                                                                                           tCAGCCGGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGa                                           <  2:2902710/150‑1 (MQ=255)
                                                                                                            cAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaa                                          <  1:4065014/150‑1 (MQ=255)
                                                                                                            cAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaa                                          <  2:6801502/150‑1 (MQ=255)
                                                                                                                cTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaacaac                                      >  2:1753413/1‑150 (MQ=255)
                                                                                                                 tGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAAcc                                     <  1:472244/150‑1 (MQ=255)
                                                                                                                 tGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAAcc                                     >  1:6389116/1‑150 (MQ=255)
                                                                                                                    cccTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGc                                  <  2:35541/150‑1 (MQ=255)
                                                                                                                     ccTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGcc                                 >  1:6601021/1‑150 (MQ=255)
                                                                                                                       tCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCaa                               >  1:3707336/1‑150 (MQ=255)
                                                                                                                           cgcCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGc                           <  2:7682592/150‑1 (MQ=255)
                                                                                                                              cGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCctt                        <  1:8036629/150‑1 (MQ=255)
                                                                                                                                 ccGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCcttctt                     <  2:3283709/150‑1 (MQ=255)
                                                                                                                                    ggACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGaa                  <  1:4005416/150‑1 (MQ=255)
                                                                                                                                    ggACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGaa                  >  2:1183405/1‑150 (MQ=255)
                                                                                                                                    ggACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGaa                  <  1:2046497/150‑1 (MQ=255)
                                                                                                                                     gACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAg                 <  2:7262808/150‑1 (MQ=255)
                                                                                                                                       cTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTTAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGTCGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGtt               <  2:5822636/150‑1 (MQ=255)
                                                                                                                                        tCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTg              >  1:93873/1‑150 (MQ=255)
                                                                                                                                        gtccagctgctCGCCACCCGGACGGGTG‑‑TCG‑CGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCGTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTg              >  1:8797849/6‑150 (MQ=38)
                                                                                                                                         cGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc             <  2:1504490/150‑1 (MQ=255)
                                                                                                                                          gctgctGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGcc            >  1:3609129/1‑150 (MQ=255)
                                                                                                                                            agctgctCGCCACCCGGACGGGTGTCGCGCAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCGTGCTGGAAGGTGc             >  2:7947984/2‑150 (MQ=37)
                                                                                                                                            agctgctCGCCACCCGGACGGGTG‑‑TCG‑CGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCt          >  1:1141798/2‑150 (MQ=255)
                                                                                                                                              ctgcGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTg        <  1:4882527/150‑1 (MQ=255)
                                                                                                                                               tgcTCGCCACCCGGACGGGTG‑‑TCG‑CGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGc       >  1:6439473/1‑150 (MQ=255)
                                                                                                                                               tgcGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGc       <  2:3629716/150‑1 (MQ=255)
                                                                                                                                                   cgACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGAc   >  2:6283829/1‑150 (MQ=255)
                                                                                                                                                   cgACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGAc   <  2:1568362/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTG‑‑GCG‑AGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACG  >  NC_002947/5187687‑5187984

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: ‑