Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,982 G→C 6.1% intergenic (+28/+32) yhcG → / ← yhcH DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily/DUF386 family protein, cupin superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,9820GC 6.1% 346.7 / 11.5 98intergenic (+28/+32)yhcG/yhcHDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily/DUF386 family protein, cupin superfamily
Reads supporting (aligned to +/- strand):  ref base G (53/39);  new base C (4/2);  total (57/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

TTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCGG  >  NC_000913/3368837‑3369130
                                                                                                                                                 |                                                                                                                                                    
tttGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtgt                                                                                                                                                   <  3:117470/149‑1 (MQ=255)
 ttGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTc                                                                                                                                                  <  1:185468/149‑1 (MQ=255)
    ccTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGc                                                                                                                                               >  4:332812/1‑149 (MQ=255)
     cTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCa                                                                                                                                              <  1:168802/149‑1 (MQ=255)
     cTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCa                                                                                                                                              >  2:470088/1‑149 (MQ=255)
     cTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCa                                                                                                                                              >  2:569140/1‑149 (MQ=255)
       aaCACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGa                                                                                                                                            >  2:149970/1‑149 (MQ=255)
          acCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGGTGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGAAtt                                                                                                                                         <  3:254483/149‑2 (MQ=255)
          acCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATc                                                                                                                                         >  2:724969/1‑149 (MQ=255)
            ctcacTGGCAAGCGAATATAGGAAGCAACTACCTGATGAGAAACGACACGCAGATGAGCTCGTTCGAACACAGGCGTTGCTAGAGGAAGGCTAAAACCGCCGATAATTTCAGCCGTTTCGCATTCTGCAGACAGTGTCTGCAGAATCtt                                                                                                                                       <  4:337435/145‑1 (MQ=255)
            cGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCtt                                                                                                                                       <  1:724969/149‑1 (MQ=255)
             gTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCttt                                                                                                                                      >  2:69191/1‑149 (MQ=255)
               tCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTaa                                                                                                                                    <  4:26740/149‑1 (MQ=255)
                     aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTaaaa                                                                                                                                  <  1:449792/145‑1 (MQ=255)
                     aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTaaaa                                                                                                                                  >  2:449792/1‑145 (MQ=255)
                     aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                              >  3:746340/1‑149 (MQ=255)
                     aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                              <  2:160945/149‑1 (MQ=255)
                     aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                              >  1:1551/1‑149 (MQ=255)
                     aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                              >  1:868560/1‑149 (MQ=255)
                     aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                              <  1:403592/149‑1 (MQ=255)
                       gCGAATATAAGATGCCACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGTAGACAGTGTCTGCAGAATCTTTAAAACGAACt                                                                                                                            >  2:349920/1‑149 (MQ=255)
                       gCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTGGAGGAAGGCTATAGACGCCGTTAATTTTAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTCACAACGAACt                                                                                                                            >  4:419725/1‑149 (MQ=255)
                       gCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACt                                                                                                                            >  2:822178/1‑149 (MQ=255)
                       gCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACt                                                                                                                            >  3:704654/1‑149 (MQ=255)
                         gAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTaa                                                                                                                          <  4:513267/149‑1 (MQ=255)
                               aaGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTa                                                                                                                    >  3:441054/1‑149 (MQ=255)
                                aGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                                                                   <  1:285753/149‑1 (MQ=255)
                                aGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                                                                   <  2:317263/149‑1 (MQ=255)
                                 gATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAg                                                                                                                  >  3:859019/1‑149 (MQ=255)
                                    gCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCa                                                                                                               >  2:352043/1‑149 (MQ=255)
                                      aaCTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCAtt                                                                                                             >  4:752693/1‑149 (MQ=255)
                                      aaCTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCAtt                                                                                                             >  1:318896/1‑149 (MQ=255)
                                      aaCTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCa                                                                                                                                              <  1:140216/116‑1 (MQ=255)
                                      aaCTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCa                                                                                                                                              >  2:140216/1‑116 (MQ=255)
                                          aCCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAAt                                                                                                         >  1:782455/1‑149 (MQ=255)
                                            cTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCa                                                                                                       <  1:441684/149‑1 (MQ=255)
                                            cTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCa                                                                                                       <  3:677915/149‑1 (MQ=255)
                                                   gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                                                                   <  2:756973/130‑1 (MQ=255)
                                                   gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTaa                                                                                                                   >  1:756973/1‑130 (MQ=255)
                                                   gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTaa                                                                                                <  1:352043/149‑1 (MQ=255)
                                                   gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTaa                                                                                                <  3:43184/149‑1 (MQ=255)
                                                   gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTaa                                                                                                >  4:484661/1‑149 (MQ=255)
                                                   gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTaa                                                                                                >  4:520106/1‑149 (MQ=255)
                                                    aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                               <  2:253079/149‑1 (MQ=255)
                                                    aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                               <  4:147641/149‑1 (MQ=255)
                                                      aCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCt                                                                                             >  4:490834/1‑149 (MQ=255)
                                                         aCTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTca                                                                                          <  2:1551/149‑1 (MQ=255)
                                                           tCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcacc                                                                                        >  1:602633/1‑149 (MQ=255)
                                                             gCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTGTTAAGCCATTAAATCAGCCTTAACCTTcaccac                                                                                      >  4:583464/1‑149 (MQ=255)
                                                                  tGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCt                                                                                 >  3:375728/1‑149 (MQ=255)
                                                                    aaCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttt                                                                               <  3:332812/149‑1 (MQ=255)
                                                                     aCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAAATCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                              <  3:583464/149‑1 (MQ=255)
                                                                      cTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttttt                                                                             >  4:387709/1‑149 (MQ=255)
                                                                      cTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttttt                                                                             <  1:420465/149‑1 (MQ=255)
                                                                       tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTa                                                                            >  2:506027/1‑149 (MQ=255)
                                                                       tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTa                                                                            >  1:863500/1‑149 (MQ=255)
                                                                        cGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGTCAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTaa                                                                           <  3:18351/149‑1 (MQ=255)
                                                                        cGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTaa                                                                           <  2:670309/149‑1 (MQ=255)
                                                                        cGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTaa                                                                           <  1:867803/149‑1 (MQ=255)
                                                                         gTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAt                                                                          >  2:99570/1‑149 (MQ=255)
                                                                             gAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtctct                                                                      >  4:219288/1‑149 (MQ=255)
                                                                                 acaGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTTGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAg                                                                  >  4:206924/1‑149 (MQ=255)
                                                                                   aGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGc                                                                >  4:858420/1‑149 (MQ=255)
                                                                                     gCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTc                                                              <  3:520106/149‑1 (MQ=255)
                                                                                          gCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGa                                                         >  3:640995/1‑149 (MQ=255)
                                                                                            tAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACa                                                       >  2:662531/1‑149 (MQ=255)
                                                                                             agagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTaa                                                                                                                          >  2:486603/1‑81 (MQ=255)
                                                                                             agagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTaa                                                                                                                          <  1:486603/81‑1 (MQ=255)
                                                                                              gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttttt                                                                             >  4:466125/1‑125 (MQ=255)
                                                                                              gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttttt                                                                             <  3:466125/125‑1 (MQ=255)
                                                                                              gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAAc                                                     <  2:14489/149‑1 (MQ=255)
                                                                                               tggACGGTTATAGACGCCGATCATTTCTGCCGATTCGAATTCTGTAAACAGTGTCTGCAGCATCTCTAAAACGAACTAAGTTTTAAGACATAAAATCCGGCTTGACCTTCACCACAAGCTTTTTTATCTCTCCAGGCTGGCCGACAACg                                                    >  2:564753/2‑149 (MQ=255)
                                                                                               agGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACg                                                    >  3:447496/1‑149 (MQ=255)
                                                                                                      cTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTg                                             >  2:20774/1‑149 (MQ=255)
                                                                                                       tataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTgg                                            >  1:86135/1‑149 (MQ=255)
                                                                                                        ataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGt                                           <  4:519743/149‑1 (MQ=255)
                                                                                                          aGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGttt                                         >  3:194968/1‑149 (MQ=255)
                                                                                                           gACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTa                                        <  2:185770/149‑1 (MQ=255)
                                                                                                            aCGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTAt                                       <  2:698348/149‑1 (MQ=255)
                                                                                                             cGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATg                                      >  1:885952/1‑149 (MQ=255)
                                                                                                              gCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                               <  4:61097/91‑1 (MQ=255)
                                                                                                              gCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                               >  3:61097/1‑91 (MQ=255)
                                                                                                                  ttAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGtt                                 >  2:681659/1‑149 (MQ=255)
                                                                                                                      tttCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCAcc                             <  3:308768/149‑1 (MQ=255)
                                                                                                                       ttCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCt                            <  1:812935/149‑1 (MQ=255)
                                                                                                                        tCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTg                           >  1:884064/1‑149 (MQ=255)
                                                                                                                          aGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGc                         <  1:822178/149‑1 (MQ=255)
                                                                                                                             ggTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATa                      >  3:152133/1‑149 (MQ=255)
                                                                                                                               tttCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATaaa                    >  2:756129/1‑149 (MQ=255)
                                                                                                                                 tCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAaca                  >  4:689935/1‑149 (MQ=255)
                                                                                                                                  cGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAacac                 >  1:228217/1‑149 (MQ=255)
                                                                                                                                   gAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACg                <  3:747302/149‑1 (MQ=255)
                                                                                                                                   gAATTCTGCAGACAGGGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACg                >  2:459889/1‑149 (MQ=255)
                                                                                                                                   agaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACg                >  2:249528/3‑149 (MQ=255)
                                                                                                                                    gaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACgg               >  4:509533/2‑149 (MQ=255)
                                                                                                                                    gaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACgg               >  4:287256/2‑149 (MQ=255)
                                                                                                                                    aaTTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACgg               <  4:375728/149‑1 (MQ=255)
                                                                                                                                      ttCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAATTCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCg             >  2:568427/1‑149 (MQ=255)
                                                                                                                                               cAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTccc    <  2:782455/149‑1 (MQ=255)
                                                                                                                                                 gtgtCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCgg  >  3:382079/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCGG  >  NC_000913/3368837‑3369130

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