Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,509,440 T→C 11.6% 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,4400TC11.6% 227.2 / 16.8 86S108G (AGC→GGC) PP_1324endonuclease
Reads supporting (aligned to +/- strand):  ref base T (39/37);  new base C (5/5);  total (44/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 = 7.80e-01

TCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGT  >  NC_002947/1509300‑1509589
                                                                                                                                            |                                                                                                                                                     
tCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAACCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCAGGGTAGCCACCATGCAg                                                                                                                                              <  1:1945166/150‑1 (MQ=255)
    gTGACCTTAAATGGGGCGGGGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCg                                                                                                                                          >  1:1772069/1‑150 (MQ=255)
          ttAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACg                                                                                                                                    >  2:2215245/1‑150 (MQ=255)
            aaaTGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCg                                                                                                                                 >  1:728888/1‑150 (MQ=255)
                 gggCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAg                                                                                                                              <  1:2232566/149‑1 (MQ=255)
                 gggCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAg                                                                                                                              >  2:2232566/1‑149 (MQ=255)
                       tGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAgg                                                                                                                       <  2:1186124/150‑1 (MQ=255)
                       tGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAgg                                                                                                                       >  2:1682302/1‑150 (MQ=255)
                          ccTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGt                                                                                                                    >  2:1276195/1‑150 (MQ=255)
                           cTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTa                                                                                                                   <  2:1495856/150‑1 (MQ=255)
                            tttGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAg                                                                                                                  >  1:91822/1‑150 (MQ=255)
                            tttGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGCGAATCCTTATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAg                                                                                                                  >  2:64906/1‑150 (MQ=255)
                              tGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGtt                                                                                                                <  1:2384470/150‑1 (MQ=255)
                              tGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGtt                                                                                                                <  1:2215245/150‑1 (MQ=255)
                                cTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGc                                                                                                              <  2:1415269/150‑1 (MQ=255)
                                cTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGc                                                                                                              <  1:1094117/150‑1 (MQ=255)
                                    gcgcAACAACGACCAGCCCGCGATACATGGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCa                                                                                                          >  1:1646635/1‑150 (MQ=255)
                                    gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTc                                                                                                                                  >  1:1342724/1‑126 (MQ=255)
                                    gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTc                                                                                                                                  <  2:1342724/126‑1 (MQ=255)
                                    gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                        >  2:1183586/1‑136 (MQ=255)
                                    gcgcAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                        <  1:1183586/136‑1 (MQ=255)
                                           aacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAc                                                                                                                                        >  2:1081700/1‑113 (MQ=255)
                                           aacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGAc                                                                                                                                        <  1:1081700/113‑1 (MQ=255)
                                           aacGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGcc                                                                                                   <  2:1772069/150‑1 (MQ=255)
                                                ccAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCa                                                                                                          >  1:1281920/1‑138 (MQ=255)
                                                ccAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCa                                                                                                          <  2:1281920/138‑1 (MQ=255)
                                                ccAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCtt                                                                                              >  2:652434/1‑150 (MQ=255)
                                                    cccGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCAAGGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGc                                                                                          >  1:739844/1‑150 (MQ=255)
                                                      cgcgATACATAGAACAAAAAAAACGAAGTGGGTTCAGTATTCGAAGGCGTTTTTTAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCCTCTGCAg                                                                                        <  2:1645768/150‑1 (MQ=255)
                                                        cgATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGa                                                                                      <  2:91822/150‑1 (MQ=255)
                                                         gATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGTTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGaa                                                                                     <  1:167589/150‑1 (MQ=255)
                                                           tACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTCGCCCCAGTGGGCCTCCTTCTGCAGGAACa                                                                                   <  1:64906/150‑1 (MQ=255)
                                                              aTAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAg                                                                                                                                              <  1:811013/88‑1 (MQ=255)
                                                              aTAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAg                                                                                                                                              >  2:811013/1‑88 (MQ=255)
                                                              aTAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAggg                                                                                >  1:1599775/1‑150 (MQ=255)
                                                               tAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGTCTACCCTGCAGTGGGACCACATCGAATCGGCAGGGGTAGTTACCCCAGTGGGCCTCCTGCTGCAGGAACAgggg                                                                               >  2:1707111/1‑149 (MQ=255)
                                                                agagCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAg                                                                                                                  <  2:128262/114‑1 (MQ=255)
                                                                agagCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAg                                                                                                                  >  1:128262/1‑114 (MQ=255)
                                                                 gagCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                        >  2:258868/1‑107 (MQ=255)
                                                                 gagCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAg                                                                                                                        <  1:258868/107‑1 (MQ=255)
                                                                 gagCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGATCCCTAGCCTGCAAAGAGACCCCGCCGCACCTGAAGGTCGAGTTGGCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTgg                                                                             >  1:1143012/1‑150 (MQ=255)
                                                                 gagCAAAAAAAACGAAGGGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTgg                                                                             >  2:1235377/1‑150 (MQ=255)
                                                                   gCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTggcg                                                                           <  1:652434/150‑1 (MQ=255)
                                                                   caaaaaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGt                                                                                                        >  1:997852/3‑121 (MQ=255)
                                                                   caaaaaaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGt                                                                                                        <  2:997852/119‑1 (MQ=255)
                                                                        aaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCa                                                                                                                                               <  1:1209567/77‑1 (MQ=255)
                                                                        aaaaaCGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCa                                                                                                                                               >  2:1209567/1‑77 (MQ=255)
                                                                              gAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCaa                                                                <  1:1235377/150‑1 (MQ=255)
                                                                              gAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCGCCCTGCAGGGTAGCCACCATGCa                                                                                                                                               <  2:2116440/71‑1 (MQ=18)
                                                                              gAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCGCCCTGCAGGGTAGCCACCATGCa                                                                                                                                               >  1:2116440/1‑71 (MQ=18)
                                                                               aaGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAAc                                                               <  1:1682302/150‑1 (MQ=255)
                                                                                 gTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGc                                                                                          >  2:1773385/1‑121 (MQ=255)
                                                                                 gTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGc                                                                                          <  1:1773385/121‑1 (MQ=255)
                                                                                      ttCAGCATTCGAAGGCGGTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCt                                                        >  1:2268335/1‑150 (MQ=255)
                                                                                        cAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGc                                                      <  2:368568/150‑1 (MQ=255)
                                                                                          gCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCtt                                                    <  2:728888/150‑1 (MQ=255)
                                                                                                            cAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTc                                  <  2:1519396/150‑1 (MQ=255)
                                                                                                             aGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCg                                 >  1:2060102/1‑150 (MQ=255)
                                                                                                              gCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTg                                                                              >  1:81289/1‑104 (MQ=255)
                                                                                                              gCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTg                                                                              <  2:81289/104‑1 (MQ=255)
                                                                                                              gCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCga                                >  2:2378226/1‑150 (MQ=255)
                                                                                                               ccATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCgag                               >  2:1964424/1‑150 (MQ=255)
                                                                                                               ccATTGAAGCGGTCTGCCTGCATGGTGGCTACACTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCAACAACCGTTTCTGCTTGCGCCCGTAGATAATGTCgag                               >  1:79505/1‑150 (MQ=255)
                                                                                                                cATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCgaga                              <  2:1646635/150‑1 (MQ=255)
                                                                                                                  ttGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCa                                                                                         <  2:1430414/89‑1 (MQ=255)
                                                                                                                  ttGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCa                                                                                         >  1:1430414/1‑89 (MQ=255)
                                                                                                                        cGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTcc                      >  1:2198289/1‑150 (MQ=255)
                                                                                                                          gTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGa                    >  2:9556/1‑150 (MQ=255)
                                                                                                                             tgcctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAgg                                                                                 <  1:1811727/86‑1 (MQ=255)
                                                                                                                             tgcctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAgg                                                                                 >  2:1811727/1‑86 (MQ=255)
                                                                                                                              gcctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGcc                >  1:735254/1‑150 (MQ=255)
                                                                                                                               cctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGccc               >  1:2191391/1‑150 (MQ=255)
                                                                                                                               cctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCaa                                                                       <  2:1654293/94‑1 (MQ=38)
                                                                                                                               cctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCaa                                                                       >  1:1654293/1‑94 (MQ=38)
                                                                                                                               cctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGttt                                                          <  2:478638/107‑1 (MQ=255)
                                                                                                                               cctgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGttt                                                          >  1:478638/1‑107 (MQ=255)
                                                                                                                               cccgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGccc               >  2:2265339/4‑150 (MQ=255)
                                                                                                                                ctgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCcg              >  1:1310724/1‑150 (MQ=255)
                                                                                                                                 tgcATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCcgc             <  2:2040786/150‑1 (MQ=255)
                                                                                                                                 tgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGttt                                                          >  1:1247333/1‑105 (MQ=255)
                                                                                                                                 tgcAGGGTAGCCACCATGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGa                                         >  1:799838/1‑122 (MQ=255)
                                                                                                                                 tgcAGGGTAGCCACCATGCAGGGCGACCACGGCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGa                                         <  2:799838/122‑1 (MQ=255)
                                                                                                                                   cATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGt           <  1:1964424/150‑1 (MQ=255)
                                                                                                                                        tggCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAAc      <  1:9556/150‑1 (MQ=255)
                                                                                                                                          gCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGCCAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCgg   <  1:1707111/150‑1 (MQ=255)
                                                                                                                                           cTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCgg   <  2:735254/150‑1 (MQ=255)
                                                                                                                                            tACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCa                                                                 <  2:1167093/87‑1 (MQ=255)
                                                                                                                                            tACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCa                                                                 >  1:1167093/1‑87 (MQ=255)
                                                                                                                                            tACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAgc                            >  1:434930/1‑124 (MQ=255)
                                                                                                                                            tACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAgc                            <  2:434930/124‑1 (MQ=255)
                                                                                                                                            tACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGt  >  2:1515034/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
TCTGGTGACCTTAAATGGGGCGGTGGCCTTTGCTACGCGCAACAACGACCAGCCCGCGATACATAGAGCAAAAAAAACGAAGTGGGTTCAGCATTCGAAGGCGTTTTTCAGCCATTGAAGCGGTCTGCCTGCATGGTGGCTACCCTGCAGGGCGACCACGTCGAAGCGGCAGGGGTAGTTGCCCCAGTGGGCCTCCTTCTGCAGGAACAGGGTGGCGGCAAGCACCAACCGTTTCTGCTTGCGCCCGTCGATACTGTCGAGAGCGCCTCCGAAGCCCGCGTGCAACCGGT  >  NC_002947/1509300‑1509589

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