Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,137,911 A→T 5.7% intergenic (‑254/+21) 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,9110AT 5.7% 341.0 / 7.2 121intergenic (‑254/+21)yijF/gldADUF1287 family protein/glycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (48/64);  new base T (2/5);  total (50/71)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 122 total aligned reads displayed.

GGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTG  >  NC_000913/4137762‑4138059
                                                                                                                                                     |                                                                                                                                                    
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)
             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                                                                                                                                     <  2:1216003/150‑5 (MQ=18)
                 ggCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGt                                                                                                                                     >  2:1395832/1‑150 (MQ=255)
                  gCACAGAGTGCCGGATGCGGAGTTAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGAGTAGGCCTGATAAGGGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTa                                                                                                                                    <  2:1302616/150‑1 (MQ=255)
                    acaGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAgg                                                                                                                                  <  2:1667784/150‑1 (MQ=255)
                     caGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGt                                                                                                                                 >  2:1707123/1‑150 (MQ=255)
                     caGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGt                                                                                                                                 >  1:1591832/1‑150 (MQ=255)
                      agagTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGtt                                                                                                                                <  2:616761/150‑1 (MQ=255)
                      agagTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGtt                                                                                                                                <  2:1114258/150‑1 (MQ=255)
                            ccGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTccc                                                                                                                          >  2:1220999/1‑150 (MQ=255)
                              ggATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCAGGAGTTTGGAGTAGGTTATTCCCAc                                                                                                                        <  1:644252/150‑1 (MQ=255)
                                   cGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTg                                                                                                                   <  1:804017/150‑1 (MQ=255)
                                    gtaGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGGTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGc                                                                                                                  <  2:1138630/148‑1 (MQ=255)
                                          aaCGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGaaa                                                                                                            <  2:471595/150‑1 (MQ=255)
                                          aaCGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGaaa                                                                                                            >  1:75093/1‑150 (MQ=255)
                                          aaCGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGaaa                                                                                                            >  2:1403583/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)
                                               cTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTGCGTTCAAACTCCCGGACAAACAGGGAGTTTTGAGTATGTTATTCCCACTCTTGCAGTAAACGCTg                                                                                                       >  2:31605/1‑150 (MQ=25)
                                                 tATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGTCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGAc                                                                                                     <  1:961453/150‑1 (MQ=255)
                                                 tATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGAc                                                                                                     >  2:1703644/1‑150 (MQ=255)
                                                      ggCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTAc                                                                                                >  2:943189/1‑150 (MQ=255)
                                                      ggCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTAc                                                                                                >  1:1103536/1‑150 (MQ=255)
                                                         cTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTgg                                                                                             <  2:1123810/150‑1 (MQ=255)
                                                           aCAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTc                                                                                           >  2:306310/1‑150 (MQ=255)
                                                            cAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCg                                                                                          <  2:117963/150‑1 (MQ=255)
                                                                aGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTa                                                                                      >  2:513108/1‑150 (MQ=255)
                                                                 gCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCCATTTTGCGTTCAAACTCCCGGCCAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTAc                                                                                     >  1:833702/1‑150 (MQ=255)
                                                                 gCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTAc                                                                                     >  1:382772/1‑150 (MQ=255)
                                                                    cGCAAATTCAATTTATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGGATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGTGTAGGTTATTCCGACTCTTGCAGGAAACGCTGACCGTAATGGTCGGCTACcag                                                                                  >  2:864682/1‑150 (MQ=38)
                                                                      cAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagca                                                                                >  1:1429722/1‑150 (MQ=255)
                                                                      cAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagca                                                                                >  1:369575/1‑150 (MQ=255)
                                                                      cAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagca                                                                                >  1:1000085/1‑150 (MQ=255)
                                                                       aaaTTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagcag                                                                               <  2:948511/150‑1 (MQ=255)
                                                                       aaaTTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagcag                                                                               >  1:1693210/1‑150 (MQ=255)
                                                                         aTTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCagag                                                                             <  2:431967/150‑1 (MQ=255)
                                                                              atatatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgcggcg                                                                        <  1:1500782/150‑1 (MQ=255)
                                                                              atatatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgcggcg                                                                        >  1:1720708/1‑150 (MQ=255)
                                                                              atatatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgcggcg                                                                        >  2:1541191/1‑150 (MQ=255)
                                                                                atatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTa                                                                      <  1:1604940/150‑1 (MQ=255)
                                                                                  atTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTaaa                                                                    <  2:1159046/150‑1 (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)
                                                                                                gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTc                                                      <  1:1210669/150‑1 (MQ=255)
                                                                                                    gCCTGATAAGCTTAGCGCATCAGTCAATTTTGATTTAAAACTAAAGGACAAGCAGTGAGTTTTGAGTAGGTTATTAACACTCTTGCAGGAAACGCTGACCGTACTGGTAGTCTACCAGCAGAGCGGCGTAAAACTGATCTGGCGTAgcgc                                                  >  2:1245498/1‑150 (MQ=255)
                                                                                                     ccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCAACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTTGCGTCGcccc                                                 >  2:1273350/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                                              <  2:75093/150‑1 (MQ=255)
                                                                                                        gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                              <  2:94936/150‑1 (MQ=255)
                                                                                                        gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                              <  1:513108/150‑1 (MQ=255)
                                                                                                             gCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCTTACAAGCCGGGAGTTTGGAGTAGGTTATTCCCTCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatc                                         <  2:437966/150‑4 (MQ=255)
                                                                                                             gCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatc                                         <  2:1072330/150‑4 (MQ=255)
                                                                                                              cGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatca                                        <  2:930468/150‑5 (MQ=255)
                                                                                                              cGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGCTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatca                                        <  1:65447/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)
                                                                                                                  gcgcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgttg                                    <  1:943189/150‑1 (MQ=255)
                                                                                                                          ggCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGt                            >  1:51744/1‑150 (MQ=255)
                                                                                                                          ggCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGt                            >  1:1412627/1‑150 (MQ=255)
                                                                                                                           gCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTTTTCCCACTCTTGCAGGAAACGCTGCCCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGtt                           >  1:65691/1‑150 (MQ=255)
                                                                                                                            cAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGttt                          >  1:1731118/1‑150 (MQ=255)
                                                                                                                                  tGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCAGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCtt                    >  2:1408015/1‑150 (MQ=255)
                                                                                                                                  tGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCAGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCtt                    >  2:1401000/1‑150 (MQ=255)
                                                                                                                                    cGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCt                  >  2:1232603/1‑150 (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)
                                                                                                                                                    cAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTg  <  1:1395832/150‑1 (MQ=255)
                                                                                                                                                    cAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAGGTTGTGAATGGTTTCACCTACTGCACATGCCGCTTCTg  >  1:1453116/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCTG  >  NC_000913/4137762‑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: