Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,137,909 A→C 6.4% intergenic (‑252/+23) 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,9090AC 6.4% 335.8 / 8.2 118intergenic (‑252/+23)yijF/gldADUF1287 family protein/glycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (46/63);  new base C (3/5);  total (49/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

GTGGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCT  >  NC_000913/4137760‑4138058
                                                                                                                                                     |                                                                                                                                                     
gTGGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGa                                                                                                                                                       >  2:1059739/1‑150 (MQ=255)
 ggggAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGAc                                                                                                                                                      <  2:454786/149‑1 (MQ=255)
  gggAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGcaca                                                                                                                                                     >  2:1445475/1‑150 (MQ=255)
        ttCAGATTACTGGCACAGAGTGCCGGATGCGTAGTGAACGCCTTATCCGGACTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCgg                                                                                                                                               <  2:264869/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)
                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)
                    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:1114258/150‑1 (MQ=255)
                        agagTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGtt                                                                                                                               <  2:616761/150‑1 (MQ=255)
                              ccGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTccc                                                                                                                         >  2:1220999/1‑150 (MQ=255)
                                ggATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCAGGAGTTTGGAGTAGGTTATTCCCAc                                                                                                                       <  1:644252/150‑1 (MQ=255)
                                 gATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCAct                                                                                                                      <  2:1017680/150‑1 (MQ=255)
                                     cGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTg                                                                                                                  <  1:804017/150‑1 (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)
                                                 cTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTGCGTTCAAACTCCCGGACAAACAGGGAGTTTTGAGTATGTTATTCCCACTCTTGCAGTAAACGCTg                                                                                                      >  2:31605/1‑150 (MQ=25)
                                                  ttATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGa                                                                                                     <  1:1406010/150‑1 (MQ=255)
                                                   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)
                                                         gCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACt                                                                                              <  2:110572/150‑1 (MQ=255)
                                                         gCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACt                                                                                              <  1:833990/150‑1 (MQ=255)
                                                          ccTACATAAGCAAGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGGGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGttgat                                                                                             <  1:151343/150‑5 (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)
                                                                   gCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTAc                                                                                    >  1:382772/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                                                                              <  2:948511/150‑1 (MQ=255)
                                                                         aaaTTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagcag                                                                              >  1:1693210/1‑150 (MQ=255)
                                                                           aTTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCagag                                                                            <  2:431967/150‑1 (MQ=255)
                                                                            ttCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgc                                                                           <  1:553127/150‑1 (MQ=255)
                                                                                atatatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgcggcg                                                                       >  2:1541191/1‑150 (MQ=255)
                                                                                atatatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgcggcg                                                                       <  1:1500782/150‑1 (MQ=255)
                                                                                  atatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTa                                                                     <  1:1604940/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)
                                                                                                       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:94936/150‑1 (MQ=255)
                                                                                                          gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                             <  2:75093/150‑1 (MQ=255)
                                                                                                               gCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCTTACAAGCCGGGAGTTTGGAGTAGGTTATTCCCTCTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGGCGCGACGCCAgatc                                        <  2:437966/150‑4 (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)
                                                                                                                   aGCGCATCAGGCAATTTTGCGTTCAAACTCCCAGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACAGGTCGGCTACCAGCAGAGCTGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgtt                                    >  2:563100/1‑150 (MQ=255)
                                                                                                                            ggCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGt                           >  1:51744/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)
                                                                                                                                        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)
                                                                                                                                           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)
                                                                                                                                                     aCAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCt  <  2:369575/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTTCT  >  NC_000913/4137760‑4138058

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: