Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,509,440 T→C 5.7% S108G (AGC→GGC)  PP_1324 ← endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,509,4400TC 5.7% 259.6 / 4.3 88S108G (AGC→GGC) PP_1324endonuclease
Reads supporting (aligned to +/- strand):  ref base T (41/42);  new base C (3/2);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

TGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGT  >  NC_002947/1509294‑1509589
                                                                                                                                                  |                                                                                                                                                     
tGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTAcc                                                                                                                                                    <  2:2076843/150‑1 (MQ=255)
    cATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGc                                                                                                                                                <  1:1038465/150‑1 (MQ=255)
    cATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGc                                                                                                                                                >  2:1038465/1‑150 (MQ=255)
    cATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGGTCTTCAGCCATTGAACCGGTCTGCCTGCATGCTGGCTAcccccc                                                                                                                                                >  1:174362/1‑147 (MQ=255)
      tCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTAcc                                                                                                                                                    <  1:566961/144‑1 (MQ=255)
      tCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTAcc                                                                                                                                                    >  2:566961/1‑144 (MQ=255)
       cTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAgg                                                                                                                                             <  1:2236109/150‑1 (MQ=255)
          gTGACCTTAAATGGGGCGGTGGCCTTTGCTACGACCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAAAGCAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCCTTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCg                                                                                                                                          >  1:1939534/1‑150 (MQ=255)
           tGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGa                                                                                                                                         <  2:1903177/150‑1 (MQ=255)
            gACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAc                                                                                                                                        >  2:1818610/1‑150 (MQ=255)
             aCCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAcc                                                                                                                                       >  2:2048629/1‑150 (MQ=255)
             aCCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGGTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAcc                                                                                                                                       >  1:480479/1‑150 (MQ=255)
                ttAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACg                                                                                                                                    >  1:1277658/1‑150 (MQ=255)
                  aaaTGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTc                                                                                                                                  <  2:2096848/150‑1 (MQ=255)
                      ggggCGGTGGCCTTTGCTACTCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCCAAg                                                                                                                              >  2:514603/1‑150 (MQ=255)
                       gggCGGTGGCCTTTGCTAGGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGc                                                                                                                             <  2:1590044/150‑1 (MQ=255)
                              ggCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAggg                                                                                                                      <  2:480479/150‑1 (MQ=255)
                               gCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAGGCGGCAgggg                                                                                                                     >  2:806085/1‑150 (MQ=255)
                               gCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAgggg                                                                                                                     <  2:1277658/150‑1 (MQ=255)
                                         cgcgCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTTAAGCGGTCTGCCTGCATGGTGGCTACGCTGCAGGGCGACCACGTCGGAGCGGCAGAGGTAGTTGcccc                                                                                                           >  2:529150/1‑150 (MQ=255)
                                          gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGtt                                                                                                                >  2:1537609/1‑144 (MQ=255)
                                          gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGtt                                                                                                                <  1:1537609/144‑1 (MQ=255)
                                             caacaaCGACCAGCCCGCAATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTg                                                                                                       <  1:264413/150‑1 (MQ=255)
                                              aacaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTgg                                                                                                      <  1:2048629/150‑1 (MQ=255)
                                              aacaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTgg                                                                                                      <  1:434255/150‑1 (MQ=255)
                                              aacaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTgg                                                                                                      <  2:1993838/150‑1 (MQ=255)
                                               acaacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTggg                                                                                                     >  2:1319995/1‑150 (MQ=255)
                                                caacGACCAGCCCTCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                    <  1:529150/150‑1 (MQ=255)
                                                 aacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGcc                                                                                                   >  2:235813/1‑150 (MQ=255)
                                                    gACCAGCCCGCGATACATAGAGCAAAAAAAACGAATTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGcctcc                                                                                                >  1:1436291/1‑150 (MQ=255)
                                                    gACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAcc                                                                                                                                       <  1:1351468/111‑1 (MQ=255)
                                                    gACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAcc                                                                                                                                       >  2:1351468/1‑111 (MQ=255)
                                                        aGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCt                                                                                            >  1:1612884/1‑150 (MQ=255)
                                                          cccGCGATACATAGAGC‑AAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                        <  2:1836497/119‑1 (MQ=255)
                                                          cccGCGATACATAGAGC‑AAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                        >  1:1836497/1‑119 (MQ=255)
                                                           ccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                    >  2:803490/1‑139 (MQ=255)
                                                           ccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                    <  1:803490/139‑1 (MQ=255)
                                                           ccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCa                                                                                         <  1:1512605/150‑1 (MQ=255)
                                                                  aCATAGAGCAAAAAAAACGAATTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAg                                                                                  <  2:1436291/150‑1 (MQ=255)
                                                                  aCATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAg                                                                                  >  1:1771330/1‑150 (MQ=255)
                                                                   cATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCGCAGTGGGCGTCCTTCTGAAGGaaaggg                                                                                 >  2:916328/1‑146 (MQ=255)
                                                                   cATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAgg                                                                                 <  2:1771330/150‑1 (MQ=255)
                                                                       gagCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTGTGTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTgg                                                                             <  2:187071/150‑1 (MQ=255)
                                                                          cAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAGGGTAGCCACCATGCAg                                                                                                                                              >  1:467683/1‑82 (MQ=38)
                                                                          cAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAGGGTAGCCACCATGCAg                                                                                                                                              <  2:467683/82‑1 (MQ=38)
                                                                           aaaaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTggcggc                                                                         >  2:1059041/1‑150 (MQ=255)
                                                                           aaaaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGCTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTggcggc                                                                         <  1:1319995/150‑1 (MQ=255)
                                                                             aaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCaa                                                                       <  1:128694/150‑1 (MQ=255)
                                                                                 aaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTCGCTGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCAc                                                                   <  1:614726/150‑1 (MQ=255)
                                                                                  aCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCAcc                                                                  <  2:1612884/150‑1 (MQ=255)
                                                                                    gAAGTGGGTTCAGCATTCGAAGGCGTTTTTTTTTTTTTTTTTTGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCaa                                                                >  2:976992/1‑150 (MQ=255)
                                                                                    gAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCGGCAGGAACAGGGTGGCGGCAAGCACCac                                                                >  2:1775269/1‑149 (MQ=255)
                                                                                     aaGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAAc                                                               >  1:1869100/1‑150 (MQ=255)
                                                                                       gTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCg                                                             >  2:340275/1‑150 (MQ=255)
                                                                                                    tCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTgcgc                                                >  2:1974842/1‑150 (MQ=255)
                                                                                                      gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                    >  2:1848163/1‑96 (MQ=255)
                                                                                                      gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                    <  1:1848163/96‑1 (MQ=255)
                                                                                                      gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTgg                                                                             >  2:257859/1‑119 (MQ=255)
                                                                                                      gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTgg                                                                             <  1:257859/119‑1 (MQ=255)
                                                                                                       aaGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAgg                                                                                 >  2:690624/1‑114 (MQ=255)
                                                                                                       aaGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAgg                                                                                 <  1:690624/114‑1 (MQ=255)
                                                                                                           cGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGa                                         <  1:1137456/150‑1 (MQ=255)
                                                                                                            gTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGAt                                        >  1:1298787/1‑150 (MQ=255)
                                                                                                               tttCAGCCATTGAAGCGGTCGGCCATCATGATGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACt                                     <  1:514603/150‑1 (MQ=255)
                                                                                                                  cAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTc                                  <  2:1869100/150‑1 (MQ=255)
                                                                                                                    gCCATTGAAGCTGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCga                                <  1:916328/150‑1 (MQ=255)
                                                                                                                    gCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCATGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCga                                >  1:2071449/1‑150 (MQ=255)
                                                                                                                     ccATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAGGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCgag                               <  1:806085/150‑1 (MQ=255)
                                                                                                                          gAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGAGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACAAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGCGAgcgc                          >  2:2230647/1‑150 (MQ=255)
                                                                                                                            aGCGGTCTGCCTGCATGGTGGCTACCCTGCATGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCt                        <  2:2071449/150‑1 (MQ=255)
                                                                                                                              cGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTcc                      <  1:1818610/150‑1 (MQ=255)
                                                                                                                                gTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGa                    >  1:290503/1‑150 (MQ=255)
                                                                                                                                gTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGCAGATGCGTTGGTGAGCGCCTCCGa                    >  2:540962/1‑150 (MQ=255)
                                                                                                                                 tCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGaa                   <  1:976992/150‑1 (MQ=255)
                                                                                                                                     ccttcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGccc               <  1:588313/150‑1 (MQ=255)
                                                                                                                                      ctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCcg              >  2:1424517/1‑150 (MQ=255)
                                                                                                                                      ctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTc                                                                                             <  2:1412188/71‑1 (MQ=37)
                                                                                                                                      ctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTc                                                                                             >  1:1412188/1‑71 (MQ=37)
                                                                                                                                      ctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGTGCGCCTCCGAAGcccc              >  2:1045913/1‑149 (MQ=255)
                                                                                                                                          aTGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTg          >  1:1284778/1‑150 (MQ=255)
                                                                                                                                           tggtggCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGccc                                              <  2:868753/113‑1 (MQ=255)
                                                                                                                                           tggtggCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGGGGCGGCAAGCACCAACCGTTTCTGCTTGCGccc                                              >  1:868753/1‑113 (MQ=255)
                                                                                                                                              tggCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAAc      <  2:174362/150‑1 (MQ=255)
                                                                                                                                                gCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCg    >  2:2245399/1‑150 (MQ=255)
                                                                                                                                                gCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCg    >  2:2080526/1‑150 (MQ=255)
                                                                                                                                                gCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCg    >  2:2023385/1‑150 (MQ=255)
                                                                                                                                                 cTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCgg   <  1:235813/150‑1 (MQ=255)
                                                                                                                                                 cTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCgg   <  1:1424517/150‑1 (MQ=255)
                                                                                                                                                  tACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGt  >  1:1971993/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGTCCATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGT  >  NC_002947/1509294‑1509589

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