Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,248,223 A→G 6.0% Y104H (TAC→CAC)  lysP ← lysine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,248,2230AG 6.0% 202.7 / 2.5 66Y104H (TAC→CAC) lysPlysine transporter
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base G (2/2);  total (34/32)
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.09e-01

CAAAGCCACGAACTGAGATGTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATAT  >  NC_000913/2248076‑2248369
                                                                                                                                                   |                                                                                                                                                  
cAAAGCCACGAACTGAGATGTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTAc                                                                                                                                                   >  8:109466/1‑149 (MQ=255)
  aaGCCACGAACTGAGATGTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTAc                                                                                                                                                   >  7:374550/1‑147 (MQ=255)
        cGAACTGAGATGTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGtt                                                                                                                                           <  7:4150/149‑1 (MQ=255)
        cGAACTGAGATGTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGtt                                                                                                                                           <  4:18167/149‑1 (MQ=255)
         gAACTGAGATGTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTAACCGCCCAGTTGTACCAGTAGTTc                                                                                                                                          <  7:283197/149‑1 (MQ=255)
          aaCTGAGATGTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTcc                                                                                                                                         <  7:304725/149‑1 (MQ=255)
             tGAGATGTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTccagc                                                                                                                                      >  3:353516/1‑149 (MQ=255)
               agaTGTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTccagcc                                                                                                                                     >  7:208931/1‑148 (MQ=255)
                   gTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAgcg                                                                                                                                >  2:224776/1‑149 (MQ=255)
                          agcagGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGaa                                                                                                                            <  3:22402/146‑1 (MQ=255)
                          agcagGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGaa                                                                                                                            >  4:22402/1‑146 (MQ=255)
                            cagGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGCGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGCAGCCa                                                                                                                        >  4:73574/1‑148 (MQ=255)
                               gAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTACCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGc                                                                                                                    >  1:235174/1‑149 (MQ=255)
                                aaGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGcc                                                                                                                   <  1:164382/149‑1 (MQ=255)
                                   aTAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCCCCCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCcttc                                                                                                                >  8:134470/1‑149 (MQ=255)
                                   aTAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCcttc                                                                                                                <  4:38543/149‑1 (MQ=255)
                                           gAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACa                                                                                                         >  2:225225/1‑148 (MQ=255)
                                           gAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACa                                                                                                         >  2:218988/1‑148 (MQ=255)
                                            aGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATa                                                                                                       <  1:332995/149‑1 (MQ=255)
                                                       aCTCCAGATCCAGCCCGGTGTATCCTGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGt                                                                                            <  3:36999/149‑1 (MQ=255)
                                                                ccAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGa                                                                                                <  5:136714/136‑1 (MQ=255)
                                                                ccAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGa                                                                                                >  6:136714/1‑136 (MQ=255)
                                                                    cccGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTc                                                                                                                                          >  8:49540/1‑90 (MQ=255)
                                                                    cccGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTc                                                                                                                                          <  7:49540/90‑1 (MQ=255)
                                                                      cGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAAc                                                                             <  7:109466/149‑1 (MQ=255)
                                                                       ggTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAAcc                                                                            >  7:177702/1‑149 (MQ=255)
                                                                        gtgtATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCg                                                                           <  2:5522/149‑1 (MQ=255)
                                                                         tgtATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCg                                                                           >  3:157459/1‑148 (MQ=255)
                                                                            aTCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCATCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGaaa                                                                       >  8:402419/1‑149 (MQ=255)
                                                                             tCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGaaa                                                                       <  6:163308/148‑1 (MQ=255)
                                                                             tCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAAc                                                                      >  3:234429/1‑149 (MQ=255)
                                                                                 ggAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCgg                                                                   >  8:216946/1‑148 (MQ=255)
                                                                                 ggAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCgg                                                                   >  7:97332/1‑148 (MQ=255)
                                                                                  gAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGc                                                                  >  2:116198/1‑148 (MQ=255)
                                                                                  gAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCa                                                                 <  1:224776/149‑1 (MQ=255)
                                                                                       accaGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCatatat                                                            >  1:277579/1‑149 (MQ=255)
                                                                                        ccaGCTCATCACCAGCTTAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAATTGGCAAAAGAACCGGAAACCGGCatatata                                                           >  6:204425/1‑148 (MQ=255)
                                                                                         caGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTGCCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATg                                                           >  8:72436/1‑148 (MQ=255)
                                                                                         caGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTGCCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATg                                                           >  8:65688/1‑148 (MQ=255)
                                                                                         caGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCTTATATGc                                                          <  1:225225/149‑1 (MQ=255)
                                                                                         caGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGc                                                          <  1:218988/149‑1 (MQ=255)
                                                                                            cTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGc                                                       >  7:11388/1‑149 (MQ=255)
                                                                                            cTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGc                                                       <  8:208931/149‑1 (MQ=255)
                                                                                               aTGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCa                                                     <  2:277579/148‑1 (MQ=255)
                                                                                                        tgagctgCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCgag                                           >  4:291614/1‑149 (MQ=255)
                                                                                                              gCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGt                                     >  2:76635/1‑149 (MQ=255)
                                                                                                               cAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGt                                     >  7:26337/1‑148 (MQ=255)
                                                                                                                aaCCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAGACGAACCGGAAACCGGCATATATGCAGCCa                                                     <  7:312663/131‑1 (MQ=255)
                                                                                                                aaCCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAGACGAACCGGAAACCGGCATATATGCAGCCa                                                     >  8:312663/1‑131 (MQ=255)
                                                                                                                  ccAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGtcat                                  >  6:14877/1‑148 (MQ=255)
                                                                                                                   cAGGTCAACGGCGATAGTCACCGCCCAGTGGTACCAGTCGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTCAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGtcatca                                <  4:96421/149‑1 (MQ=255)
                                                                                                                         aaCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTGTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGACAGTTCACCGAGACTGGTCATCAGGAAg                           >  8:198519/1‑148 (MQ=255)
                                                                                                                         aaCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGt                          <  7:134470/149‑1 (MQ=255)
                                                                                                                             gCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAAc                      <  8:11388/149‑1 (MQ=255)
                                                                                                                                aTAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCAt                   >  1:364753/1‑149 (MQ=255)
                                                                                                                                  aGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATc                  <  4:157459/148‑1 (MQ=255)
                                                                                                                                   gTCACCGCCCAGTTGTACCAGTAGTTCCAGCCATGCGCGAAGCAAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGTAAACCGGCATATCTGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCa                 <  7:374087/148‑1 (MQ=255)
                                                                                                                                    tCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCTAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGCGACTGGTCATCAGGAAGTAAACCATCAgg               >  8:155039/1‑149 (MQ=255)
                                                                                                                                     cACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCa                                                         <  2:160819/106‑1 (MQ=255)
                                                                                                                                     cACCGCCCAGTTGTACCAGTAGTTACAGCCCAGCGCGAAGCCAACGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACTCACCGGAAACCGGAATATATGCAGCCAGTTCAACGAGAATGGTCATCAGGAAGTAAACAATAAGGc              >  2:311261/1‑149 (MQ=255)
                                                                                                                                      aCCGCCCAGTTGTGCCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGc              <  7:72436/148‑1 (MQ=255)
                                                                                                                                      aCCGCCCAGTTGTGCCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGc              <  7:65688/148‑1 (MQ=255)
                                                                                                                                           ccAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCa        <  8:97332/149‑1 (MQ=255)
                                                                                                                                               ttGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGTCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCat    <  1:116198/149‑1 (MQ=255)
                                                                                                                                               ttGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCat    <  4:353516/149‑1 (MQ=255)
                                                                                                                                                  tACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTATGACAGTAATTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGATAATGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCatat  <  1:311261/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAAAGCCACGAACTGAGATGTAGTTCAGCAGGAAGATAACGCCGAGGAACAACGCACTCCAGATCCAGCCCGGTGTATCCGGGAACCACCAGCTCATGACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATAT  >  NC_000913/2248076‑2248369

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: ‑