Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,502,271 C→T 5.7% V739I (GTT→ATT)  fryA ← putative PTS enzyme: Hpr, enzyme I and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,502,2710CT 5.7% 227.6 / 2.4 71V739I (GTT→ATT) fryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base T (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATT  >  NC_000913/2502125‑2502415
                                                                                                                                                  |                                                                                                                                                
gTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACg                                                                                                                                                 <  7:76457/148‑1 (MQ=255)
    gCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCgg                                                                                                                                             <  1:300778/148‑1 (MQ=255)
      aCCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACAAACGCCggtgg                                                                                                                                          <  8:258297/149‑1 (MQ=255)
          aGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTaa                                                                                                                                       >  6:249144/1‑148 (MQ=255)
             gTCAGCATGATCACCAGTGCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAAt                                                                                                                                   <  7:271395/149‑1 (MQ=255)
                aGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAAtctc                                                                                                                                <  5:70649/149‑1 (MQ=255)
                   aTGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTc                                                                                                                             <  6:164991/149‑1 (MQ=255)
                   aTGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTc                                                                                                                             <  3:53342/149‑1 (MQ=255)
                    tGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTcc                                                                                                                            >  1:209216/1‑149 (MQ=255)
                     gATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTccc                                                                                                                           >  8:82643/1‑149 (MQ=255)
                       tCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGCTACTGGAATGACGGAGCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGGAACAATCTCTTCCCGc                                                                                                                         >  2:237798/1‑149 (MQ=255)
                           cAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGc                                                                                                                      >  5:33162/1‑148 (MQ=255)
                            agatCGACTTCGCCCATATCTGACTGCCAGCCACTCGGTGTCGCCAGCCGGGCAATGCTGATACTCGACTCACGCATCCACTGAGATTTGGTCTGCGGGATCTCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCa                                                                                                                    <  4:191715/146‑1 (MQ=255)
                              ttCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTaa                                                                                                                                       <  6:5940/128‑1 (MQ=255)
                              ttCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTaa                                                                                                                                       >  5:5940/1‑128 (MQ=255)
                              ttCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGa                                                                                                                  >  2:94846/1‑149 (MQ=255)
                              ttCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGAAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCTGTGGTAACAATCTCTTCCCGCTGCCAGa                                                                                                                  <  2:93871/149‑1 (MQ=255)
                                 gACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCggtg                                                                                                                                           >  3:48090/1‑121 (MQ=255)
                                 gACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCggtg                                                                                                                                           <  4:48090/121‑1 (MQ=255)
                                      cGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACAtcttc                                                                                                           <  2:189467/148‑1 (MQ=255)
                                         ccATTTCTGACTGCCAGCCTATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATCACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTcc                                                                                                       <  5:70755/149‑1 (MQ=255)
                                          cATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCa                                                                                                      <  3:211524/149‑1 (MQ=255)
                                                  aCTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCaaa                                                                                               <  5:249144/148‑1 (MQ=255)
                                                      ccagccaATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACgg                                                                                            <  8:30358/147‑1 (MQ=255)
                                                         gccaATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGtt                                                                                       <  6:33162/149‑1 (MQ=255)
                                                          ccaATCGGTTACGCCAGCCGGGCAATGCAGATACTGGAATGACGGATCCACAGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTATTCCCGCTGCCAGACATATTCTGCCAGCTCAAACGGATCtt                                                                                       >  6:105320/1‑148 (MQ=255)
                                                            aaTCGGTTTAGCCAGCAGGGCAATGCTGATACTGGAATGACTAAGCCACTGAGATTTGGTGTTCGGTCTCGCAACTCAAAAACCAAAGCCGGGGGTAACAAGCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCATACGGATCTTCAg                                                                                    >  4:209027/1‑149 (MQ=255)
                                                             aTCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGt                                                                                   <  8:89667/149‑1 (MQ=255)
                                                             aTCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGt                                                                                   <  1:242124/149‑1 (MQ=255)
                                                                      gCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGt                                                                           >  3:200043/1‑148 (MQ=255)
                                                                            cGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTTGTAACAATGTCTTCCCGCTGCCAGACAGCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGcc                                                                    <  8:329797/149‑1 (MQ=255)
                                                                                            tGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAg                                                     >  2:25726/1‑148 (MQ=255)
                                                                                              gAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGcc                                                                              >  7:100617/1‑121 (MQ=255)
                                                                                              gAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGcc                                                                              <  8:100617/121‑1 (MQ=255)
                                                                                                 tGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAATGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCACACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTg                                               >  6:167241/1‑149 (MQ=255)
                                                                                                 tGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAATGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTg                                               >  6:181219/1‑149 (MQ=255)
                                                                                                   aCGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGcaca                                                      >  4:323860/1‑140 (MQ=255)
                                                                                                   aCGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGcaca                                                      <  3:323860/140‑1 (MQ=255)
                                                                                                   aCGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGa                                             >  7:53645/1‑149 (MQ=255)
                                                                                                    cGGATCCACTGAGATTTTGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGAt                                            <  1:324144/149‑1 (MQ=255)
                                                                                                    cGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGAt                                            <  1:94846/149‑1 (MQ=255)
                                                                                                          cACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCt                                       <  7:82643/148‑1 (MQ=255)
                                                                                                           aCTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCTAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGc                                     >  7:38538/1‑149 (MQ=255)
                                                                                                           aCTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGc                                     >  6:298567/1‑149 (MQ=255)
                                                                                                            cTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGct                                    >  7:84096/1‑149 (MQ=255)
                                                                                                             tGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGctc                                   >  1:285810/1‑149 (MQ=255)
                                                                                                             tGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGctc                                   >  1:152478/1‑149 (MQ=255)
                                                                                                             tGAGATTTGGTGTGAGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGctc                                   >  5:324178/1‑149 (MQ=255)
                                                                                                               agaTTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCtt                                 >  6:1382/1‑149 (MQ=255)
                                                                                                                 aTTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCTAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTa                               <  8:38538/149‑1 (MQ=255)
                                                                                                                 aTTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTa                               >  4:195743/1‑149 (MQ=255)
                                                                                                                  tttGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTAt                              >  2:54747/1‑149 (MQ=255)
                                                                                                                    tGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTAtt                             <  2:285810/148‑1 (MQ=255)
                                                                                                                       tgtgCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCt                          >  7:283633/1‑148 (MQ=255)
                                                                                                                         tgCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCACTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTaa                        >  2:213964/1‑148 (MQ=255)
                                                                                                                            gggATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATc                    <  5:298567/149‑1 (MQ=255)
                                                                                                                             ggATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGcc                                        >  2:61004/1‑128 (MQ=255)
                                                                                                                             ggATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGcc                                        <  1:61004/128‑1 (MQ=255)
                                                                                                                             ggATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATcc                   >  6:121998/1‑149 (MQ=255)
                                                                                                                              gATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCt                  <  4:200043/149‑1 (MQ=255)
                                                                                                                                  gCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGAt               <  8:84096/148‑1 (MQ=255)
                                                                                                                                    tACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATcaa            >  8:204005/1‑149 (MQ=255)
                                                                                                                                    tACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATcaa            >  7:324014/1‑149 (MQ=255)
                                                                                                                                        ccAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACaaa        <  5:296643/149‑1 (MQ=255)
                                                                                                                                         cAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATcaaca          >  8:178905/1‑146 (MQ=255)
                                                                                                                                         cAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAg       <  8:324014/149‑1 (MQ=255)
                                                                                                                                          aaaaCCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATcaaca          >  1:193049/1‑145 (MQ=255)
                                                                                                                                          aaaaCCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATcaaca          >  7:163896/1‑145 (MQ=255)
                                                                                                                                            aaCCAACGCCGGTGGTAACAATCTATTACCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATCTTCAGTGCGCCCGTTAACGCCGAGGTTGCTGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGata    >  4:109985/1‑149 (MQ=255)
                                                                                                                                              ccAATGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATAtt  <  5:167241/149‑1 (MQ=255)
                                                                                                                                              ccAATGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATAtt  <  5:181219/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATGACGGATCCACTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATT  >  NC_000913/2502125‑2502415

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