Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,137,910 C→T 6.8% intergenic (‑253/+22) yijF ← / ← gldA DUF1287 family protein/glycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,137,9100CT 6.8% 369.0 / 10.0 118intergenic (‑253/+22)yijF/gldADUF1287 family protein/glycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (47/63);  new base T (2/6);  total (49/69)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 123 total aligned reads displayed.

TGGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTG  >  NC_000913/4137761‑4138059
                                                                                                                                                     |                                                                                                                                                     
ggggAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGAc                                                                                                                                                       <  2:454786/149‑1 (MQ=255)
 gggAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGcaca                                                                                                                                                      >  2:1445475/1‑150 (MQ=255)
       ttCAGATTACTGGCACAGAGTGCCGGATGCGTAGTGAACGCCTTATCCGGACTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCgg                                                                                                                                                <  2:264869/150‑1 (MQ=255)
        tCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGATCGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCggg                                                                                                                                               <  1:367248/150‑1 (MQ=255)
        tCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCggg                                                                                                                                               <  2:1061834/150‑1 (MQ=255)
        tCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCggg                                                                                                                                               <  1:361954/150‑1 (MQ=255)
        tCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCggg                                                                                                                                               >  2:1273143/1‑150 (MQ=255)
        tCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATAGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCggg                                                                                                                                               >  2:65447/1‑150 (MQ=255)
         cAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCCGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGa                                                                                                                                              >  1:1311467/1‑150 (MQ=255)
            aTTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGtt                                                                                                                                           <  2:1249231/150‑1 (MQ=255)
              tACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTg                                                                                                                                         <  2:1287694/150‑1 (MQ=255)
               aCTGGCACAGAGTGCCGGATGCGGAGTGAACTCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTgg                                                                                                                                        <  1:1119091/150‑1 (MQ=255)
               aCTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCAGGCCTACAAAAGCACGCAAATTCAATATATTTCAGAGATTGCGTAGGCCTGATAATCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTgg                                                                                                                                        <  2:58931/150‑1 (MQ=255)
                cTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGCTTGTCCGGGAGTTTgaa                                                                                                                                       <  2:1114902/150‑3 (MQ=18)
                 tgtcacaGAGTGCCGGATGCGGAGTGAACGCATTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGTTTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGCTTGTCCGGGAGTTTgaac                                                                                                                                      <  2:1083135/147‑4 (MQ=17)
                 tGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAg                                                                                                                                      <  1:44090/150‑1 (MQ=255)
                  ggCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGCTTGTCCGGGAGTTTgaacg                                                                                                                                     <  2:1398758/150‑5 (MQ=18)
                  ggCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGCTTGTCCGGGAGTTTgaacg                                                                                                                                     <  1:1666241/150‑5 (MQ=18)
                  ggCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGCTTGTCCGGGAGTTTgaacg                                                                                                                                     <  2:1216003/150‑5 (MQ=18)
                  ggCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGt                                                                                                                                     >  2:1395832/1‑150 (MQ=255)
                  ggCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGt                                                                                                                                     <  2:1683276/150‑1 (MQ=255)
                   gCACAGAGTGCCGGATGCGGAGTTAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGAGTAGGCCTGATAAGGGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTa                                                                                                                                    <  2:1302616/150‑1 (MQ=255)
                     acaGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAgg                                                                                                                                  <  2:1667784/150‑1 (MQ=255)
                      caGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGt                                                                                                                                 >  1:1591832/1‑150 (MQ=255)
                       agagTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGtt                                                                                                                                <  2:1114258/150‑1 (MQ=255)
                       agagTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGtt                                                                                                                                <  2:616761/150‑1 (MQ=255)
                               ggATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCAGGAGTTTGGAGTAGGTTATTCCCAc                                                                                                                        <  1:644252/150‑1 (MQ=255)
                                gATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCAct                                                                                                                       <  2:1017680/150‑1 (MQ=255)
                                    cGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTg                                                                                                                   >  2:733087/1‑150 (MQ=255)
                                    cGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTg                                                                                                                   <  1:804017/150‑1 (MQ=255)
                                        gTGAACGCCTTATCCCGCCTACAAAAGCACGCAAATTCTATATATTGCAGAGTTTGCTTAGGCCTGATAAGCGTTGCGCATCAGGCAATTTTGGGTTCTAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTGATGCAgg                                                                                                               <  2:391022/150‑1 (MQ=21)
                                           aaCGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGaaa                                                                                                            >  2:1403583/1‑150 (MQ=255)
                                           aaCGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGaaa                                                                                                            <  2:471595/150‑1 (MQ=255)
                                           aaCGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGaaa                                                                                                            >  1:75093/1‑150 (MQ=255)
                                            aCGCCTTATCCGGCCTACAAAAGCAAGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACATGACGGGAGTTTGGAGTATGTTATTAACACTCTTGCTGGAAAc                                                                                                           >  1:1065841/1‑150 (MQ=21)
                                               ccTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTATCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCt                                                                                                        <  1:56934/150‑1 (MQ=255)
                                               ccTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGGGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCt                                                                                                        <  1:900293/150‑1 (MQ=255)
                                                 ttATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGa                                                                                                      <  1:1406010/150‑1 (MQ=255)
                                                  tATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGTCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGAc                                                                                                     <  1:961453/150‑1 (MQ=255)
                                                  tATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGAc                                                                                                     >  2:1703644/1‑150 (MQ=255)
                                                       ggCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTAc                                                                                                >  1:1103536/1‑150 (MQ=255)
                                                       ggCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTAc                                                                                                >  2:943189/1‑150 (MQ=255)
                                                        gCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACt                                                                                               <  1:833990/150‑1 (MQ=255)
                                                        gCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACt                                                                                               <  2:110572/150‑1 (MQ=255)
                                                          cTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTgg                                                                                             <  2:1123810/150‑1 (MQ=255)
                                                            aCAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTc                                                                                           >  2:306310/1‑150 (MQ=255)
                                                             cAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCg                                                                                          <  2:117963/150‑1 (MQ=255)
                                                               aaaGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGc                                                                                        <  1:261584/150‑1 (MQ=255)
                                                                 aGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTa                                                                                      >  2:513108/1‑150 (MQ=255)
                                                                     cGCAAATTCAATTTATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGGATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGTGTAGGTTATTCCGACTCTTGCAGGAAACGCTGACCGTAATGGTCGGCTACcag                                                                                  >  2:864682/1‑150 (MQ=38)
                                                                       cAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagca                                                                                >  1:1000085/1‑150 (MQ=255)
                                                                       cAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagca                                                                                >  1:369575/1‑150 (MQ=255)
                                                                       cAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagca                                                                                >  1:1429722/1‑150 (MQ=255)
                                                                        aaaTTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagcag                                                                               >  1:1693210/1‑150 (MQ=255)
                                                                        aaaTTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagcag                                                                               <  2:948511/150‑1 (MQ=255)
                                                                          aTTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCagag                                                                             <  2:431967/150‑1 (MQ=255)
                                                                           ttCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgc                                                                            <  1:553127/150‑1 (MQ=255)
                                                                               atatatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgcggcg                                                                        >  1:1720708/1‑150 (MQ=255)
                                                                               atatatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgcggcg                                                                        <  1:1500782/150‑1 (MQ=255)
                                                                                   atTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTaaa                                                                    <  2:1162295/150‑1 (MQ=255)
                                                                                   atTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTaaa                                                                    <  2:1159046/150‑1 (MQ=255)
                                                                                           gaTTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTTTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCt                                                            >  1:409275/1‑150 (MQ=255)
                                                                                           gaTTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCt                                                            >  1:55707/1‑150 (MQ=255)
                                                                                             ttGAGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTgg                                                          <  1:15304/150‑1 (MQ=255)
                                                                                              tGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGc                                                         <  2:1608382/150‑1 (MQ=255)
                                                                                                cGTAGGCCTGTTAAGCGTAGCGCTTCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGt                                                       <  1:335025/150‑1 (MQ=255)
                                                                                                 gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTTATCTGGCGTc                                                      >  2:1357736/1‑150 (MQ=255)
                                                                                                 gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTc                                                      <  1:1210669/150‑1 (MQ=255)
                                                                                                     gCCTGATAAGCTTAGCGCATCAGTCAATTTTGATTTAAAACTAAAGGACAAGCAGTGAGTTTTGAGTAGGTTATTAACACTCTTGCAGGAAACGCTGACCGTACTGGTAGTCTACCAGCAGAGCGGCGTAAAACTGATCTGGCGTAgcgc                                                  >  2:1245498/1‑150 (MQ=255)
                                                                                                       cTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGcgccg                                                <  2:85711/150‑1 (MQ=255)
                                                                                                        tGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGcgccgc                                               >  1:1104874/1‑150 (MQ=255)
                                                                                                         gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                              <  1:1237878/150‑1 (MQ=255)
                                                                                                         gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                              <  1:513108/150‑1 (MQ=255)
                                                                                                         gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                              <  2:94936/150‑1 (MQ=255)
                                                                                                               cGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatca                                        <  2:930468/150‑5 (MQ=255)
                                                                                                                gaaGCGCATCAGGCAATTTTGCGGTCAAACTCCCGGACCCGCCGGGAGTTTGGAGTAGGGTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatcag                                       <  1:151818/148‑6 (MQ=255)
                                                                                                                gTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtg                                       <  2:641316/150‑1 (MQ=255)
                                                                                                                 tAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGATAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgt                                      <  2:756268/150‑1 (MQ=255)
                                                                                                                 tAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgt                                      <  1:1707123/150‑1 (MQ=255)
                                                                                                                 tAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatcagg                                      <  2:1720708/150‑7 (MQ=255)
                                                                                                                     gcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAACCGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGTCGTCGCGCCGCCAGGCATGTtgtg                                  >  2:281090/1‑150 (MQ=255)
                                                                                                                           ggCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGt                            >  1:1412627/1‑150 (MQ=255)
                                                                                                                            gCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTTTTCCCACTCTTGCAGGAAACGCTGCCCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGtt                           >  1:65691/1‑150 (MQ=255)
                                                                                                                                   tGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCAGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCtt                    >  2:1408015/1‑150 (MQ=255)
                                                                                                                                     cGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCt                  >  2:1232603/1‑150 (MQ=255)
                                                                                                                                       ttCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGc                <  1:674757/150‑1 (MQ=255)
                                                                                                                                       tccAAACTCCCGGCTTGTCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGGGAATGGTTTCACCtgctgc                >  1:90236/3‑150 (MQ=255)
                                                                                                                                         cAAACTCCCGGCTTGTCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGcac              >  2:253303/1‑150 (MQ=255)
                                                                                                                                         cAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGcac              <  1:775663/150‑1 (MQ=255)
                                                                                                                                          aaaCTGCCGGTCCAGCCGGGAGTTTGGAGTAGGTTATTCCCTCTCTTGCAGGAAACGCTGACCGTTCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTTATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGcaca             <  1:1586636/150‑1 (MQ=255)
                                                                                                                                          aaaCTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGcaca             >  1:454590/1‑150 (MQ=255)
                                                                                                                                           aaCTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACAt            >  2:1018786/1‑150 (MQ=255)
                                                                                                                                             cTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCGCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGc          >  1:889972/1‑150 (MQ=255)
                                                                                                                                             cTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAACCCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGc          >  2:843057/1‑150 (MQ=255)
                                                                                                                                                ccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCATACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGc       <  2:332129/150‑1 (MQ=255)
                                                                                                                                                 cGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCt      <  1:1273143/150‑1 (MQ=255)
                                                                                                                                                   gACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTc    <  2:55707/150‑1 (MQ=255)
                                                                                                                                                    aCAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCt   <  2:369575/150‑1 (MQ=255)
                                                                                                                                                     cAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTTCTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGACAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTg  <  2:1433329/150‑1 (MQ=255)
                                                                                                                                                     cAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAGGTTGTGAATGGTTTCACCTACTGCACATGCCGCTTCTg  >  1:1453116/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTG  >  NC_000913/4137761‑4138059

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