Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,480,610 C→A 7.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 7.5% 379.0 / 9.6 119A145A (GCC→GCA) PP_3972short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (54/55);  new base A (4/6);  total (58/61)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

ATCATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGC  >  NC_002947/4480461‑4480759
                                                                                                                                                     |                                                                                                                                                     
atcatcGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGcc                                                                                                                                                       >  1:1049464/1‑150 (MQ=255)
   atcGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGc                                                                                                                                                    <  2:1384380/150‑1 (MQ=255)
    tcGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCa                                                                                                                                                   >  2:914663/1‑150 (MQ=255)
         cccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgc                                                                                                                                              <  2:1281328/150‑1 (MQ=255)
         cccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgc                                                                                                                                              <  1:1844472/150‑1 (MQ=255)
         cccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgc                                                                                                                                              <  1:402267/150‑1 (MQ=255)
          ccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgcg                                                                                                                                             >  2:1645460/1‑150 (MQ=255)
          ccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAAGCCGCGCCGCCGCGCACGCCTTTACAGGCCAAGgcg                                                                                                                                             >  2:1061532/1‑150 (MQ=255)
           cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgcgc                                                                                                                                            <  1:2322612/150‑1 (MQ=255)
           cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgc                                                                                                                                              >  1:1043624/1‑148 (MQ=255)
           cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgc                                                                                                                                              <  2:1043624/148‑1 (MQ=255)
            ttCGAACACCGTAGCAGCGGCACAGAGAAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                           >  2:564554/1‑150 (MQ=255)
            tccGAACACCGTAGCAGCGGCACAGAGCAGGAAGCGCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgcgtg                                                                                                                                           <  1:1061532/148‑3 (MQ=255)
              cGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACg                                                                                                                                         >  2:2004793/1‑150 (MQ=255)
                   acCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTgg                                                                                                                                    >  2:706276/1‑150 (MQ=255)
                   acCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTgg                                                                                                                                    >  1:418201/1‑150 (MQ=255)
                    ccGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGt                                                                                                                                   <  2:2228284/150‑1 (MQ=255)
                        agcagcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCGGCGTTGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTc                                                                                                                               <  2:2132604/150‑1 (MQ=255)
                           agcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCa                                                                                                                            <  2:1520232/150‑1 (MQ=255)
                                cacaGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                                     <  2:1579196/136‑1 (MQ=255)
                                cacaGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                                     >  1:1579196/1‑136 (MQ=255)
                                 acaGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCcgc                                                                                                                      >  1:490314/1‑150 (MQ=255)
                                   agagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGa                                                                                                                    <  2:1208099/150‑1 (MQ=255)
                                     agCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATg                                                                                                                  >  1:1926494/1‑150 (MQ=255)
                                       cAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaa                                                                                                                >  2:1542723/1‑150 (MQ=255)
                                       cAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaa                                                                                                                >  1:1241221/1‑150 (MQ=255)
                                        aagaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaag                                                                                                               <  1:995300/150‑1 (MQ=255)
                                          gaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaagaa                                                                                                             <  2:959267/150‑1 (MQ=255)
                                          gaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaagaa                                                                                                             >  2:2303590/1‑150 (MQ=255)
                                           aagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaagaa                                                                                                             <  2:1471969/149‑1 (MQ=255)
                                           aagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGaagaag                                                                                                            >  2:1013914/1‑150 (MQ=255)
                                           aagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGaagaag                                                                                                            <  1:914663/150‑1 (MQ=255)
                                            agCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGc                                                                                                           <  2:920942/150‑1 (MQ=255)
                                               tCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAACTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCatt                                                                                                        >  2:1750296/1‑150 (MQ=255)
                                                  ccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                           >  1:516863/1‑112 (MQ=255)
                                                  ccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCattatt                                                                                                     <  1:600190/150‑1 (MQ=255)
                                                  ccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgcgt                                                                                                                                            <  2:1044119/111‑2 (MQ=255)
                                                  ccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgcgt                                                                                                                                            >  1:1044119/1‑110 (MQ=255)
                                                   cGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGt                                                                                                                                   <  2:563159/119‑1 (MQ=255)
                                                   cGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTAttt                                                                                                    <  1:547236/150‑1 (MQ=255)
                                                       ggggATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACg                                                                                                <  1:1415307/150‑1 (MQ=255)
                                                         ggATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGtt                                                                                              >  1:1774707/1‑150 (MQ=255)
                                                          gATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTc                                                                                             >  1:139905/1‑150 (MQ=255)
                                                              tCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                                     >  1:12245/1‑106 (MQ=255)
                                                              tCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                                     <  2:12245/106‑1 (MQ=255)
                                                               cccGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGAt                                                                                        <  1:1386912/150‑1 (MQ=255)
                                                               cccGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGAt                                                                                        <  1:1595262/150‑1 (MQ=255)
                                                                ccGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACTGTCTGGATg                                                                                       >  1:1572879/1‑150 (MQ=255)
                                                                ccGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATg                                                                                       >  1:1533421/1‑150 (MQ=255)
                                                                ccGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgcgtgc                                                                                                                                          >  1:1796724/1‑96 (MQ=255)
                                                                  gTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGaa                                                                                     <  1:2004793/150‑1 (MQ=255)
                                                                      gggggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCa                                                                                 <  2:436145/150‑1 (MQ=255)
                                                                         ggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCg                                                                              >  2:550498/1‑150 (MQ=255)
                                                                          gCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGCCGTTCTGGATGAATGCACCGa                                                                             <  2:1241221/150‑1 (MQ=255)
                                                                          gCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGa                                                                             <  2:970559/150‑1 (MQ=255)
                                                                            cGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGAtt                                                                           <  2:187468/150‑1 (MQ=255)
                                                                            cGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGAtt                                                                           <  2:1533421/150‑1 (MQ=255)
                                                                            cGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGAtt                                                                           <  1:1013914/150‑1 (MQ=255)
                                                                               tGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgtg                                                                        <  2:1926494/150‑1 (MQ=255)
                                                                               tGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgtg                                                                        <  1:1589691/150‑1 (MQ=255)
                                                                               tGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgtg                                                                        <  1:1791066/150‑1 (MQ=255)
                                                                               tGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgtg                                                                        <  2:139905/150‑1 (MQ=255)
                                                                                 cGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGt                                                                      >  1:1147484/1‑150 (MQ=255)
                                                                                      cAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGa                                                                 >  2:2159571/1‑150 (MQ=255)
                                                                                       aaaGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAAt                                                                                    >  1:838630/1‑130 (MQ=255)
                                                                                       aaaGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAAt                                                                                    <  2:838630/130‑1 (MQ=255)
                                                                                         aGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGc                                                              <  2:1147484/150‑1 (MQ=255)
                                                                                             gtgGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATa                                                          <  1:706276/150‑1 (MQ=255)
                                                                                              tgGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAg                                                         >  2:169436/1‑150 (MQ=255)
                                                                                                   aTCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTc                                                    >  2:2121510/1‑150 (MQ=255)
                                                                                                    tCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTAttt                                                                                                    >  2:481789/1‑101 (MQ=255)
                                                                                                    tCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTAttt                                                                                                    <  1:481789/101‑1 (MQ=255)
                                                                                                    tCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTca                                                   >  2:1126426/1‑150 (MQ=255)
                                                                                                      acaGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgc                                                                                                                                              >  2:1914222/1‑57 (MQ=255)
                                                                                                      acaGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgc                                                                                                                                              <  1:1914222/57‑1 (MQ=255)
                                                                                                      acaGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTcacc                                                 >  2:1491067/1‑150 (MQ=255)
                                                                                                      acaGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTcacc                                                 >  1:1417578/1‑150 (MQ=255)
                                                                                                          gggCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGCTGCCATAGCGGTCACCACCg                                             >  2:86010/1‑150 (MQ=255)
                                                                                                          gggCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTc                                                    <  2:1080834/143‑1 (MQ=255)
                                                                                                              aTCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGGTGAATGCACCGATTGTGGTACGGATGCCATCGCGGGCACCACCGACGt                                         >  1:567883/1‑150 (MQ=255)
                                                                                                                ggcgCGGCCTATGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTAc                                                                    <  2:1682888/118‑1 (MQ=255)
                                                                                                                         tAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCg                              <  1:2159571/150‑1 (MQ=255)
                                                                                                                            gccgcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCACGATGAAGAAGCATTAATGGACTTTCTGTATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGAcc                           >  1:294195/1‑150 (MQ=255)
                                                                                                                              cgcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAg                         <  1:564554/150‑1 (MQ=255)
                                                                                                                              cgcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCCATCCCCACCGACGTCACCCACAGCGACCAg                         <  1:1172278/150‑1 (MQ=255)
                                                                                                                              cgcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGCATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAg                         >  1:1599668/1‑150 (MQ=255)
                                                                                                                               gcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAgg                        <  1:2121510/150‑1 (MQ=255)
                                                                                                                                cgccgccgcGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGt                       <  1:1645460/150‑1 (MQ=255)
                                                                                                                                    gccgcGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGAACAGGTTCaa                   >  1:2194107/1‑150 (MQ=255)
                                                                                                                                     ccgcGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAg                  >  1:2096068/1‑150 (MQ=255)
                                                                                                                                           aCGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTgg            >  2:1038964/1‑150 (MQ=255)
                                                                                                                                           aCGCCTGTGCAGGCAAGGGCGACCGCTTGATCCTCGCAGCCGGCGATGAAGAAGCATTATTTGTCGTTCTGGATGAATGCACCGATTGTGATACTGATGCCATAGCGGTCACCAACGACGTCACCCACAGCGACCAGGTTCAAGCGCTgg            <  2:567883/150‑1 (MQ=255)
                                                                                                                                            cGCCTTTGCAGGCAAAGGCGTGCGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgc           <  2:1812123/150‑1 (MQ=255)
                                                                                                                                             gCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgcc          >  2:350802/1‑150 (MQ=255)
                                                                                                                                              ccTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCAGACATGCAGAATCATTATTGGACGTTCTGGAGTAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccg         >  1:1854972/1‑150 (MQ=255)
                                                                                                                                                   gCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAg    <  1:598413/150‑1 (MQ=255)
                                                                                                                                                    ccTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTAGCTGAATGAACGGATGGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAgg   <  1:1495597/150‑1 (MQ=255)
                                                                                                                                                     cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGAc                                                                                                 >  1:865326/1‑55 (MQ=255)
                                                                                                                                                     cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGAc                                                                                                 <  2:865326/55‑1 (MQ=255)
                                                                                                                                                     cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCACGTTCAAGCGCTGGCCGCCCAGGc  <  1:169436/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGC  >  NC_002947/4480461‑4480759

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