Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 243,508 T→A 12.2% *221C (TGA→TGT) ‡ fkl ← FKBP‑type peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947243,5080TA12.2% 243.2 / 24.2 90*221C (TGA→TGT) ‡fklFKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base T (37/42);  new base A (4/7);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01

GGGAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGAC  >  NC_002947/243360‑243654
                                                                                                                                                    |                                                                                                                                                  
gggAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAAAGAACTCAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTc                                                                                                                                                   >  1:1843076/1‑150 (MQ=255)
 ggAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCANAACGCACACNAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCa                                                                                                                                                  <  1:478346/150‑1 (MQ=255)
 ggAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCa                                                                                                                                                  <  1:1424667/150‑1 (MQ=255)
       cAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGGGCAGGCCTGAGGTTGGCGTCAGTCGcc                                                                                                                                            >  2:1275932/1‑150 (MQ=255)
           gACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGACTGACGCCAAcc                                                                                                                                        <  1:1235654/150‑1 (MQ=255)
            acacCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGACTGACGCCAAcct                                                                                                                                       <  2:1124742/150‑2 (MQ=255)
              acCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGcc                                                                                                                                     >  1:301029/1‑150 (MQ=255)
              acCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGACTGACGCCAAcctca                                                                                                                                     <  2:163684/150‑4 (MQ=255)
               ccGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGACTGACGCCAAcctcag                                                                                                                                    <  1:188881/150‑5 (MQ=255)
                      cTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTc                                                                                                                             <  2:2362352/150‑1 (MQ=255)
                      cTTGGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTc                                                                                                                             <  2:316334/150‑1 (MQ=255)
                       tttGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGcc                                                                                                                                     >  1:210202/1‑141 (MQ=255)
                       tttGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGcc                                                                                                                                     <  2:210202/141‑1 (MQ=255)
                       tttGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCg                                                                                                                            <  2:860131/150‑1 (MQ=255)
                       tttGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCg                                                                                                                            <  1:2194311/150‑1 (MQ=255)
                       tttGCCGAGAAGCTGGCCAAATAGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCg                                                                                                                            >  2:1027937/1‑150 (MQ=255)
                            cGAGAAGCTGGCCAAGTTGAGAAGGTTGTTGCATGAGCAATTCGAGCTTGAAGGCTGCCTGGTCGAAGTGCTGCACAACGCACACAATGGAGAGGGTACAGTGCAGGCCTGAGGGTGGCGTCAATAGCCAACCGCCAGCAGTTCGAtctc                                                                                                                       >  2:977489/1‑150 (MQ=255)
                            cGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGAtctc                                                                                                                       <  1:1904468/150‑1 (MQ=255)
                             gagaAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCAAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCg                                                                                                                      >  2:275043/1‑150 (MQ=255)
                               gaAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGaa                                                                                                                    >  1:1263055/1‑150 (MQ=255)
                                 aGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATa                                                                                                                  <  2:1793597/150‑1 (MQ=255)
                                    tGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGtt                                                                                                                              <  1:1175252/135‑1 (MQ=255)
                                    gggCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGtt                                                                                                                              >  2:1175252/2‑135 (MQ=255)
                                       ccAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCGACCGCCAGCAGTTCGATCTCGAATACCAAgg                                                                                                            >  2:1249303/1‑150 (MQ=255)
                                       ccAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAgg                                                                                                            <  2:1823328/150‑1 (MQ=255)
                                        cAAATTTGGCGGGTTGTTGCCTGCGCGATTAGCGATTGGGCGCTGCTTGGCCTAATCGCCGCACATCTGACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAggg                                                                                                           <  2:1843076/150‑1 (MQ=255)
                                        cAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACACGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAggg                                                                                                           <  1:161230/150‑1 (MQ=255)
                                        cAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATAc                                                                                                                 <  1:1181275/144‑1 (MQ=255)
                                        cAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATAc                                                                                                                 >  2:1181275/1‑144 (MQ=255)
                                        cAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAggg                                                                                                           >  2:1291897/1‑150 (MQ=255)
                                        cAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAggg                                                                                                           <  2:1383983/150‑1 (MQ=255)
                                         aaaTTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAgggg                                                                                                          <  2:316409/150‑1 (MQ=255)
                                           aTTGGGCGGGTAGTTGCAGGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGgt                                                                                                        <  2:60623/150‑1 (MQ=255)
                                              gggCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTa                                                                                                     <  1:1027937/150‑1 (MQ=255)
                                              gggCGGGTTGTGGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTa                                                                                                     <  2:301029/150‑1 (MQ=255)
                                                    gttgttGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCaa                                                                                                                                          >  2:2037517/1‑107 (MQ=255)
                                                    gttgttGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCaa                                                                                                                                          <  1:2037517/107‑1 (MQ=255)
                                                        ttgCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATc                                                                                           >  2:2340508/1‑150 (MQ=255)
                                                        ttgCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATc                                                                                           >  2:417790/1‑150 (MQ=255)
                                                          gCATGAGCGCGTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAg                                                                                         <  2:119610/150‑1 (MQ=255)
                                                           cATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAg                                                                                         <  1:2346104/149‑1 (MQ=255)
                                                           cATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAg                                                                                         >  2:2346104/1‑149 (MQ=255)
                                                            aTGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACAAAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGaa                                                                                                                    >  1:1058267/1‑121 (MQ=255)
                                                            aTGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACAAAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGaa                                                                                                                    <  2:1058267/121‑1 (MQ=255)
                                                                   cTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGc                                                                                >  2:2114169/1‑150 (MQ=255)
                                                                    ttCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAGGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAgggg                                                                                                          <  1:1735933/123‑1 (MQ=255)
                                                                    ttCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAGGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAgggg                                                                                                          >  2:1735933/1‑123 (MQ=255)
                                                                     tCGAGCTTGAAGATTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGccc                                                                              >  2:1268796/1‑150 (MQ=255)
                                                                              aaGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAg                                                                                                             <  2:1201693/110‑1 (MQ=255)
                                                                              aaGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAg                                                                                                             >  1:1201693/1‑110 (MQ=255)
                                                                                gACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAgggg                                                                                                          >  2:995836/1‑111 (MQ=255)
                                                                                gACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAgggg                                                                                                          <  1:995836/111‑1 (MQ=255)
                                                                                   tGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAg                                                                <  1:1829004/150‑1 (MQ=255)
                                                                                   tGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAg                                                                <  2:1010766/150‑1 (MQ=255)
                                                                                   tGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAg                                                                >  2:1829004/1‑150 (MQ=255)
                                                                                        gATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTa                                                                                                     <  2:1667275/108‑1 (MQ=255)
                                                                                        gATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTa                                                                                                     >  1:1667275/1‑108 (MQ=255)
                                                                                        gATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTg                                                           >  2:2324899/1‑150 (MQ=255)
                                                                                          tCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGACTCCCGCCCCCTCGGCGCCATAGGCTTgcg                                                         >  2:289164/1‑150 (MQ=255)
                                                                                             aaGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTg                                                      >  2:1789277/1‑150 (MQ=255)
                                                                                               gTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATGAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGaa                                                    >  2:1413034/1‑150 (MQ=255)
                                                                                                tgctgcACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAg                                                   <  1:1874578/150‑1 (MQ=255)
                                                                                                    gcACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAAt                                               <  1:2181670/150‑1 (MQ=255)
                                                                                                     cacaACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCa                                                                   <  1:1674211/129‑1 (MQ=255)
                                                                                                     cacaACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCa                                                                   >  2:1674211/1‑129 (MQ=255)
                                                                                                       caACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTAc                                            <  2:1783481/150‑1 (MQ=255)
                                                                                                        aaCGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTAcc                                           >  2:954230/1‑150 (MQ=255)
                                                                                                             acacaAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAgg                                                                                        >  1:2289586/1‑100 (MQ=255)
                                                                                                             acacaAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAgg                                                                                        <  2:2289586/100‑1 (MQ=255)
                                                                                                                caAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCa                                   >  1:2014654/1‑150 (MQ=255)
                                                                                                                   ggAAGAGGGTGCAGTGCAGGCATGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGTAATCAGGTCTCCCGCCCCCTCGCCGCCATAGGCTTGCGCTGAAGGGATTAGCAGGCGCCAttt                                <  2:937840/150‑1 (MQ=255)
                                                                                                                    gAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAg                                         >  1:282515/1‑140 (MQ=255)
                                                                                                                    gAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAg                                         <  2:282515/140‑1 (MQ=255)
                                                                                                                    gAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGAATTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTg                               <  1:1341393/150‑1 (MQ=255)
                                                                                                                           gtgcagtgcaGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCgg                        >  2:136541/1‑150 (MQ=255)
                                                                                                                           gtgcagtgcaGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCAAATGGCGcctg                        <  2:1563785/150‑3 (MQ=255)
                                                                                                                           gtgcagtgcaGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCAAATGGCGcctg                        <  1:275043/150‑3 (MQ=255)
                                                                                                                           gtgcagtgcaGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGAGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTAAAGGAATTACCAGGCGCCAAATGGCGcctg                        <  2:1775362/150‑3 (MQ=255)
                                                                                                                             gcagtgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGgtg                      >  1:968379/1‑150 (MQ=255)
                                                                                                                                     tgctggcggTTGGCGACTGACGCCAACCTCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCAtttt              >  1:2190185/8‑150 (MQ=255)
                                                                                                                                      gctggcggTTGGCGACTGACGCCAACCTCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTc             <  2:2190185/144‑1 (MQ=255)
                                                                                                                                       ccTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATATAGTATATAATGGGGGTGTAGGGTGCAAACAGGGCTCCCGCCACCTAGGCGCCATGGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCg            >  1:1863391/1‑150 (MQ=255)
                                                                                                                                        tggcggTTGGCGACTGACGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTgc                                                          <  2:2083934/99‑1 (MQ=255)
                                                                                                                                        tggcggTTGGCGACTGACGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTgc                                                          >  1:2083934/5‑103 (MQ=255)
                                                                                                                                        cTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCgg           <  2:2429435/150‑1 (MQ=255)
                                                                                                                                          gcggTTGGCGACTGACGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGaa                                                    >  1:2033593/3‑107 (MQ=255)
                                                                                                                                          gcggTTGGCGACTGACGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGaa                                                    <  2:2033593/105‑1 (MQ=255)
                                                                                                                                           cggTTGGCGACTGACGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCaa        >  1:1335113/2‑150 (MQ=255)
                                                                                                                                               tGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCg    >  2:2074031/1‑150 (MQ=255)
                                                                                                                                                 gCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGAc  <  1:1413034/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGGAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGAC  >  NC_002947/243360‑243654

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