Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,657,895 A→T 5.6% intergenic (+25/‑174) ynfD → / → ynfE DUF1161 family periplasmic protein/putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,657,8950AT 5.6% 456.1 / 13.0 124intergenic (+25/‑174)ynfD/ynfEDUF1161 family periplasmic protein/putative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (59/58);  new base T (3/4);  total (62/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 130 total aligned reads displayed.

TTCCCAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTAT  >  NC_000913/1657747‑1658041
                                                                                                                                                    |                                                                                                                                                  
ttCCCAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGata                                                                                                                                                    >  1:1494175/1‑149 (MQ=255)
ttCCCAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGata                                                                                                                                                    <  4:1216229/149‑1 (MQ=255)
 tCCCAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGatat                                                                                                                                                   <  1:553834/149‑1 (MQ=255)
    cAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTgtc                                                                                                                                                <  2:1338571/149‑3 (MQ=255)
     aGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGa                                                                                                                                                 >  2:825835/1‑147 (MQ=255)
     aGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGa                                                                                                                                                 <  1:825835/147‑1 (MQ=255)
     aGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACa                                                                                                                                               >  1:197839/1‑149 (MQ=255)
       gtcgtcGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCTGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACaaa                                                                                                                                             <  1:1350430/149‑1 (MQ=255)
       gtcgtcGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACaaa                                                                                                                                             <  2:197839/149‑1 (MQ=255)
            cgGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                                                                                                        <  2:1193544/149‑1 (MQ=255)
                cATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAt                                                                                                                                    <  2:81640/149‑1 (MQ=255)
                cATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAt                                                                                                                                    <  1:186999/149‑1 (MQ=255)
                  ttGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAttt                                                                                                                                  <  1:1447575/149‑1 (MQ=255)
                     cgcTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGaa                                                                                                                               <  2:1494175/149‑1 (MQ=255)
                       cTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                                                                                             >  2:1270275/1‑149 (MQ=255)
                       cTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                                                                                             >  2:1138032/1‑149 (MQ=255)
                            gATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                                                                      >  2:8382/1‑135 (MQ=255)
                              tACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                                                                                                      >  4:1521450/1‑149 (MQ=255)
                                 gCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATaa                                                                                                                   >  2:1762581/1‑149 (MQ=255)
                                 gCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATaa                                                                                                                   <  3:570698/149‑1 (MQ=255)
                                  cATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATaa                                                                                                                   >  4:640364/1‑148 (MQ=255)
                                      aaaTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGatat                                                                                                                                                   >  2:623427/1‑112 (MQ=255)
                                      aaaTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGatat                                                                                                                                                   <  1:623427/112‑1 (MQ=255)
                                      aaaTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGaa                                                                                                                               <  3:1730313/132‑1 (MQ=255)
                                      aaaTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCa                                                                                                              >  3:866653/1‑149 (MQ=255)
                                        aTTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAcc                                                                                                            <  4:866653/149‑1 (MQ=255)
                                         ttCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCt                                                                                                           >  4:1335190/1‑149 (MQ=255)
                                            tctATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCgg                                                                                                                                        >  4:582243/1‑117 (MQ=255)
                                            tctATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCgg                                                                                                                                        <  3:582243/117‑1 (MQ=255)
                                                taCCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt                                                                                                    >  4:143740/1‑149 (MQ=255)
                                                   ccGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc                                                                                                                                       >  2:1758351/1‑111 (MQ=255)
                                                   ccGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc                                                                                                                                       <  1:1758351/111‑1 (MQ=255)
                                                   ccGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGa                                                                                                                                                 <  3:906800/101‑1 (MQ=255)
                                                   ccGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGa                                                                                                                                                 >  4:906800/1‑101 (MQ=255)
                                                   ccGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGACTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCAt                                                                                                 >  2:1613287/1‑149 (MQ=255)
                                                        ccaccaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCt                                                                                            >  4:482596/1‑149 (MQ=255)
                                                          accaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc                                                                                          <  1:1138032/149‑1 (MQ=255)
                                                          accaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc                                                                                          <  1:1270275/149‑1 (MQ=255)
                                                                ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGa                                                                                    <  1:1762581/149‑1 (MQ=255)
                                                                      gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                                                                                                      >  3:336936/1‑109 (MQ=255)
                                                                      gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                                                                                                      <  4:336936/109‑1 (MQ=255)
                                                                      gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat                                                                                   >  4:752199/1‑144 (MQ=255)
                                                                      gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat                                                                                   <  3:609434/144‑1 (MQ=255)
                                                                      gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat                                                                                   <  3:752199/144‑1 (MQ=255)
                                                                      gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGa                                                                              >  1:974642/1‑149 (MQ=255)
                                                                      gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGa                                                                              >  3:245957/1‑149 (MQ=255)
                                                                      gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAAAACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGa                                                                              >  1:1214784/1‑149 (MQ=255)
                                                                       ttactgcggttcggcaggcgcgccatcctggctagactgcGCGGGAGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                             >  2:451475‑M2/41‑149 (MQ=255)
                                                                       tctcTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAt                                                                                       <  1:831347/139‑1 (MQ=255)
                                                                       tctcTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                             <  2:1214784/149‑1 (MQ=255)
                                                                        ctctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAGTGATATGGAtt                                                                            <  3:1521450/149‑1 (MQ=255)
                                                                         tctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa                                                                                        >  4:1520887/1‑136 (MQ=255)
                                                                         tctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa                                                                                        <  3:1520887/136‑1 (MQ=255)
                                                                         tctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATta                                                                           >  4:1845641/1‑149 (MQ=255)
                                                                          ctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaa                                                                          <  3:1335190/149‑1 (MQ=255)
                                                                           tGCTCCCGCGCAGTCTAGCCAGCATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATCTGGATtaat                                                                         <  1:1613287/149‑1 (MQ=255)
                                                                            gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaata                                                                        >  2:842857/1‑149 (MQ=255)
                                                                              tCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACCCATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATACCACCGTATTACTCATACCCTTCTATTGATATGGATTaataat                                                                      >  3:1229851/1‑149 (MQ=255)
                                                                              tCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                             >  4:419106/1‑142 (MQ=255)
                                                                              tCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTaataat                                                                      >  1:1736281/1‑149 (MQ=255)
                                                                              tCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTaataat                                                                      >  4:763846/1‑149 (MQ=255)
                                                                                ttcggcaggcgcgccatcctggctagactgcCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCt                                                                                         <  4:449083‑M2/97‑1 (MQ=255)
                                                                                ttcggcaggcgcgccatcctggctagactgcCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCt                                                                                         >  3:449083‑M2/32‑128 (MQ=255)
                                                                                ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTc                                                                    >  1:900996/1‑149 (MQ=255)
                                                                                  ggcaggcgcgccatcctggctagactgcgCGGGAGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTa                                                                 <  1:451475‑M2/120‑1 (MQ=255)
                                                                                   cgcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                             >  4:1717285/1‑137 (MQ=255)
                                                                                    gcaggcgcgccatcctggctagactgcCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAc                                                                                                             <  4:410076‑M2/77‑1 (MQ=255)
                                                                                    gcaggcgcgccatcctggctagactgcCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAc                                                                                                             >  3:410076‑M2/28‑104 (MQ=255)
                                                                                     CAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                                                                        <  4:469764‑M1/92‑1 (MQ=255)
                                                                                     CAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                                                                        >  3:469764‑M1/1‑92 (MQ=255)
                                                                                     CAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAc                                                               >  2:1687415‑M1/1‑149 (MQ=255)
                                                                                      AGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGa                                                                                                                                                 >  3:432879‑M1/1‑66 (MQ=255)
                                                                                       ggcgcgccatcctggctagactgcCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                                                                <  3:1631234‑M2/122‑1 (MQ=255)
                                                                                       ggcgcgccatcctggctagactgcCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                                                                >  4:1631234‑M2/25‑146 (MQ=255)
                                                                                          tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                             >  2:13327/1‑130 (MQ=255)
                                                                                          tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                             <  1:13327/130‑1 (MQ=255)
                                                                                           aGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaaa                                                         <  3:1845641/149‑1 (MQ=255)
                                                                                            gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                                                                <  4:884792/141‑1 (MQ=255)
                                                                                                ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGt                                                    <  4:1440891/149‑1 (MQ=255)
                                                                                                ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGt                                                    >  3:1440891/1‑149 (MQ=255)
                                                                                                       gcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTcc                                             <  3:143740/149‑1 (MQ=255)
                                                                                                         gcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCct                                           <  1:1100533/149‑1 (MQ=255)
                                                                                                             gcCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAt                                                                                                                                    >  2:183190‑M2/3‑56 (MQ=255)
                                                                                                             gcCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAt                                                                                                                                    <  1:183190‑M2/54‑1 (MQ=255)
                                                                                                              cCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTAtt                                      <  3:763846‑M2/148‑1 (MQ=255)
                                                                                                              cCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTAtt                                      <  1:842857‑M2/148‑1 (MQ=255)
                                                                                                                 aaCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGatat                                                                                 >  2:1764956/1‑103 (MQ=255)
                                                                                                                 aaCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGatat                                                                                 <  1:1764956/103‑1 (MQ=255)
                                                                                                                  aCCGCAGTGACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGt                                    >  3:1741021/1‑147 (MQ=255)
                                                                                                                  aCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGt                                    <  4:1741021/147‑1 (MQ=255)
                                                                                                                  aCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTa                                  <  2:900996/149‑1 (MQ=255)
                                                                                                                         tAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATCTGAACATACAACATAATCCCACCTTATTACTCa                                                                                                  <  3:98532/78‑1 (MQ=255)
                                                                                                                               aTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                             <  4:1775722/93‑1 (MQ=255)
                                                                                                                               aTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt                                                                             >  3:1775722/1‑93 (MQ=255)
                                                                                                                                ttATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATAtt                    >  1:732363/1‑149 (MQ=255)
                                                                                                                                   tCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACCACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTccc                                            <  2:1344431/122‑1 (MQ=255)
                                                                                                                                   tCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTccc                                            >  1:1344431/1‑122 (MQ=255)
                                                                                                                                      aCTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGatat                                                                                 <  3:790766/82‑1 (MQ=255)
                                                                                                                                             gCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCaa                 <  1:490913/139‑1 (MQ=255)
                                                                                                                                             gCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCaa                 >  2:490913/1‑139 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCCCAGGTCGTCGGCCATTGCGCTAATGATACGCATAAAATTCTCTATACCCGCACCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTAT  >  NC_000913/1657747‑1658041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: