Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,654,984 C→A 5.8% G309C (GGC→TGC)  PP_4113 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,654,9840CA 5.8% 764.5 / 24.0 227G309C (GGC→TGC) PP_4113hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (97/116);  new base A (5/8);  total (102/125)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 237 total aligned reads displayed.

ATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTCAA  >  NC_002947/4654835‑4655133
                                                                                                                                                     |                                                                                                                                                     
aTCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGcc                                                                                                                                                       >  1:3801580/1‑150 (MQ=255)
      cgcgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGGGCAAGCACAGccc                                                                                                                                                 <  2:2542923/150‑1 (MQ=255)
       gcgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGcccc                                                                                                                                                <  2:1204011/150‑1 (MQ=255)
       gcgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGGGCAAGCACAGccct                                                                                                                                                <  1:5141996/150‑2 (MQ=255)
          ccGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGc                                                                                                                                             <  1:1271356/150‑1 (MQ=255)
          ccGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGc                                                                                                                                             <  2:3505381/150‑1 (MQ=255)
               ccGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCaaaaa                                                                                                                                        >  1:3385294/1‑150 (MQ=255)
               ccGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCaaaaa                                                                                                                                        >  1:3725133/1‑150 (MQ=255)
                cGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCGGGGCTGGGCAAGCACAGCCCTGGcttcgag                                                                                                                                       <  2:1553945/150‑7 (MQ=255)
                  gTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGcc                                                                                                                                     <  1:1010745/150‑1 (MQ=255)
                  gTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGcc                                                                                                                                     <  1:2805788/150‑1 (MQ=255)
                   tGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCGCCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGccc                                                                                                                                    <  1:4976382/150‑1 (MQ=255)
                   tGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGccc                                                                                                                                    <  1:4658071/150‑1 (MQ=255)
                    gAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCt                                                                                                                                   >  1:2987026/1‑150 (MQ=255)
                          cccAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCt                                                                                                                             >  2:2837746/1‑150 (MQ=255)
                           ccAGCCCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTg                                                                                                                            >  2:391767/1‑150 (MQ=255)
                            cAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTgc                                                                                                                           <  1:2321585/150‑1 (MQ=255)
                                tCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAg                                                                                                                       >  2:3970969/1‑150 (MQ=255)
                                   aGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCga                                                                                                                    <  1:907623/150‑1 (MQ=255)
                                   aGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCga                                                                                                                    >  1:540333/1‑150 (MQ=255)
                                    gTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCgag                                                                                                                   >  1:113862/1‑150 (MQ=255)
                                         gCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGGTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGc                                                                                                              >  2:3209043/1‑150 (MQ=255)
                                           aCCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCt                                                                                                            <  1:2641886/150‑1 (MQ=255)
                                              cGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCg                                                                                                         <  1:4436497/150‑1 (MQ=255)
                                              cGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCg                                                                                                         >  1:917812/1‑150 (MQ=255)
                                               gCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGa                                                                                                        <  2:1719629/150‑1 (MQ=255)
                                                  cGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACAt                                                                                                     <  2:2573078/150‑1 (MQ=255)
                                                   gCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATg                                                                                                    <  1:995553/150‑1 (MQ=255)
                                                   gCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATg                                                                                                    <  2:2378327/150‑1 (MQ=255)
                                                   gCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATg                                                                                                    <  2:2378493/150‑1 (MQ=255)
                                                     cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGc                                                                                                  >  2:5460707/1‑150 (MQ=255)
                                                     cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGc                                                                                                  >  1:3631812/1‑150 (MQ=255)
                                                     cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCACAGCGAGAATGCCTGCGACATGGc                                                                                                  >  2:4496271/1‑150 (MQ=255)
                                                             aCCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCGCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTgg                                                                                          <  1:190570/150‑1 (MQ=255)
                                                             aCCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTgg                                                                                          >  1:2947485/1‑150 (MQ=255)
                                                                aGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCtt                                                                                       <  1:5594016/150‑1 (MQ=255)
                                                                  gTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGc                                                                                     <  2:4037695/150‑1 (MQ=255)
                                                                  gTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGc                                                                                     <  2:5516715/150‑1 (MQ=255)
                                                                       ccGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGc                                                                                <  1:4491856/150‑1 (MQ=255)
                                                                        cGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCGGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                               <  2:4322644/150‑1 (MQ=255)
                                                                        cGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                               <  1:622436/150‑1 (MQ=255)
                                                                        cGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCg                                                                               <  1:330471/150‑1 (MQ=255)
                                                                            gCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTa                                                                           >  2:2364413/1‑150 (MQ=255)
                                                                               gCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCg                                                                        >  2:757080/1‑150 (MQ=255)
                                                                                 gATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTc                                                                      >  1:4667632/1‑150 (MQ=255)
                                                                                 gATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTc                                                                      >  2:1296793/1‑150 (MQ=255)
                                                                                   tGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCGACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGTTTACCGTcgc                                                                    <  2:5421588/150‑1 (MQ=255)
                                                                                     ttGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGa                                                                  <  1:4076495/150‑1 (MQ=255)
                                                                                     ttGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGa                                                                  <  2:2196871/150‑1 (MQ=255)
                                                                                       ggCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCAAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGAGGCTTGCCGGGCGGTTACCGTCGCGAGc                                                                <  2:1305555/150‑1 (MQ=255)
                                                                                          ggTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGACTGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTa                                                             <  2:187852/150‑1 (MQ=255)
                                                                                            tGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGTGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTaga                                                           <  2:5273798/150‑1 (MQ=255)
                                                                                             gATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTagag                                                          <  2:3163416/150‑1 (MQ=255)
                                                                                               tACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                        <  1:1555036/150‑1 (MQ=255)
                                                                                                 ccGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCAt                                                      <  1:1241657/150‑1 (MQ=255)
                                                                                                 ccGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCAt                                                      <  2:4752846/150‑1 (MQ=255)
                                                                                                  ctggcttcGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATg                                                     >  2:4841016/8‑150 (MQ=255)
                                                                                                  cGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATg                                                     >  1:5191482/1‑150 (MQ=255)
                                                                                                  cGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATg                                                     >  1:2765878/1‑150 (MQ=255)
                                                                                                  cGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCAGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATg                                                     >  2:5586108/1‑150 (MQ=255)
                                                                                                  cGGCATCCGAGTTGCGCTCTACCATCACCGCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTCCCGGCGCGAGCGTAGAGGCATg                                                     >  2:1652953/1‑150 (MQ=255)
                                                                                                   ggCATCCGAGATGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGt                                                    >  2:2669751/1‑150 (MQ=255)
                                                                                                      aTCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTg                                                 <  1:3048137/150‑1 (MQ=255)
                                                                                                        ccGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGc                                               <  1:1630934/150‑1 (MQ=255)
                                                                                                          gAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCt                                             >  1:1473462/1‑150 (MQ=255)
                                                                                                           aGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCtt                                            <  1:2278936/150‑1 (MQ=255)
                                                                                                             ttGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGt                                          >  2:2058367/1‑150 (MQ=255)
                                                                                                                 gcTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGc                                      >  2:1522540/1‑150 (MQ=255)
                                                                                                                  ctctACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCa                                     <  2:306721/150‑1 (MQ=255)
                                                                                                                  ctctACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCa                                     <  1:367725/150‑1 (MQ=255)
                                                                                                                  ctctACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCa                                     >  2:1595117/1‑150 (MQ=255)
                                                                                                                  ctctACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCa                                     >  2:1594765/1‑150 (MQ=255)
                                                                                                                      aCCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGccc                                 >  1:1737233/1‑150 (MQ=255)
                                                                                                                       ccATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcccc                                >  2:4204555/1‑150 (MQ=255)
                                                                                                                       ccATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcccc                                >  2:1446885/1‑150 (MQ=255)
                                                                                                                       ccATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcccc                                <  1:5485476/150‑1 (MQ=255)
                                                                                                                        cATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGa                                                                  >  1:773992/1‑115 (MQ=255)
                                                                                                                        cATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcc                                  <  2:4917548/147‑1 (MQ=255)
                                                                                                                        cATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcc                                  >  1:4917548/1‑147 (MQ=255)
                                                                                                                          tCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTc                             <  1:1948124/150‑1 (MQ=255)
                                                                                                                               tCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATaa                        >  2:2750171/1‑150 (MQ=255)
                                                                                                                               tCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATaa                        >  1:2225357/1‑150 (MQ=255)
                                                                                                                               tCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATaa                        >  2:4069057/1‑150 (MQ=255)
                                                                                                                                 tttgCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATg                      >  1:247476/4‑150 (MQ=255)
                                                                                                                                 gAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATg                      <  1:3332795/150‑1 (MQ=255)
                                                                                                                                   atccAGGGCTGTGATTGACCATCCACGTCAAAAAGCCCTCACCCTGCTCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGc                    <  1:176533/148‑1 (MQ=255)
                                                                                                                                   aGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGc                    <  1:1172326/150‑1 (MQ=255)
                                                                                                                                    gCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt                   <  2:909711/150‑1 (MQ=255)
                                                                                                                                      cggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTAcc                                                                         <  1:599804/92‑1 (MQ=38)
                                                                                                                                      cggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTAcc                                                                         >  2:599804/3‑94 (MQ=38)
                                                                                                                                            tgtgCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCg           >  1:1972007/1‑150 (MQ=255)
                                                                                                                                            tgtgCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCg           >  1:3136219/1‑150 (MQ=255)
                                                                                                                                             gtgCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCgg          <  2:4171399/150‑1 (MQ=255)
                                                                                                                                             gtgCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCgg          <  1:4331137/150‑1 (MQ=255)
                                                                                                                                              tgCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                        >  1:1639997/1‑103 (MQ=255)
                                                                                                                                                cTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGc       <  2:3725133/150‑1 (MQ=255)
                                                                                                                                                  tGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGt     <  1:4835106/150‑1 (MQ=255)
                                                                                                                                                     ccAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTCaa  <  2:3801580/150‑1 (MQ=255)
                                                                                                                                                     ccAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTCaa  >  2:4682352/1‑150 (MQ=255)
                                                                                                                                                     ccAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTCaa  >  2:3570963/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTCAA  >  NC_002947/4654835‑4655133

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: