Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 824,572 T→A 5.7% intergenic (‑75/‑58) ybhP ← / → ybhQ endo/exonuclease/phosphatase family protein/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913824,5720TA 5.7% 637.4 / 19.6 175intergenic (‑75/‑58)ybhP/ybhQendo/exonuclease/phosphatase family protein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (88/77);  new base A (7/3);  total (95/80)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

Only 100 of 181 total aligned reads displayed.

GCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTG  >  NC_000913/824430‑824719
                                                                                                                                              |                                                                                                                                                   
    tAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTaaa                                                                                                                                           <  3:1733513/149‑1 (MQ=255)
    tAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTaaa                                                                                                                                           <  2:1142591/149‑1 (MQ=255)
          cTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCgg                                                                                                                                     >  3:1020749/1‑149 (MQ=255)
                gAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTActct                                                                                                                               >  3:1117631/1‑149 (MQ=255)
                   tatTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGa                                                                                                                             <  4:896194/148‑1 (MQ=255)
                   tatTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCGTATCTCATTGAAATAGTGTAGTCGGCGTTACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGag                                                                                                                            <  2:233684/149‑1 (MQ=255)
                     ttGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagta                                                                                                                          >  4:847330/1‑149 (MQ=255)
                     ttGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagta                                                                                                                          >  4:1295768/1‑149 (MQ=255)
                       gATGGTGAGCACCTTGAACGAAAATTGCTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAtttt                                                                                                                                                  >  1:160669/1‑123 (MQ=255)
                               gCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATa                                                                                                                >  2:92415/1‑149 (MQ=255)
                                   cTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAAATTCCGTaa                                                                                                                                            >  2:1612007/1‑117 (MQ=255)
                                   cTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAAATTCCGTaa                                                                                                                                            <  1:1612007/117‑1 (MQ=255)
                                   cTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAAATTCCGTaa                                                                                                                                            <  1:1342272/117‑1 (MQ=255)
                                    ttGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttc                                                                                                           >  3:1676626/1‑149 (MQ=255)
                                      gAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttctt                                                                                                         <  4:735358/149‑1 (MQ=255)
                                        aCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                                                                                                       <  4:748181/149‑1 (MQ=255)
                                          gAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAgg                                                                                                     >  4:1480831/1‑149 (MQ=255)
                                           aaaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGga                                                                                                    >  4:1201301/1‑149 (MQ=255)
                                            aaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTaaa                                                                                                                                           >  1:901460/1‑109 (MQ=255)
                                                  ttgtGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAg                                                                                                      >  1:1154897/1‑140 (MQ=255)
                                                  ttgtGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgagag                                                                                                 <  4:1061282/145‑1 (MQ=255)
                                                      gtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATg                                                                                         <  3:1295768/149‑1 (MQ=255)
                                                       tttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaa                                                                                                                         >  3:1616490/1‑116 (MQ=255)
                                                       tttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaa                                                                                                                         <  4:1616490/116‑1 (MQ=255)
                                                       tttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGaa                                                                                               <  3:187722/142‑1 (MQ=255)
                                                       tttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGa                                                                                        >  3:1195084/1‑149 (MQ=255)
                                                       tttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGa                                                                                        >  2:1290622/1‑149 (MQ=255)
                                                                     aTTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACgt                                                                          >  4:276569/1‑149 (MQ=255)
                                                                      ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTaaa                                                                                                                                           >  2:1112508/1‑83 (MQ=255)
                                                                      ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTaaa                                                                                                                                           <  1:1112508/83‑1 (MQ=255)
                                                                      ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATTAAGTGGCAACAACgtg                                                                         >  3:988468/1‑149 (MQ=255)
                                                                      ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACgtg                                                                         >  1:451795/1‑149 (MQ=255)
                                                                           cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                                                                                                       <  3:513089/114‑1 (MQ=255)
                                                                           cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                                                                                                       >  2:1417143/1‑114 (MQ=255)
                                                                           cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgt                                                                    >  3:1733494/1‑149 (MQ=255)
                                                                            tGTCTTTGCCTCTTATCTCTTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgtg                                                                   >  3:1500907/1‑149 (MQ=255)
                                                                              tCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                                                                                                       <  3:401205/111‑1 (MQ=255)
                                                                              tCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTc                                                                 >  3:1060457/1‑149 (MQ=255)
                                                                                tttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGc                                                               >  2:830062/1‑149 (MQ=255)
                                                                                    ccTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAtt                                                                                                                                                    <  4:1428413/60‑1 (MQ=255)
                                                                                    ccTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAtt                                                                                                                                                    >  3:1428413/1‑60 (MQ=255)
                                                                                     ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                          >  3:907317/1‑149 (MQ=255)
                                                                                      tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGcaac                                                                               <  3:826824/127‑1 (MQ=255)
                                                                                      tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACggg                                                         >  4:1225834/1‑149 (MQ=255)
                                                                                       ctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACgtg                                                                         <  4:573244/132‑1 (MQ=255)
                                                                                           tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                          >  4:450957/1‑143 (MQ=255)
                                                                                           tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                          <  3:450957/143‑1 (MQ=255)
                                                                                           tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTaa                                                    <  1:1290622/149‑1 (MQ=255)
                                                                                             tcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGt                                                                                     >  1:1464195/1‑114 (MQ=255)
                                                                                               aTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATaaa                                                                                                              >  1:520898/1‑87 (MQ=255)
                                                                                               aTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGcc                                                                                            >  3:577109/1‑105 (MQ=255)
                                                                                               aTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                          >  4:705813/1‑139 (MQ=255)
                                                                                               aTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTg                                                >  1:1284710/1‑149 (MQ=255)
                                                                                               aTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTg                                                >  3:263325/1‑149 (MQ=255)
                                                                                                ttGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTc                                                                      <  2:752796/126‑1 (MQ=255)
                                                                                                ttGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGc                                               <  4:263325/149‑1 (MQ=255)
                                                                                                 tGAAATATTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCt                                              >  2:1026048/1‑149 (MQ=255)
                                                                                                    aaTAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGc                                           >  4:1154089/1‑149 (MQ=255)
                                                                                                       aGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTc                                                                 >  2:1660155/1‑124 (MQ=255)
                                                                                                           tAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTagag                                                                                                                   >  2:1811732/1‑70 (MQ=255)
                                                                                                           tAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTc                                                                 <  3:676481/120‑1 (MQ=255)
                                                                                                           tAGTCGGCGTCACAAAAAGGTGCAGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTagag                                                                                                                   <  1:1811732/70‑1 (MQ=255)
                                                                                                            aGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATg                                   <  4:589301/149‑1 (MQ=255)
                                                                                                             gTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGt                                  <  1:1026048/149‑1 (MQ=255)
                                                                                                             gTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGt                                  <  2:1574072/149‑1 (MQ=255)
                                                                                                              tCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGtt                                 >  1:700834/1‑149 (MQ=255)
                                                                                                              tCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGtt                                 <  3:1154089/149‑1 (MQ=255)
                                                                                                              tCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGtt                                 <  3:1225834/149‑1 (MQ=255)
                                                                                                                  cGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATTAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCAt                             <  4:988468/149‑1 (MQ=255)
                                                                                                                      acaAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTAc                         <  4:817052/149‑1 (MQ=255)
                                                                                                                       caAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACt                        <  3:81334/149‑1 (MQ=255)
                                                                                                                       caAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACt                        >  2:1562282/1‑149 (MQ=255)
                                                                                                                        aaaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGAt                                                                                                                 <  2:1271198/59‑1 (MQ=255)
                                                                                                                        aaaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTg                       >  3:1344976/1‑149 (MQ=255)
                                                                                                                        aaaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTg                       >  1:1660722/1‑149 (MQ=255)
                                                                                                                        aaaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTg                       <  2:451795/149‑1 (MQ=255)
                                                                                                                        aaaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTg                       >  2:540255/1‑149 (MQ=255)
                                                                                                                               tGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTgg                >  4:127479/1‑149 (MQ=255)
                                                                                                                               tGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTgg                <  4:1733494/149‑1 (MQ=255)
                                                                                                                                 cGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGgcg              >  3:712101/1‑149 (MQ=255)
                                                                                                                                   gTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGa                                                                                                                             >  2:340237/1‑36 (MQ=255)
                                                                                                                                   gTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGa                                                                                                                             <  1:340237/36‑1 (MQ=255)
                                                                                                                                   gTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgct            <  1:1562282/149‑1 (MQ=255)
                                                                                                                                    tCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctg           >  4:1698117/1‑149 (MQ=255)
                                                                                                                                     tttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATg                                   >  3:882057/2‑124 (MQ=255)
                                                                                                                                     tttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgc          >  4:1415396/2‑149 (MQ=255)
                                                                                                                                     tttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgc          >  1:1020676/2‑149 (MQ=255)
                                                                                                                                     tttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgc          >  4:1572612/2‑149 (MQ=255)
                                                                                                                                     cTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgc          <  3:727633/149‑1 (MQ=255)
                                                                                                                                      ttACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                          <  3:1464221/100‑1 (MQ=255)
                                                                                                                                      ttACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgct         >  4:1027386/1‑149 (MQ=255)
                                                                                                                                        aCGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtgg       <  4:1344976/149‑1 (MQ=255)
                                                                                                                                          ggAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAg                   >  2:1145949/1‑135 (MQ=255)
                                                                                                                                          ggAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAg                   <  1:1145949/135‑1 (MQ=255)
                                                                                                                                          ggAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtg     >  3:1790607/1‑149 (MQ=255)
                                                                                                                                          ggAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtg     <  2:700834/149‑1 (MQ=255)
                                                                                                                                            aaTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtggt   >  3:67834/1‑149 (MQ=255)
                                                                                                                                            aaTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtggt   <  4:1864749/149‑1 (MQ=255)
                                                                                                                                             aTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtggtg  <  4:1676626/149‑1 (MQ=255)
                                                                                                                                              ttttCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGAtt                                      <  4:481371/112‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTG  >  NC_000913/824430‑824719

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: