Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,589,666 G→C 7.2% S39S (TCC→TCG) ‡ PP_4063 ← long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,589,6660GC 7.2% 889.2 / 40.7 269S39S (TCC→TCG) ‡PP_4063long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (123/125);  new base C (7/13);  total (130/139)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

GCCAGCTGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACG  >  NC_002947/4589517‑4589814
                                                                                                                                                     |                                                                                                                                                    
gCCAGCTGCCAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCgg                                                                                                                                                      <  2:1363748/150‑1 (MQ=255)
gCCAGCTGCACGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCgg                                                                                                                                                      <  1:4359700/150‑1 (MQ=255)
gCCAGCTGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCgg                                                                                                                                                      <  2:8674749/150‑1 (MQ=255)
gCCAGCTGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCgg                                                                                                                                                      <  2:8556170/150‑1 (MQ=255)
gCCAGCTGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCgg                                                                                                                                                      <  2:7072092/150‑1 (MQ=255)
gCCAGCTGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCgg                                                                                                                                                      <  2:3952774/150‑1 (MQ=255)
gCCAGCTGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCgg                                                                                                                                                      <  1:5875430/150‑1 (MQ=255)
   aGCTGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGac                                                                                                                                                   >  2:8236860/1‑150 (MQ=255)
     cTGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGacac                                                                                                                                                 <  2:294572/150‑1 (MQ=255)
      tGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACAcc                                                                                                                                                >  2:691799/1‑150 (MQ=255)
      tGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCGGTAACGCAGGCCCTGGTGGCGGGACAcc                                                                                                                                                >  1:50215/1‑150 (MQ=255)
       gCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCa                                                                                                                                               <  1:1326482/150‑1 (MQ=255)
       gCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCa                                                                                                                                               >  1:3788481/1‑150 (MQ=255)
       gCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCa                                                                                                                                               >  1:3711830/1‑150 (MQ=255)
        cAAGATGCACCACTGGGCGCAGTTGGACGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTAGATTTCGACCTGTTCGGCCAACTGCCGCCAGCGGTAACGCAGGCCCTGGTGGCGGGACACCAg                                                                                                                                              >  2:9116960/1‑150 (MQ=255)
         aaGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGTCCCGCCACCAgg                                                                                                                                             <  2:3260842/150‑1 (MQ=255)
           gATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGTCCCGCCACCAGGGc                                                                                                                                           <  2:6611034/150‑1 (MQ=255)
            aTGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCAGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGcc                                                                                                                                          <  2:246404/150‑1 (MQ=255)
              gCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTc                                                                                                                                        >  1:3957120/1‑150 (MQ=255)
                accacTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCTGTGCGATCAGGGCTCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCATCTGTAACTCAGGCCCTGGTGGCGGGACACCAGGGCCTCGc                                                                                                                                      <  1:4241494/150‑1 (MQ=255)
                accacTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGc                                                                                                                                      >  2:3119642/1‑150 (MQ=255)
                    cTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTc                                                                                                                                  <  1:1741894/150‑1 (MQ=255)
                    cTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTc                                                                                                                                  <  2:8013148/150‑1 (MQ=255)
                            aGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGAGGCGGGACACCAGGGCCTCGCTGTCGCAACAGc                                                                                                                          >  2:772319/1‑150 (MQ=255)
                             gTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCg                                                                                                                         <  2:9013257/150‑1 (MQ=255)
                              ttGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCgg                                                                                                                        <  2:4586870/150‑1 (MQ=255)
                              ttGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCgg                                                                                                                        <  2:1753598/150‑1 (MQ=255)
                               tGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCggg                                                                                                                       <  2:2081293/150‑1 (MQ=255)
                                     aCCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACgg                                                                                                                 >  1:5754581/1‑150 (MQ=255)
                                       cAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTg                                                                                                               >  1:4674645/1‑150 (MQ=255)
                                             ccAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCTTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCg                                                                                                         >  1:3274095/1‑150 (MQ=255)
                                             ccAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCg                                                                                                         <  2:2207830/150‑1 (MQ=255)
                                               aaCACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGaa                                                                                                       <  1:8072302/150‑1 (MQ=255)
                                                    ggTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCt                                                                                                  >  1:1086832/1‑150 (MQ=255)
                                                        gCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGcc                                                                                              >  1:5709491/1‑150 (MQ=255)
                                                        gCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGcc                                                                                              >  1:6590593/1‑150 (MQ=255)
                                                            gtgtTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATg                                                                                          <  1:5175486/150‑1 (MQ=255)
                                                               ttCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTc                                                                                       >  2:5585473/1‑150 (MQ=255)
                                                                tCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCt                                                                                      >  1:8192477/1‑150 (MQ=255)
                                                                 cACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTg                                                                                     <  2:3739598/150‑1 (MQ=255)
                                                                  aCGCCCAGTGTGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTgg                                                                                    <  1:2894817/150‑1 (MQ=255)
                                                                  aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGc                                                                                                                                           <  2:3786681/95‑1 (MQ=255)
                                                                  aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTgg                                                                                    <  1:8674363/150‑1 (MQ=255)
                                                                  aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTgg                                                                                    >  1:3517317/1‑150 (MQ=255)
                                                                   cGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTggg                                                                                   <  1:1488107/150‑1 (MQ=255)
                                                                     cccAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTc                                                                                 <  2:5752826/150‑1 (MQ=255)
                                                                       cAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAg                                                                               >  1:6483710/1‑150 (MQ=255)
                                                                          tGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCaa                                                                            <  1:5416131/150‑1 (MQ=255)
                                                                             gATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGAGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAggg                                                                         >  2:4346221/1‑150 (MQ=255)
                                                                                cAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGggtcg                                                                      <  2:4393733/150‑1 (MQ=255)
                                                                                   ggCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGCAGGCCTGGCCGCTGGGCTGGGTCAGCAAGgggcgggc                                                                   >  1:8879120/1‑148 (MQ=255)
                                                                                   ggCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGggtcggtc                                                                   <  2:340428/150‑1 (MQ=255)
                                                                                      gcgGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTg                                                                <  2:3573649/150‑1 (MQ=255)
                                                                                      gcgGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTg                                                                <  1:628467/150‑1 (MQ=255)
                                                                                      gcgGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTg                                                                <  2:1091213/150‑1 (MQ=255)
                                                                                          gCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGa                                                            <  2:3861092/150‑1 (MQ=255)
                                                                                           cGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGAc                                                           >  1:5552541/1‑150 (MQ=255)
                                                                                              atatTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACcgc                                                        >  1:3438819/1‑150 (MQ=255)
                                                                                              atatTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACcgc                                                        >  1:9086056/1‑150 (MQ=255)
                                                                                               tatTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACcgcg                                                       <  1:8598249/150‑1 (MQ=255)
                                                                                                atTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGCCACCAGGGCCTCGCTGTCGCAACCGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCgg                                                      <  1:7389694/150‑1 (MQ=255)
                                                                                                 tttCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCggg                                                     >  1:4460111/1‑150 (MQ=255)
                                                                                                 tttCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCggg                                                     >  1:8044357/1‑150 (MQ=255)
                                                                                                 aaaagACCTGTTCGGCCAACTGCCGCCAGCAGTAACGCAGGCCCTGCTGTCCCGCCACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCggg                                                     >  2:5172428/5‑150 (MQ=255)
                                                                                                       ccTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGc                                               >  1:7774483/1‑150 (MQ=255)
                                                                                                         tGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGACACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTc                                             >  2:5473392/1‑150 (MQ=255)
                                                                                                          gTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCg                                            <  1:2562141/150‑1 (MQ=255)
                                                                                                           ttCGGACAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCgg                                           >  1:4146834/1‑150 (MQ=255)
                                                                                                             cGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGtt                                         <  2:8108463/150‑1 (MQ=255)
                                                                                                              ggCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCACCAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTg                                        >  2:6181794/1‑150 (MQ=255)
                                                                                                              ggCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTg                                        >  2:2548267/1‑150 (MQ=255)
                                                                                                              ggCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCAACGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTg                                        >  1:1391704/1‑150 (MQ=255)
                                                                                                               gCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGa                                       >  1:2238449/1‑150 (MQ=255)
                                                                                                               gCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGa                                       <  1:6931680/150‑1 (MQ=255)
                                                                                                                ccAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGAc                                      <  1:2627391/150‑1 (MQ=255)
                                                                                                                  aaCTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTc                                    >  2:7872109/1‑150 (MQ=255)
                                                                                                                   aCTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCa                                   <  1:7251842/150‑1 (MQ=255)
                                                                                                                     tGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATg                                 >  2:3364223/1‑150 (MQ=255)
                                                                                                                      gccgccAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATgg                                <  2:3031133/150‑1 (MQ=255)
                                                                                                                       ccgccAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATggg                               <  1:4645263/150‑1 (MQ=255)
                                                                                                                           cAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTc                           <  1:7308005/150‑1 (MQ=255)
                                                                                                                            aGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTcc                          <  1:7526056/150‑1 (MQ=255)
                                                                                                                               tGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCctt                       >  2:1510895/1‑150 (MQ=255)
                                                                                                                               tGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCctt                       >  2:59544/1‑150 (MQ=255)
                                                                                                                               tGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCctt                       >  2:4947664/1‑150 (MQ=255)
                                                                                                                               tGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCctt                       >  1:2440146/1‑150 (MQ=255)
                                                                                                                                gTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCctgc                      >  1:1911046/1‑148 (MQ=255)
                                                                                                                                 tAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCcttct                     >  2:1282709/1‑150 (MQ=255)
                                                                                                                                     gCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTtgtg                 <  1:6112828/150‑1 (MQ=255)
                                                                                                                                        ggCCCTGGTGTCCCGCCACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGgtt              >  2:5299478/1‑150 (MQ=255)
                                                                                                                                          cccTGGTGTCCCGCCACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATgg                                <  1:5270228/130‑1 (MQ=255)
                                                                                                                                          cccTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGgttgt            <  1:604379/150‑1 (MQ=255)
                                                                                                                                            cTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGttt          <  2:2849855/150‑1 (MQ=255)
                                                                                                                                             tggtggCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCcttc                      <  2:737909/137‑1 (MQ=255)
                                                                                                                                             tggtggCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGtttt         >  2:6392010/1‑150 (MQ=255)
                                                                                                                                              ggtggCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTg        <  1:2097991/150‑1 (MQ=255)
                                                                                                                                               gtggCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTgg       <  2:1964354/150‑1 (MQ=255)
                                                                                                                                               gtggCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTgg       <  1:4372962/150‑1 (MQ=255)
                                                                                                                                                   cGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCAc   <  2:6067544/150‑1 (MQ=255)
                                                                                                                                                    gggACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACg  >  1:7896635/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCAGCTGCAAGATGCACCACTGGGCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACG  >  NC_002947/4589517‑4589814

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: