Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,982 G→C 11.2% 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,9820GC11.2% 328.7 / 27.1 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 (46/41);  new base C (6/5);  total (52/46)
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.98e-01

Only 100 of 105 total aligned reads displayed.

TTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTC  >  NC_000913/3368837‑3369126
                                                                                                                                                 |                                                                                                                                                
tttGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtgt                                                                                                                                               <  1:349101/149‑1 (MQ=255)
     cTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtgt                                                                                                                                               <  4:1765521/144‑1 (MQ=255)
     cTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtgt                                                                                                                                               >  3:1765521/1‑144 (MQ=255)
     cTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCa                                                                                                                                          >  1:1670081/1‑149 (MQ=255)
      tAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAg                                                                                                                                         <  4:1066304/149‑1 (MQ=255)
          acCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATc                                                                                                                                     >  2:1327873/1‑149 (MQ=255)
                 tGGGAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCt                                                                                                                                    >  2:1047319/1‑143 (MQ=255)
                 tGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCt                                                                                                                                    <  1:1047319/143‑1 (MQ=255)
                    cAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa                                                                                                                                       <  2:413345/137‑1 (MQ=255)
                    cAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa                                                                                                                                       >  1:413345/1‑137 (MQ=255)
                     aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                          <  3:438034/149‑1 (MQ=255)
                      aGCGAATATAAGATGCAACTACCTGATGAGAAACGTCTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGTCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGa                                                                                                                                        <  4:1578372/134‑1 (MQ=255)
                      aGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGa                                                                                                                                        >  3:1578372/1‑134 (MQ=255)
                         gAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTaa                                                                                                                      <  1:1327873/149‑1 (MQ=255)
                             ataAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                          <  4:185891/141‑1 (MQ=255)
                             ataAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                          >  3:185891/1‑141 (MQ=255)
                             ataAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCATTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAAt                                                                                                                                      >  4:748989/1‑129 (MQ=255)
                             ataAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCATTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAAt                                                                                                                                      <  3:748989/129‑1 (MQ=255)
                               aaGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTa                                                                                                                <  1:1473487/149‑1 (MQ=255)
                               aaGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTa                                                                                                                >  2:1473487/1‑149 (MQ=255)
                                  aTGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCt                                                                                                                                             <  4:1851877/117‑1 (MQ=255)
                                  aTGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCt                                                                                                                                             >  3:1851877/1‑117 (MQ=255)
                                     cAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCAt                                                                                                          <  2:1670081/149‑1 (MQ=255)
                                      aaCTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCAtt                                                                                                         >  3:5078/1‑149 (MQ=255)
                                          aCCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAAt                                                                                                     >  2:1679527/1‑149 (MQ=255)
                                           ccTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgt                                                                                                                                                 <  4:358217/104‑1 (MQ=255)
                                           ccTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgt                                                                                                                                                 >  3:358217/1‑104 (MQ=255)
                                            cTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCt                                                                                                                                             <  3:145489/107‑1 (MQ=255)
                                            cTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCt                                                                                                                                             >  4:145489/1‑107 (MQ=255)
                                            cTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGaa                                                                                                                                       <  1:1489459/113‑1 (MQ=255)
                                            cTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGaa                                                                                                                                       >  2:1489459/1‑113 (MQ=255)
                                            cTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCa                                                                                                   >  3:1643608/1‑149 (MQ=255)
                                                tgaGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                          <  1:291021/122‑1 (MQ=255)
                                                tgaGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                          >  2:291021/1‑122 (MQ=255)
                                                 gagaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCtt                                                                                              >  3:613775/1‑149 (MQ=255)
                                                 gagaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa                                                                                                                                       >  4:913719/1‑108 (MQ=255)
                                                 gagaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa                                                                                                                                       <  3:913719/108‑1 (MQ=255)
                                                  agaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTa                                                                                             >  4:811775/1‑149 (MQ=255)
                                                   gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTaa                                                                                            >  3:868235/1‑149 (MQ=255)
                                                    aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGaa                                                                                                                                       <  1:1525270/105‑1 (MQ=255)
                                                    aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGaa                                                                                                                                       >  2:1525270/1‑105 (MQ=255)
                                                    aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                           <  4:5078/149‑1 (MQ=255)
                                                     aaCGACTCGCAGATGAACTCGCTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAcc                                                                                          >  3:469455/1‑149 (MQ=255)
                                                        gACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGc                                                                                                 <  2:852530/139‑1 (MQ=255)
                                                         aCTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTca                                                                                      <  4:868235/149‑1 (MQ=255)
                                                                 aTGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCt                                                                                               <  3:553822/132‑1 (MQ=255)
                                                                 aTGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCt                                                                                               >  4:553822/1‑132 (MQ=255)
                                                                   gAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTaaaa                                                                                                                              <  1:1700999/99‑1 (MQ=255)
                                                                   gAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTaaaa                                                                                                                              >  2:1700999/1‑99 (MQ=255)
                                                                     aCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                          <  4:613775/149‑1 (MQ=255)
                                                                     aCTCGTTCGAACACAGGCGGTGCTAGAGAAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa                                                                                                                                       <  4:1770016/88‑1 (MQ=255)
                                                                     aCTCGTTCGAACACAGGCGGTGCTAGAGAAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACACTGTCTGCAGaa                                                                                                                                       >  3:1770016/1‑88 (MQ=255)
                                                                      cTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                          >  4:1401198/1‑100 (MQ=255)
                                                                      cTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                          <  3:1401198/100‑1 (MQ=255)
                                                                       tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTca                                                                                      <  3:180488/135‑1 (MQ=255)
                                                                       tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTca                                                                                      >  4:180488/1‑135 (MQ=255)
                                                                       tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTa                                                                        >  2:1001751/1‑149 (MQ=255)
                                                                       tCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGCTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTa                                                                        >  3:1616762/1‑149 (MQ=255)
                                                                         gTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCt                                                                                                                                             <  3:340281/78‑1 (MQ=255)
                                                                         gTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCt                                                                                                                                             >  4:340281/1‑78 (MQ=255)
                                                                              aaCACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtctctc                                                                 >  2:1396160/1‑149 (MQ=255)
                                                                                 acaGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAg                                                              <  3:811775/149‑1 (MQ=255)
                                                                                       ggTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAAt                                                                                                     <  1:1281100/104‑1 (MQ=255)
                                                                                       ggTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCAGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAAt                                                                                                     >  2:1281100/1‑104 (MQ=255)
                                                                                             agagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcacca                                                                                   >  4:1070935/1‑116 (MQ=255)
                                                                                             agagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcacca                                                                                   <  3:1070935/116‑1 (MQ=255)
                                                                                                    ggCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                          >  1:998841/1‑70 (MQ=255)
                                                                                                    ggCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATcc                                           >  4:1324513/1‑149 (MQ=255)
                                                                                                      cTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGc                                                        >  3:500151/1‑134 (MQ=255)
                                                                                                      cTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGc                                                        <  4:500151/134‑1 (MQ=255)
                                                                                                            aCGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttt                                                                           >  1:337042/1‑109 (MQ=255)
                                                                                                            aCGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttt                                                                           <  2:337042/109‑1 (MQ=255)
                                                                                                            aCGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTAt                                   <  1:1396160/149‑1 (MQ=255)
                                                                                                                cGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCgg                               <  1:1679527/149‑1 (MQ=255)
                                                                                                                cGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCgg                               <  3:1324513/149‑1 (MQ=255)
                                                                                                                 gTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAcc                                                                                          <  3:695664/89‑1 (MQ=255)
                                                                                                                 gTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAcc                                                                                          >  4:695664/1‑89 (MQ=255)
                                                                                                                  ttAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTg                                         <  4:1345057/137‑1 (MQ=255)
                                                                                                                  ttAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTg                                         >  3:1345057/1‑137 (MQ=255)
                                                                                                                          aGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGc                     <  1:1001751/149‑1 (MQ=255)
                                                                                                                           gCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCa                    >  3:741798/1‑149 (MQ=255)
                                                                                                                           gCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCa                    <  4:741798/149‑1 (MQ=255)
                                                                                                                                 tCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAaca              >  1:1540852/1‑149 (MQ=255)
                                                                                                                                   gAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                          <  1:1711424/87‑1 (MQ=255)
                                                                                                                                   gAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                          >  2:1711424/1‑87 (MQ=255)
                                                                                                                                     aTTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGc          <  4:1616762/149‑1 (MQ=255)
                                                                                                                                     aTTCTGCAGACACTCTCAGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATGAAATCAGCCTTGACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGc          >  2:1727745/1‑149 (MQ=255)
                                                                                                                                      ttCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACACCCTTTCTAATCTCTCCAGGCTCGCCGACAACGCATc                                            <  3:327087/114‑1 (MQ=255)
                                                                                                                                      ttCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATc                                            >  4:327087/1‑114 (MQ=255)
                                                                                                                                      ttCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGt                              <  4:493102/128‑1 (MQ=255)
                                                                                                                                       tCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGa        >  1:1610188/1‑149 (MQ=255)
                                                                                                                                         tGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGc                                                                                                 <  3:978581/58‑1 (MQ=255)
                                                                                                                                         tGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAAc      >  3:639537/1‑149 (MQ=255)
                                                                                                                                         tGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAAc      >  1:695580/1‑149 (MQ=255)
                                                                                                                                             gACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTc  >  3:1229161/1‑149 (MQ=255)
                                                                                                                                             gACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTc  >  2:744408/1‑149 (MQ=255)
                                                                                                                                              aCAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATc                                                                                                    <  3:1737809/50‑1 (MQ=255)
                                                                                                                                              aCAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATc                                                                                                    >  4:1737809/1‑50 (MQ=255)
                                                                                                                                                 gtgtCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGaa       <  1:1759358/140‑1 (MQ=255)
                                                                                                                                                 gtgtCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGaa       >  2:1759358/1‑140 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTC  >  NC_000913/3368837‑3369126

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: