Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,485,560 T→G 6.8% C146G (TGC→GGC)  PP_3972 → short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,485,5600TG 6.8% 1327.0 / 59.6 393C146G (TGC→GGC) PP_3972short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (179/186);  new base G (12/15);  total (191/202)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 399 total aligned reads displayed.

TCATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGCC  >  NC_002947/4485411‑4485709
                                                                                                                                                     |                                                                                                                                                     
    cGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCaa                                                                                                                                                   >  2:1315571/1‑150 (MQ=255)
     gCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCaaa                                                                                                                                                  >  1:294181/1‑150 (MQ=255)
     gCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCaaa                                                                                                                                                  <  2:5969257/150‑1 (MQ=255)
       aCCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAgg                                                                                                                                                >  1:2708523/1‑150 (MQ=255)
       aCCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCGTTGCCTGccatgg                                                                                                                                                >  1:4489708/1‑145 (MQ=255)
        cccTTCGAACCCCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgc                                                                                                                                               <  2:529362/150‑1 (MQ=255)
        cccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGTATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgc                                                                                                                                               <  1:3904789/150‑1 (MQ=255)
          cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgcgc                                                                                                                                             <  1:6515407/150‑1 (MQ=255)
          cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgcgt                                                                                                                                             <  1:741273/150‑2 (MQ=255)
          cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgcgt                                                                                                                                             <  1:6978768/150‑2 (MQ=255)
           ttCGAACACCGTAGTAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                            >  2:7999660/1‑150 (MQ=255)
           ttCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                            >  1:5409034/1‑150 (MQ=255)
           ttCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                            >  1:1771901/1‑150 (MQ=255)
           ttCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                            >  2:2643626/1‑150 (MQ=255)
           ttCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                            >  2:2014618/1‑150 (MQ=255)
               aaCACCGTAGCGGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCt                                                                                                                                        >  1:4932033/1‑150 (MQ=255)
               aaCACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCGCGCCTTTGCAGGCAAAGgcgtg                                                                                                                                            <  1:8013167/146‑3 (MQ=255)
                acacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCtt                                                                                                                                       <  1:1263495/150‑1 (MQ=255)
                 cacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                                      <  1:8608431/150‑1 (MQ=255)
                     gTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGAAATCCCGTCAGGGGGCCGCTGCCTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGc                                                                                                                                  >  2:333411/1‑150 (MQ=255)
                        gcagcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCg                                                                                                                               <  1:6460894/150‑1 (MQ=255)
                          agcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCa                                                                                                                             <  2:353623/150‑1 (MQ=255)
                            cGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGc                                                                                                                           <  2:3405037/150‑1 (MQ=255)
                             ggCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGcc                                                                                                                          >  1:7038767/1‑150 (MQ=255)
                                acaGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCcgc                                                                                                                       <  2:8071595/150‑1 (MQ=255)
                                  cgagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAACGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGa                                                                                                                     <  2:5146187/149‑1 (MQ=255)
                                  agagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGa                                                                                                                     <  1:3303241/150‑1 (MQ=255)
                                  agagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGa                                                                                                                     <  1:5123762/150‑1 (MQ=255)
                                   gagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGAt                                                                                                                    >  1:2014080/1‑150 (MQ=255)
                                   gagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGAt                                                                                                                    >  1:1107963/1‑150 (MQ=255)
                                    agCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATg                                                                                                                   >  1:4110488/1‑150 (MQ=255)
                                       aagaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaag                                                                                                                >  1:7639820/1‑150 (MQ=255)
                                         gaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaagaa                                                                                                              <  2:7787264/150‑1 (MQ=255)
                                         gaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaagaa                                                                                                              <  2:954460/150‑1 (MQ=255)
                                           agCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGc                                                                                                            <  1:3769780/150‑1 (MQ=255)
                                              tCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCatt                                                                                                         >  1:4290737/1‑150 (MQ=255)
                                              tCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCatt                                                                                                         >  1:153338/1‑150 (MQ=255)
                                               ccccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCatta                                                                                                        <  2:44103/150‑1 (MQ=255)
                                                cccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCattat                                                                                                       <  2:7471661/150‑1 (MQ=255)
                                                 ccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCattatt                                                                                                      >  1:1403819/1‑150 (MQ=255)
                                                         gATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTc                                                                                              <  2:2075027/150‑1 (MQ=255)
                                                           tatCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTg                                                                                            <  2:8700401/150‑1 (MQ=255)
                                                           tatCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTg                                                                                            <  2:6686587/150‑1 (MQ=255)
                                                           tatCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTg                                                                                            <  2:4268477/150‑1 (MQ=255)
                                                            atCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTgg                                                                                           >  2:5688297/1‑150 (MQ=255)
                                                             tCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGa                                                                                          <  1:3139109/150‑1 (MQ=255)
                                                             tCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGa                                                                                          <  2:1768874/150‑1 (MQ=255)
                                                                      ggggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCAc                                                                                 >  2:4569756/1‑150 (MQ=255)
                                                                      ggggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCAc                                                                                 >  1:6862221/1‑150 (MQ=255)
                                                                      ggggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCAc                                                                                 <  1:2870598/150‑1 (MQ=255)
                                                                           cGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGAtt                                                                            <  1:7163751/150‑1 (MQ=255)
                                                                              tGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgtg                                                                         >  2:104920/1‑150 (MQ=255)
                                                                              tGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgtg                                                                         <  1:944563/150‑1 (MQ=255)
                                                                                 gTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTa                                                                      <  1:5459385/150‑1 (MQ=255)
                                                                                 gTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTa                                                                      <  1:7910033/150‑1 (MQ=255)
                                                                                   ggCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACg                                                                    <  2:6523124/150‑1 (MQ=255)
                                                                                      aaaGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGAt                                                                 >  2:7179388/1‑150 (MQ=255)
                                                                                      aaaGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGAt                                                                 <  2:8761998/150‑1 (MQ=255)
                                                                                      aaaGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGAt                                                                 <  2:2994554/150‑1 (MQ=255)
                                                                                      aaaGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGAt                                                                 >  2:6095612/1‑150 (MQ=255)
                                                                                            gtgGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATa                                                           >  2:4057352/1‑150 (MQ=255)
                                                                                            gtgGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATa                                                           >  1:3672853/1‑150 (MQ=255)
                                                                                               gTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCg                                                        >  2:2240364/1‑150 (MQ=255)
                                                                                               gTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCg                                                        <  2:1205741/150‑1 (MQ=255)
                                                                                                 aaTCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGt                                                      >  2:1860493/1‑150 (MQ=255)
                                                                                                   tCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTca                                                    >  1:6148250/1‑150 (MQ=255)
                                                                                                     acaGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTcacc                                                  >  1:4808321/1‑150 (MQ=255)
                                                                                                       aGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTcaccac                                                <  2:3736453/150‑1 (MQ=255)
                                                                                                        ggcgCATCCAGCGGCAGAGGCCGCGCCGCCGCGCACGCCCTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCCGGATGAATGCACCGATAGTGGTACGGATGCCATAGCGGTCaccacc                                               <  2:4700602/147‑1 (MQ=255)
                                                                                                         gggCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCg                                              >  1:101929/1‑150 (MQ=255)
                                                                                                           gCATCCAGCGGCATCGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCCTTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGAc                                            <  2:7221534/150‑1 (MQ=255)
                                                                                                           gCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGAc                                            <  1:5947763/150‑1 (MQ=255)
                                                                                                             aTCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGt                                          >  1:8692297/1‑150 (MQ=255)
                                                                                                              tCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTc                                         >  1:5830123/1‑150 (MQ=255)
                                                                                                               ccAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCa                                        >  1:8523614/1‑150 (MQ=255)
                                                                                                                cAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCAc                                       <  1:6727345/150‑1 (MQ=255)
                                                                                                                 aGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCAcc                                      >  1:2507266/1‑150 (MQ=255)
                                                                                                                    ggCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCcac                                   <  2:2931163/150‑1 (MQ=255)
                                                                                                                     gCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCcaca                                  <  1:1471312/150‑1 (MQ=255)
                                                                                                                      cATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAg                                 <  2:4994127/150‑1 (MQ=255)
                                                                                                                       aTAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTCCGGATGCCATAGCGGTCACCACCGACGTCACCCACAGc                                >  2:1475345/1‑150 (MQ=255)
                                                                                                                       aTAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGc                                >  1:6553735/1‑150 (MQ=255)
                                                                                                                         aGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGa                              <  2:764757/150‑1 (MQ=255)
                                                                                                                           gccgcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGAcc                            >  2:2497848/1‑150 (MQ=255)
                                                                                                                             cgcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAg                          <  2:6148250/150‑1 (MQ=255)
                                                                                                                             cgcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAg                          <  1:2240364/150‑1 (MQ=255)
                                                                                                                                  cgccgcGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCa                     >  1:2081376/1‑150 (MQ=255)
                                                                                                                                   gccgcGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCaa                    <  2:6402893/150‑1 (MQ=255)
                                                                                                                                   gccgcGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCaa                    <  2:2808707/150‑1 (MQ=255)
                                                                                                                                     cgcgCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAgc                  >  1:7950520/1‑150 (MQ=255)
                                                                                                                                        gcACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCt               <  1:8228206/150‑1 (MQ=255)
                                                                                                                                         tgcGCCTTTGCAGGCAAAGGCGTGCGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTg              >  2:8512614/3‑150 (MQ=255)
                                                                                                                                           cGCCTTTGCAGGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCGGGTTCAAGCGCTGgc            >  2:5583551/1‑150 (MQ=255)
                                                                                                                                            gCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgcc           >  1:6254684/1‑150 (MQ=255)
                                                                                                                                             ccTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccg          <  2:6264618/150‑1 (MQ=255)
                                                                                                                                              cTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccgc         >  2:284696/1‑150 (MQ=255)
                                                                                                                                              cTTTGCAGGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccgc         >  1:3744039/1‑150 (MQ=255)
                                                                                                                                              cTTTGCAGGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccgc         >  1:6341900/1‑150 (MQ=255)
                                                                                                                                                    cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGc   >  2:823864/1‑150 (MQ=255)
                                                                                                                                                    cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGc   >  1:4008892/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGCC  >  NC_002947/4485411‑4485709

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