Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,509,444 C→A 5.9% Q106H (CAG→CAT)  PP_1324 ← endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,509,4440CA 5.9% 268.4 / 4.0 85Q106H (CAG→CAT) PP_1324endonuclease
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base A (3/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

CATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGTAGC  >  NC_002947/1509298‑1509592
                                                                                                                                                  |                                                                                                                                                    
cATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGc                                                                                                                                                   <  1:1038465/150‑1 (MQ=255)
cATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGc                                                                                                                                                   >  2:1038465/1‑150 (MQ=255)
cATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGGTCTTCAGCCATTGAACCGGTCTGCCTGCATGCTGGCTAcccccc                                                                                                                                                   >  1:174362/1‑147 (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                                                                                                                                          >  2:1351468/1‑111 (MQ=255)
                                                gACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAcc                                                                                                                                          <  1:1351468/111‑1 (MQ=255)
                                                    aGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCt                                                                                               >  1:1612884/1‑150 (MQ=255)
                                                      cccGCGATACATAGAGC‑AAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                           >  1:1836497/1‑119 (MQ=255)
                                                      cccGCGATACATAGAGC‑AAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                           <  2:1836497/119‑1 (MQ=255)
                                                       ccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                       <  1:803490/139‑1 (MQ=255)
                                                       ccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGc                                                                                                       >  2:803490/1‑139 (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                                                                                <  1:257859/119‑1 (MQ=255)
                                                                                                  gAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTgg                                                                                >  2:257859/1‑119 (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:1424517/150‑1 (MQ=255)
                                                                                                                                             cTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCgg      <  1:235813/150‑1 (MQ=255)
                                                                                                                                              tACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGt     >  1:1971993/1‑150 (MQ=255)
                                                                                                                                                cccTGCGGGGCGACCACGTCGAAGCGACCGGGGTAGTTGCCCCAGCGGGCCTCCTTCTGCAGGAACACGACGGCGGCAAGCACCAACCGCTTCCGCTTGCGCCCGCCGATACTGCCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGGAg   >  2:685985/1‑150 (MQ=255)
                                                                                                                                                 ccTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTc                                     >  1:159335/1‑115 (MQ=255)
                                                                                                                                                 ccTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTc                                     <  2:159335/115‑1 (MQ=255)
                                                                                                                                                 ccTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGTAGc  >  1:2267591/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CATCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGTAGC  >  NC_002947/1509298‑1509592

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