Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,454,128 G→A 6.6% R182H (CGT→CAT)  yjfF → putative sugar ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,454,1280GA 6.6% 212.6 / 3.9 78R182H (CGT→CAT) yjfFputative sugar ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base A (2/3);  total (38/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.15e-01

CTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTT  >  NC_000913/4453980‑4454274
                                                                                                                                                    |                                                                                                                                                  
cgcGTTTCGGAAGAGTCTATTCCGATCAGCCCTCCCATTTATGACACGCTCTCCAGCCTTGCGTGGAATATCCCTTGCGTCGGTCGCTTAAGTGCGATGGGACTGCTGCTGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCg                                                                                                                                                    <  2:644949/147‑1 (MQ=255)
  cGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTa                                                                                                                                                  <  1:187448/149‑1 (MQ=255)
    tttCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGGAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTAACGGCATATTCCTCGCGCATCGTAcc                                                                                                                                                <  1:375450/149‑1 (MQ=255)
     ttCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGCTGGGACTGCTGATGTTGGCGGTGGGGGGTATCGGCATATTCCTCGCGCATCGGAccc                                                                                                                                               >  1:346140/1‑149 (MQ=255)
        ggAAGAGTCGATTCCGATAAACCATCACATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACGGCTGAGGTTGGCGGTGGTGGTTATCGGCATATTACTCGGGCATCGTACCCGtt                                                                                                                                            >  1:19128/1‑149 (MQ=255)
              gTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTa                                                                                                                                      >  1:305361/1‑149 (MQ=255)
              gTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTAACGGCATATTCCTCGCGCATCGTACCCGTTTTGGTa                                                                                                                                      >  2:77882/1‑149 (MQ=255)
               tCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGCTGGCGGTGGTGGTTATCGGCATATTCCCCGCGCATCGTACCCGGTTTGGtca                                                                                                                                     >  2:461875/1‑147 (MQ=255)
                 gATTCCGATTAACCATCCCTTTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGCGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATc                                                                                                                                   <  2:148591/149‑1 (MQ=255)
                  aTTCCGCTAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAACATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGGCTGCTGATGTTGGCGGTGGTGGCTGTCGGCCTCTTCCTCCCGCCTCGTACCCGTTTCGGTACTCa                                                                                                                                  >  2:91278/1‑149 (MQ=255)
                    tCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGGATCGTACCCGTTTTGGGAATCAgg                                                                                                                                >  2:728456/1‑149 (MQ=255)
                      cGATCAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGta                                                                                                                              >  2:631511/1‑149 (MQ=255)
                       gATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGtat                                                                                                                             >  1:138504/1‑149 (MQ=255)
                       gATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGtat                                                                                                                             >  2:547752/1‑149 (MQ=255)
                            aCCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGcg                                                                                                                        >  1:625409/1‑148 (MQ=255)
                             ccATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCa                                                                                                                       <  2:548412/149‑1 (MQ=255)
                              cATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCAt                                                                                                                      >  1:707735/1‑149 (MQ=255)
                               aTCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCAtt                                                                                                                     <  1:400936/149‑1 (MQ=255)
                                   cATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCATACCCGTTTTGGTAATCAGGTATACGCCATTgg                                                                                                                   >  2:601203/1‑147 (MQ=255)
                                   cATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCATACCCGTTTTGGTAATCAGGTATACGCCATTgg                                                                                                                   >  2:601264/1‑147 (MQ=255)
                                   cATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCATACCCGTTTTGGTAATCAGGTATACGCCATTgg                                                                                                                   <  1:601264/147‑1 (MQ=255)
                                   cATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCATACCCGTTTTGGTAATCAGGTATACGCCATTgg                                                                                                                   <  1:601203/147‑1 (MQ=255)
                                   cATTTAGGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTggcg                                                                                                                 <  1:353639/149‑1 (MQ=255)
                                    ctttATGACGCGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTggcgg                                                                                                                <  2:637958/148‑1 (MQ=255)
                                       tacgACACGCTCTCAAGCCGTGCCTGGAAAGTCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCCTATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGgcaa                                                                                                             <  2:496024/146‑1 (MQ=255)
                                       tATGACACGCTCTCAAGCCGTGCGTGGAAAATCCCTGGCGGCGGGCGCTTAAGTGCGATGGGCCTGCTGATGTTGGCGGTGGTGGTTATCGGCCTATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATCCGCCATTGGCGgcaa                                                                                                             <  2:495775/149‑1 (MQ=255)
                                        aTGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGTATATTACTCGCGCATGGTACCCGTTTGGGTAGTCACGTACACGCCATTgccggtgac                                                                                                            >  2:391601/1‑145 (MQ=255)
                                        aTGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTGGGTAATCCGGTATAGGCGATTGGCGgcaag                                                                                                            >  2:391694/1‑148 (MQ=255)
                                          gACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGgcaacgc                                                                                                          <  1:448123/149‑1 (MQ=255)
                                            cacGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGgcaacgcaa                                                                                                        <  1:277899/149‑1 (MQ=255)
                                                 tctcAAGCCTTGCGTGGAAAATCCTTGGCGGCGGTCGCTTAAGTGCGATGGGGCTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCATACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCg                                                                                                   <  1:212329/149‑1 (MQ=255)
                                                    cAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCCACGTCGGCg                                                                                                >  1:32924/1‑149 (MQ=255)
                                                     taGCCTTGCGGGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGa                                                                                               <  1:446931/148‑1 (MQ=255)
                                                       gCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAAt                                                                                             <  2:183393/149‑1 (MQ=255)
                                                       gCCGTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAAt                                                                                             <  2:343839/149‑1 (MQ=255)
                                                        ccTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGTATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATc                                                                                            >  1:486909/1‑149 (MQ=255)
                                                             cGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATg                                                                                       <  2:534073/149‑1 (MQ=255)
                                                              gTGGAACATCCCTGGGGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATgg                                                                                      <  2:478608/149‑1 (MQ=255)
                                                                ggAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATgggg                                                                                    >  1:712502/1‑149 (MQ=255)
                                                                ggAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATgggg                                                                                    <  2:705580/149‑1 (MQ=255)
                                                                   aaaTCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGAtt                                                                                 >  1:7024/1‑149 (MQ=255)
                                                                    attgctgggcggcGGGCGCTTAGGTGCGCTGGGGCGGCTGATGTGGGCGGTGGTGGTTATAGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGAttt                                                                                <  2:34139/142‑1 (MQ=255)
                                                                         cTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGGTGGGGATTTCCACt                                                                           >  1:452872/1‑149 (MQ=255)
                                                                            gcggcgGCCGCTGAAGGGCGATGGGACTGTTGATGTTTGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCATGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCgc                                                                        <  2:18359/149‑1 (MQ=255)
                                                                              ggcggTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGTAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCgcag                                                                      <  1:540606/149‑1 (MQ=255)
                                                                               gcggTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCGTATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTCTACGCCATTGGCGGCAACGCAACGTGGGCGAATCTGATGGGGAGTTCCACTCgcggc                                                                     >  2:187941/1‑149 (MQ=255)
                                                                                 ggTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGcac                                                                   <  1:118812/149‑1 (MQ=255)
                                                                                        tAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTAtt                                                            >  1:242194/1‑149 (MQ=255)
                                                                                           gTGCGATGGGACTGCGGATGTTGGCAGTGGTGGATATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATCCGCCATTGGCGGCCACGCAACGTCGGCGAATCTGATGTGGATTTCCACTCGCAGCACCACTATTCGc                                                         <  1:458984/149‑1 (MQ=255)
                                                                                            tgagATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATAGGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCa                                                        <  2:122410/146‑1 (MQ=255)
                                                                                               gATGGGACGGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCAttt                                                     <  1:278539/149‑1 (MQ=255)
                                                                                                aTGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCACCGCCGGCGAATCTGATGGGGGTGTCCACTCGCAGCACCACTATTCGCATTta                                                    >  1:209565/1‑149 (MQ=255)
                                                                                                  gggACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGTTGGGGATTTCCACTCGCAGCACCACTATTCGCATTtact                                                  >  2:385244/1‑147 (MQ=255)
                                                                                                      cggcggATGTTGTCGGTGTGGTTTATCGGCATATTTCTCGCGCATCCTACCCGTTTTGGTAATCAGGTCTACGCCATTTGCGGCAACGCAACGTCGGCGAATCTGATCGTGATTTCCACCCGCAGCACCACTATTCGCATTTATATGct                                              <  2:467569/144‑1 (MQ=255)
                                                                                                           gATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGACTCTGTTCATGATGTCCCCTCGCAGCGCCACTGTTCGCATTTATATGCTCTCCa                                         >  1:545311/1‑149 (MQ=255)
                                                                                                           gATGTTGGCGGTGGTGCGTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCa                                         <  2:680430/149‑1 (MQ=255)
                                                                                                            aTGTTGGCGGTGGTGGTTATCGGCATATTACTCGCGCATCGTACCCGTTTTGGTAATCAGGTCTACACCCTTGGCGGCAACGCAACGTCGGTGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCAc                                        >  1:475592/1‑149 (MQ=255)
                                                                                                              gTTGGCGGTGTTGGTTTTCGGCACATGCGTGGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGTAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCg                                      <  2:567860/149‑1 (MQ=255)
                                                                                                               ttGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGACTCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCgg                                     >  1:596086/1‑149 (MQ=255)
                                                                                                                tGGCGTTGCTGTTTATGCGCATGTTCCTCGTGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCGACGTCGGCTAATCTGATGGGGATTGCCACTCGCAGCACCACTGTTCGTATTGTTATGCTCTCCACCGGa                                    <  2:612012/149‑1 (MQ=255)
                                                                                                                tGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGa                                    <  1:410086/149‑1 (MQ=255)
                                                                                                                  gCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACt                                  <  2:712502/149‑1 (MQ=255)
                                                                                                                   cGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTg                                 <  2:595959/149‑1 (MQ=255)
                                                                                                                    aggggtggtTAGCGTCATATTCCCCGCGCAGCGCACCCGGTTTGGTAATCAGGTCTACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTgg                                <  1:72094/146‑1 (MQ=255)
                                                                                                                          ggtTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCTTTTCTATGCTCTCCATCGGACTAGCTAcac                          >  1:256082/1‑147 (MQ=255)
                                                                                                                          ggtTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGc                          >  2:598781/1‑149 (MQ=255)
                                                                                                                           gtTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCGGCACCACTATTCGCAGTTATATGCTCTCCACCGGACGGGCACCGCt                         >  1:373381/1‑149 (MQ=255)
                                                                                                                             tATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTgg                       <  1:573135/149‑1 (MQ=255)
                                                                                                                               tCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTGTACGCCATTGGCGGCAACGCCACGTCGGCGAATCTGATGGGGATTTCCCCTCGCAGCACCACTCTTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCg                     >  2:312449/1‑149 (MQ=255)
                                                                                                                                   cATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAAAGCTGGCCGGGa                 >  2:618097/1‑149 (MQ=255)
                                                                                                                                   cATATCCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTACACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCACTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGa                 <  1:91278/149‑1 (MQ=255)
                                                                                                                                     tatTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGAtt               >  1:284269/1‑149 (MQ=255)
                                                                                                                                         ccTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCCACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTCTTCGCATTTATATGCTCTCCACCGGACGGGCAACGCTGGCGGGGTTGGtct           >  2:389235/1‑149 (MQ=255)
                                                                                                                                           tCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGAGTGgcttc         >  2:497841/1‑149 (MQ=255)
                                                                                                                                            cgcgCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGCCTGGCAACGCTGGCGGGGATTTgcttct        >  2:275888/1‑149 (MQ=255)
                                                                                                                                             gcgcATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTtctc       <  1:424206/149‑1 (MQ=255)
                                                                                                                                               gcATCGTCCCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGa     <  1:430231/149‑1 (MQ=255)
                                                                                                                                                cATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGAt    <  1:77882/149‑1 (MQ=255)
                                                                                                                                                  tCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCGACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGGCTGGCAACGCTGGCGGGGGTTGTCTTCTCggtgt  >  2:410306/1‑145 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTCGTTTCGGAAGAGTCGATTCCGATAAACCATCCCATTTATGACACGCTCTCAAGCCTTGCGTGGAAAATCCCTGGCGGCGGTCGCTTAAGTGCGATGGGACTGCTGATGTTGGCGGTGGTGGTTATCGGCATATTCCTCGCGCATCGTACCCGTTTTGGTAATCAGGTATACGCCATTGGCGGCAACGCAACGTCGGCGAATCTGATGGGGATTTCCACTCGCAGCACCACTATTCGCATTTATATGCTCTCCACCGGACTGGCAACGCTGGCGGGGATTGTCTTCTCGATTT  >  NC_000913/4453980‑4454274

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