Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,480,610 C→A 5.5% A145A (GCC→GCA)  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,6100CA 5.5% 226.1 / 2.3 73A145A (GCC→GCA) PP_3972short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base A (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

ATCATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGG  >  NC_002947/4480461‑4480749
                                                                                                                                                     |                                                                                                                                           
atcatcGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGcc                                                                                                                                             <  2:984785/150‑1 (MQ=255)
   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                                                                                                                           <  1:970335/150‑1 (MQ=255)
                  cacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                           <  2:1201371/150‑1 (MQ=255)
                  cacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                           <  2:1201536/150‑1 (MQ=255)
                           agcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCa                                                                                                                  >  2:417021/1‑150 (MQ=255)
                           agcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCa                                                                                                                  >  2:170897/1‑150 (MQ=255)
                                   agagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGa                                                                                                          <  2:295780/150‑1 (MQ=255)
                                    gagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGAt                                                                                                         >  1:1051493/1‑150 (MQ=255)
                                       cAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGc                                                                                                                >  2:480729/1‑140 (MQ=255)
                                       cAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGc                                                                                                                <  1:480729/140‑1 (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                                                                                                                    >  1:355466/1‑104 (MQ=255)
                                                                       ggggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCg                                                                                                                    <  2:355466/104‑1 (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                                                                                                 >  1:694808/1‑92 (MQ=255)
                                                                                                      acaGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGc                                                                                                 <  2:694808/92‑1 (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                                                  >  2:610989/1‑130 (MQ=255)
                                                                                                               tCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCa                                                  <  1:610989/130‑1 (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        >  2:231331/4‑146 (MQ=255)
                                                                                                                                         gtgcGCCTTTGCAGGCAAAGGCGTGCGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAg        <  1:231331/143‑1 (MQ=255)
                                                                                                                                           gcGCCTTTGCAGGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTgg  >  2:1103272/2‑150 (MQ=255)
                                                                                                                                               cTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGAc                                                                                       <  1:435690/61‑1 (MQ=255)
                                                                                                                                               cTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGAc                                                                                       >  2:435690/1‑61 (MQ=255)
                                                                                                                                                     cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTca                                         <  2:538765/101‑1 (MQ=255)
                                                                                                                                                     cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTca                                         >  1:538765/1‑101 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
ATCATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGG  >  NC_002947/4480461‑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: ‑