Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,179,260 G→T 5.2% R223M (AGG→ATG)  nagK → N‑acetyl‑D‑glucosamine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,179,2600GT 5.2% 269.3 / 2.5 77R223M (AGG→ATG) nagKN‑acetyl‑D‑glucosamine kinase
Reads supporting (aligned to +/- strand):  ref base G (35/38);  new base T (2/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCA  >  NC_000913/1179115‑1179407
                                                                                                                                                 |                                                                                                                                                   
cgccgcTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAgg                                                                                                                                                    <  6:105766/147‑1 (MQ=255)
    gcTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCa                                                                                                                                               <  7:219855/148‑1 (MQ=255)
     cTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCAc                                                                                                                                              <  1:307232/148‑1 (MQ=255)
        cGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCTCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGtt                                                                                                                                           <  4:120605/148‑1 (MQ=255)
            tgtgGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGc                                                                                                                                       >  2:206216/1‑148 (MQ=255)
              tgGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGt                                                                                                                                     <  2:124979/148‑1 (MQ=255)
                        ggCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAATGGCGCACGTTGAGCGTTATCTGGAtt                                                                                                                           >  8:38816/1‑148 (MQ=255)
                        ggCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAATGGCGCACGTTGAGCGTTATCTGGAtt                                                                                                                           >  8:31886/1‑148 (MQ=255)
                         gCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGATGTATCAACACTATTATCATCACACGTTTCAGGCTCCCGAAATTATTGCGGTTTATGATAAAGGCGATGAGCAGGCAAGGGCTCACGTTGAGCGTTATCTGGAttt                                                                                                                          <  5:134065/148‑1 (MQ=255)
                         gCTGCATTGAAAATTATCTGTCTGCTGGCGGTTTTGGGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGCCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGGAAGGGCGCACGTTGCGCGTTATCTGGAttt                                                                                                                          >  6:134065/1‑148 (MQ=255)
                           tGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttat                                                                                                                        >  3:400218/1‑148 (MQ=255)
                            gCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGGGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatt                                                                                                                       >  8:426883/1‑148 (MQ=255)
                             cATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATttatta                                                                                                                      <  7:291751/148‑1 (MQ=255)
                                tGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCg                                                                                                                   >  1:277336/1‑148 (MQ=255)
                                  aaaaTTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGAttt                                                                                                                          >  1:331301/1‑139 (MQ=255)
                                  aaaaTTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGAttt                                                                                                                          <  2:331301/139‑1 (MQ=255)
                                     aTTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGATCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTg                                                                                                              <  4:129229/148‑1 (MQ=255)
                                        aTCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGGTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTc                                                                                                            <  6:379185/147‑1 (MQ=255)
                                             tctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCg                                                                                                                   <  1:270758/135‑1 (MQ=255)
                                             tctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCg                                                                                                                   >  2:270758/1‑135 (MQ=255)
                                             tctgGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTgg                                                                                                         <  2:48359/145‑1 (MQ=255)
                                                ggTCGCGGTTTTGCGTGACTGTATCAACACTATTATCAGCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCAGGTTGCGCGTTATCTGGATTTATTAGCGGGTTGTCTGGGAAata                                                                                                   >  6:109792/1‑148 (MQ=255)
                                                  tCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATc                                                                                                 <  7:91379/148‑1 (MQ=255)
                                                      ggTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGa                                                                                             >  6:95090/1‑148 (MQ=255)
                                                         tttGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCTCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATtta                                                                                                                         >  7:95607/1‑117 (MQ=255)
                                                         tttGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCTCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATtta                                                                                                                         <  8:95607/117‑1 (MQ=255)
                                                         tttGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCa                                                                                          >  2:166902/1‑148 (MQ=255)
                                                               tGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGGTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTg                                                                                    <  8:305377/148‑1 (MQ=255)
                                                                   tGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGAGCTTTATGATCAAGGCGATGAGCAGGCAAGCGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGAccc                                                                                <  7:426883/148‑1 (MQ=255)
                                                                      aTCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGa                                                                             >  1:138029/1‑148 (MQ=255)
                                                                        cAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGAcc                                                                           >  1:446420/1‑148 (MQ=255)
                                                                         aaCACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAATGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCt                                                                          <  7:31886/148‑1 (MQ=255)
                                                                         aaCACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAATGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCt                                                                          <  7:38816/148‑1 (MQ=255)
                                                                         aaCACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGAcc                                                                           <  2:277336/147‑1 (MQ=255)
                                                                          acacTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGAc                                                                                  <  3:207952/139‑1 (MQ=255)
                                                                          acacTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGAc                                                                                  >  4:207952/1‑139 (MQ=255)
                                                                          acacTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTg                                                                         >  3:70026/1‑148 (MQ=255)
                                                                            acTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTggg                                                                       <  1:251288/148‑2 (MQ=255)
                                                                               attatCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGTCGATGAGCAGGCAAGGGCGCACTTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGAGCCTGACCTGgtcgt                                                                    >  3:243200/1‑148 (MQ=255)
                                                                                  atcatcAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGa                                                                             >  3:178652/1‑136 (MQ=255)
                                                                                  atcatcAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGa                                                                             <  4:178652/136‑1 (MQ=255)
                                                                                  atcatcAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCAt                                                                 >  5:298577/1‑148 (MQ=255)
                                                                                    catcaACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGa                                                                             <  2:393845/134‑1 (MQ=255)
                                                                                    catcaACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGa                                                                             >  1:393845/1‑134 (MQ=255)
                                                                                       caACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATtggtg                                                            >  2:212883/1‑148 (MQ=255)
                                                                                        aaCCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATtggtgg                                                           <  2:446420/148‑1 (MQ=255)
                                                                                         aCCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAg                                                                                                                     <  5:249402/89‑1 (MQ=255)
                                                                                         aCCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAg                                                                                                                     >  6:249402/1‑89 (MQ=255)
                                                                                         aCCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTg                                                                                                              >  5:103981/1‑96 (MQ=255)
                                                                                         aCCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTg                                                                                                              <  6:103981/96‑1 (MQ=255)
                                                                                         aCCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATtggtggt                                                          >  8:177873/1‑148 (MQ=255)
                                                                                          ccGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCt                                                                               >  2:26215/1‑126 (MQ=255)
                                                                                          ccGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCt                                                                               <  1:26215/126‑1 (MQ=255)
                                                                                           cGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATtggtggtgg                                                        <  8:16388/148‑1 (MQ=255)
                                                                                              tGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATtta                                                                                                                         <  4:234549/80‑1 (MQ=255)
                                                                                              tGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATtta                                                                                                                         >  3:234549/1‑80 (MQ=255)
                                                                                                cAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTAt                                                   <  6:298577/148‑1 (MQ=255)
                                                                                                    cTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGtt                                                                                                                <  3:415049/83‑1 (MQ=255)
                                                                                                    cTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGtt                                                                                                                >  4:415049/1‑83 (MQ=255)
                                                                                                       ccGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAAttt                                            >  1:12755/1‑148 (MQ=255)
                                                                                                        cGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTc                                           >  6:432187/1‑148 (MQ=255)
                                                                                                          aaaTTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGt                                                                                                                 <  8:48559/76‑1 (MQ=255)
                                                                                                          aaaTTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGt                                                                                                                 >  7:48559/1‑76 (MQ=255)
                                                                                                           aaTTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATtggt                                                             <  3:161271/127‑1 (MQ=255)
                                                                                                           aaTTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATtggt                                                             >  4:161271/1‑127 (MQ=255)
                                                                                                           aaTTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGGCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCg                                        <  4:400218/148‑1 (MQ=255)
                                                                                                                ttGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAAt                                   <  5:95090/148‑1 (MQ=255)
                                                                                                                ttGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAAt                                   <  7:177873/148‑1 (MQ=255)
                                                                                                                     cTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACaa                              <  1:206216/148‑1 (MQ=255)
                                                                                                                          tGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCaa                         <  1:166902/148‑1 (MQ=255)
                                                                                                                            aTCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGa                                                                             >  2:73400/1‑94 (MQ=255)
                                                                                                                            aTCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGa                                                                             <  1:73400/94‑1 (MQ=255)
                                                                                                                             tCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTg                      >  7:401013/1‑148 (MQ=255)
                                                                                                                                aGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCg                   >  6:378529/1‑148 (MQ=255)
                                                                                                                                   cGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGAc                <  2:12755/148‑1 (MQ=255)
                                                                                                                                     aTGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGa                 >  8:378858/1‑145 (MQ=255)
                                                                                                                                     aTGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGa                 <  7:378858/145‑1 (MQ=255)
                                                                                                                                            ggCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCt       >  1:347144/1‑148 (MQ=255)
                                                                                                                                                 gggCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCa  >  1:298320/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGCCGCTGCGGCTGTGGTCAGCATGGCTGCATTGAAAATTATCTGTCTGGTCGCGGTTTTGCGTGGCTGTATCAACACTATTATCATCAACCGTTGCAGGCTCCCGAAATTATTGCGCTTTATGATCAAGGCGATGAGCAGGCAAGGGCGCACGTTGAGCGTTATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCA  >  NC_000913/1179115‑1179407

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