Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,480,611 T→G 5.3% 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,480,6110TG 5.3% 211.1 / 4.7 75C146G (TGC→GGC) PP_3972short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (39/32);  new base G (2/2);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

ATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGG  >  NC_002947/4480464‑4480749
                                                                                                                                                   |                                                                                                                                          
atcGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGc                                                                                                                                          >  2:593333/1‑150 (MQ=255)
   gCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCaaa                                                                                                                                       >  2:42447/1‑150 (MQ=255)
      cccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgc                                                                                                                                    <  1:723086/150‑1 (MQ=255)
        cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgcgc                                                                                                                                  <  2:245800/150‑1 (MQ=255)
        cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgcgt                                                                                                                                  <  2:415407/150‑2 (MQ=255)
         ttCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                 >  1:1032386/1‑150 (MQ=255)
              acacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCGGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATAACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCtt                                                                                                                            >  1:518222/1‑150 (MQ=255)
               cacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                           <  2:1201371/150‑1 (MQ=255)
               cacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                           <  1:970335/150‑1 (MQ=255)
               cacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                           <  2:1201536/150‑1 (MQ=255)
                        agcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCa                                                                                                                  >  2:170897/1‑150 (MQ=255)
                        agcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCa                                                                                                                  >  2:417021/1‑150 (MQ=255)
                                agagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGa                                                                                                          <  2:295780/150‑1 (MQ=255)
                                 gagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGAt                                                                                                         >  1:1051493/1‑150 (MQ=255)
                                    cAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGc                                                                                                                <  1:480729/140‑1 (MQ=255)
                                    cAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGc                                                                                                                >  2:480729/1‑140 (MQ=255)
                                     aagaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaag                                                                                                     >  2:415002/1‑150 (MQ=255)
                                        aagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGaagaag                                                                                                  <  1:593333/150‑1 (MQ=255)
                                              cccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGTATCGGACGGGAATCAGCTCTGGCCTTTTGCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCattat                                                                                            >  1:414788/1‑150 (MQ=255)
                                              cccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaagaa                                                                                                   >  1:359445/1‑143 (MQ=255)
                                              cccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaagaa                                                                                                   <  2:359445/143‑1 (MQ=255)
                                              cccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCattat                                                                                            >  2:1164128/1‑150 (MQ=255)
                                                   aGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGAc                                                                                       <  1:1164128/150‑1 (MQ=255)
                                                     gggATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCat                                                                                               <  1:54448/140‑1 (MQ=255)
                                                     gggATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCat                                                                                               >  2:54448/1‑140 (MQ=255)
                                                        atatCCCGTCATGGTGCCTCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCt                                                                                  <  2:518222/150‑1 (MQ=255)
                                                          atCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTgg                                                                                >  2:869986/1‑150 (MQ=255)
                                                              cGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGa                                                                            >  1:1188274/1‑150 (MQ=255)
                                                                    ggggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCg                                                                                                                    <  2:355466/104‑1 (MQ=255)
                                                                    ggggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCg                                                                                                                    >  1:355466/1‑104 (MQ=255)
                                                                        ccGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAAt                                                                          >  1:598616/1‑142 (MQ=255)
                                                                        ccGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAAt                                                                          <  2:598616/142‑1 (MQ=255)
                                                                          gCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtg                                                                <  1:417021/150‑1 (MQ=255)
                                                                           cTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgt                                                               <  2:1188274/150‑1 (MQ=255)
                                                                                tGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTAc                                                          >  1:1037351/1‑150 (MQ=255)
                                                                                    aaaGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGAt                                                      >  2:494828/1‑150 (MQ=255)
                                                                                     aaGTTGTGCTACTCACATGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTGGCTTGGTAAGGGTCACGCATGGAGCTCGCAGACCGCGATGAAGACGAATCATTGGACGTGCTGGATGAATGCACCGATTGTGGTACGGATg                                                     <  1:415002/150‑1 (MQ=255)
                                                                                       gTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCt                                                                                  >  1:359820/1‑119 (MQ=255)
                                                                                       gTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCt                                                                                  <  2:359820/119‑1 (MQ=255)
                                                                                          gtgGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATa                                                <  2:579414/150‑1 (MQ=255)
                                                                                          gtgGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCGTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATa                                                >  2:531691/1‑150 (MQ=255)
                                                                                           tgGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgtg                                                              <  1:53764/135‑1 (MQ=255)
                                                                                           tgGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAg                                               >  1:934653/1‑150 (MQ=255)
                                                                                           tgGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAg                                               >  1:178053/1‑150 (MQ=255)
                                                                                           tgGTAATCAAAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgtg                                                              >  2:53764/1‑135 (MQ=255)
                                                                                                   acaGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGc                                                                                                 <  2:694808/92‑1 (MQ=255)
                                                                                                   acaGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGc                                                                                                 >  1:694808/1‑92 (MQ=255)
                                                                                                     aGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTcaccac                                     >  2:1143794/1‑150 (MQ=255)
                                                                                                     aGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTcaccac                                     >  2:1227/1‑150 (MQ=255)
                                                                                                       gggCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCTCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCg                                   >  1:86103/1‑150 (MQ=255)
                                                                                                         gCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGAc                                 <  1:869986/150‑1 (MQ=255)
                                                                                                            tCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCa                                                  <  1:610989/130‑1 (MQ=255)
                                                                                                            tCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCa                                                  >  2:610989/1‑130 (MQ=255)
                                                                                                             ccAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATGGTGGTACGGATGCCATAGCGGTCACCCCCGACGTcc                             <  2:700978/150‑2 (MQ=255)
                                                                                                              cAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGAt                                                                                                         <  2:411498/73‑1 (MQ=255)
                                                                                                              cAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGAt                                                                                                         >  1:411498/1‑73 (MQ=255)
                                                                                                              cAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCAc                            >  2:423694/1‑150 (MQ=255)
                                                                                                               aGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCAcc                           >  2:1166799/1‑150 (MQ=255)
                                                                                                                 cGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCca                         <  2:86103/150‑1 (MQ=255)
                                                                                                                  ggCATGGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTTCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCAc                                                                      <  1:996517/104‑1 (MQ=255)
                                                                                                                  ggCATGGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCAc                                                                      >  2:996517/1‑104 (MQ=255)
                                                                                                                    cATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAg                      >  1:96372/1‑150 (MQ=255)
                                                                                                                      tAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCg                    >  2:1155549/1‑150 (MQ=255)
                                                                                                                          ccgcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGAGGTGACCCACAGCGACCa                >  1:135895/1‑150 (MQ=255)
                                                                                                                                 gccgcGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCaa         <  2:414788/150‑1 (MQ=255)
                                                                                                                                    gcgcACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAgcg      >  1:1203476/1‑150 (MQ=255)
                                                                                                                                     cgcctgCCTTTGCCTGCAAAGGCTCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAgcgc     <  2:1037351/145‑1 (MQ=255)
                                                                                                                                     cgcACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAgcgc     <  2:178053/150‑1 (MQ=255)
                                                                                                                                      gtgcGCCTTTGCAGGCAAAGGCGTGCGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAg        <  1:231331/143‑1 (MQ=255)
                                                                                                                                      gtgcGCCTTTGCAGGCAAAGGCGTGCGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAg        >  2:231331/4‑146 (MQ=255)
                                                                                                                                        gcGCCTTTGCAGGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTgg  >  2:1103272/2‑150 (MQ=255)
                                                                                                                                            cTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGAc                                                                                       >  2:435690/1‑61 (MQ=255)
                                                                                                                                            cTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGAc                                                                                       <  1:435690/61‑1 (MQ=255)
                                                                                                                                                  cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTca                                         <  2:538765/101‑1 (MQ=255)
                                                                                                                                                  cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTca                                         >  1:538765/1‑101 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
ATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGG  >  NC_002947/4480464‑4480749

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