Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,179,425 G→A 6.5% R278H (CGT→CAT)  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,4250GA 6.5% 187.8 / 3.1 61R278H (CGT→CAT) nagKN‑acetyl‑D‑glucosamine kinase
Reads supporting (aligned to +/- strand):  ref base G (37/20);  new base A (2/2);  total (39/22)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

ATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAACGCCGC  >  NC_000913/1179278‑1179573
                                                                                                                                                   |                                                                                                                                                    
aTCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCg                                                                                                                                                      <  1:169448/148‑1 (MQ=255)
    ggATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCTGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTc                                                                                                                                                 >  7:92079/1‑149 (MQ=255)
     gATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGtt                                                                                                                                                  <  7:22143/147‑1 (MQ=255)
           ttaGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCGGCCTCGTCATCTCGTCCCTGTAGCTCATGGTCCGCGCAt                                                                                                                                          >  4:148317/1‑149 (MQ=255)
             aGCGGTTTGTCTGGGAATTATCCTGACCATCGTTGACCCTGCCCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGGCATCTCATACCTGTAGCTCGTGTTCCTCGCAtt                                                                                                                                         >  3:110372/1‑148 (MQ=255)
              gCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTg                                                                                                                                        >  6:373371/1‑148 (MQ=255)
                ggTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAAc                                                                                                                                     <  8:92079/149‑1 (MQ=255)
                ggTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAAc                                                                                                                                     <  3:304278/149‑1 (MQ=255)
                   ttGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGGGTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAAcgcg                                                                                                                                  >  5:117327/1‑149 (MQ=255)
                      tCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTTGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAAcgcgcgc                                                                                                                               >  5:321568/1‑149 (MQ=255)
                         gggAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACgcgcgcgc                                                                                                                             >  7:8938/1‑148 (MQ=255)
                          ggAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACgcgcgcgc                                                                                                                             <  1:366877/147‑1 (MQ=255)
                                 tCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGAt                                                                                                                    >  5:9053/1‑149 (MQ=255)
                                    tGACCATCGTTGACCCTGCCCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCACGTCATCTCTTACCTGTAGCTCATGTTCCGAGCATTGAACGCGCGCGCCACGGTGATGCg                                                                                                                 >  8:78772/1‑149 (MQ=255)
                                    tGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCATACCTGTAGCTCGTGTTCCGCGCATTGACCGCGCGCGCCACGGTGATGCg                                                                                                                 >  6:18394/1‑149 (MQ=255)
                                       ccATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGgga                                                                                                              <  6:9053/149‑1 (MQ=255)
                                        cATCGGTGACCCTGACCTGGGCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCGGTAGCTCATGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGggag                                                                                                             <  3:148317/149‑1 (MQ=255)
                                             ttGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATg                                                                                                        >  3:153856/1‑149 (MQ=255)
                                                aCCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCgt                                                                                                     <  3:279169/149‑1 (MQ=255)
                                                aCCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCgt                                                                                                     <  4:270432/149‑1 (MQ=255)
                                                           gtcgtcATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGcc                                                                                            <  5:373371/147‑1 (MQ=255)
                                                            tcgtcATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTc                                                                                         <  8:64354/149‑1 (MQ=255)
                                                            tcgtcATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTc                                                                                         <  4:153856/149‑1 (MQ=255)
                                                              gtcATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCgtg                                                                                                    >  6:144992/1‑136 (MQ=255)
                                                              gtcATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCgtg                                                                                                    <  5:144992/136‑1 (MQ=255)
                                                              gtcATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGGCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCt                                                                                       <  1:296396/149‑1 (MQ=255)
                                                                cATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTa                                                                                      >  5:55474/1‑148 (MQ=255)
                                                                 aTTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCATGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTAc                                                                                     <  7:78772/148‑1 (MQ=255)
                                                                        gtggCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAAcc                                                                             >  4:93446/1‑149 (MQ=255)
                                                                             ttATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGAtt                                                                         >  2:285503/1‑148 (MQ=255)
                                                                              tATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTa                                                                        >  7:94887/1‑148 (MQ=255)
                                                                                tCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTaaa                                                                      <  8:8938/148‑1 (MQ=255)
                                                                                tCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGGCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTaaa                                                                      >  3:148974/1‑148 (MQ=255)
                                                                                   aaTTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAaacaa                                                                   >  4:368359/1‑148 (MQ=255)
                                                                                         ccGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGACCGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTGCCTACATCTAACCGATTAAACACCAGAGGt                                                            >  1:157220/1‑149 (MQ=255)
                                                                                            gCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGc                                                         <  3:93446/149‑1 (MQ=255)
                                                                                              aaTCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCCCGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTa                                                       >  7:332202/1‑149 (MQ=255)
                                                                                                     aCGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTc                                                >  3:299627/1‑149 (MQ=255)
                                                                                                         aaCTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGt                                            >  3:81545/1‑149 (MQ=255)
                                                                                                         aaCTGGCGGACAGGCGGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGGTGCTATGCTGTCGCGt                                            <  6:152932/149‑1 (MQ=255)
                                                                                                           cTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCg                                          >  6:260029/1‑149 (MQ=255)
                                                                                                              gCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCggg                                        <  8:94887/148‑1 (MQ=255)
                                                                                                                      gCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGt                               >  2:54229/1‑149 (MQ=255)
                                                                                                                      gCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGt                               <  1:54229/149‑1 (MQ=255)
                                                                                                                          ccTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAg                           >  8:9738/1‑149 (MQ=255)
                                                                                                                           cTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGt                          >  7:314797/1‑149 (MQ=255)
                                                                                                                           cTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGt                          >  2:47697/1‑149 (MQ=255)
                                                                                                                           cTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGAGCGCCACGGTGATGCGGCAGGAATGCGTGGTGCGGCATTCAAAAATCTACACGATTAAACAACAGAGGTTGCAATGCTGTCGCGACGGGGTCATCGGATAAGt                          >  2:111869/1‑149 (MQ=255)
                                                                                                                               tCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGtt                      >  1:289949/1‑149 (MQ=255)
                                                                                                                               tCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGtt                      >  3:128980/1‑149 (MQ=255)
                                                                                                                                cATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTCCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGttt                     >  1:290015/1‑149 (MQ=255)
                                                                                                                                 aTCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGtttt                    >  7:53332/1‑149 (MQ=255)
                                                                                                                                  tctcTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTc                   >  6:273294/1‑149 (MQ=255)
                                                                                                                                  tctcTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTc                   >  5:139639/1‑149 (MQ=255)
                                                                                                                                    tctTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGtttt                    <  7:234128/146‑1 (MQ=255)
                                                                                                                                    tctTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGtttt                    >  8:234128/1‑146 (MQ=255)
                                                                                                                                      ttACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCgg                                         <  4:293084/123‑1 (MQ=255)
                                                                                                                                          cTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTaaaaa            >  3:31662/1‑148 (MQ=255)
                                                                                                                                             tAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATaaa        >  1:22836/1‑149 (MQ=255)
                                                                                                                                                 tCGTGTTCCGCGCATTGAACGCGCGCTCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAAAAACAGAGGTTGCAATGCTGTCGCGGCGGGGTCATCGGTTAAGTCGTGTTCGTAACAATAAAcgcc    >  2:125582/1‑149 (MQ=255)
                                                                                                                                                 tCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAAcgcc    >  6:180520/1‑149 (MQ=255)
                                                                                                                                                 gcGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGGGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGGTGCTATGCTGGCGCGGCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAAcgcc    <  2:265707/148‑1 (MQ=255)
                                                                                                                                                   gtgtTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAAcgccgc  <  3:368359/149‑1 (MQ=255)
                                                                                                                                                   gtgtTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAAcgccgc  >  7:269477/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCTGGATTTATTAGCGGTTTGTCTGGGAAATATCCTGACCATCGTTGACCCTGACCTGGTCGTCATTGGTGGTGGCTTATCGAATTTCCCGGCAATCACAACGCAACTGGCGGACAGGCTGCCTCGTCATCTCTTACCTGTAGCTCGTGTTCCGCGCATTGAACGCGCGCGCCACGGTGATGCGGGAGGAATGCGTGGTGCGGCCTTCCTACATCTAACCGATTAAACAACAGAGGTTGCTATGCTGTCGCGTCGGGGTCATCGGTTAAGTCGTTTTCGTAAAAATAAACGCCGC  >  NC_000913/1179278‑1179573

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