Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 475,972 C→G 6.4% intergenic (+21/+10) ybaY → / ← ybaZ outer membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913475,9720CG 6.4% 155.0 / 4.7 65intergenic (+21/+10)ybaY/ybaZouter membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base G (2/2);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

GCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAAC  >  NC_000913/475831‑476113
                                                                                                                                             |                                                                                                                                             
gCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCCCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATTTTAAGCGACGTTAACCCTCTCCGCCCGGagag                                                                                                                                        >  2:111558/1‑149 (MQ=255)
 ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                                                                       <  2:929024/149‑1 (MQ=255)
  cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                                                                      <  2:377651/149‑1 (MQ=255)
      aCAGCAAACCGCCGTGCCGGGTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAGGGTtaa                                                                                                                                  <  2:73941/149‑2 (MQ=255)
          cAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTagtagt                                                                                                                              >  2:384725/1‑149 (MQ=255)
            aaCCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTc                                                                                                                            <  2:372650/149‑1 (MQ=255)
            aaCCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTAGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTc                                                                                                                            <  2:4490/149‑1 (MQ=255)
                           tCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCTACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGc                                                                                                             <  2:1085785/149‑1 (MQ=255)
                             aGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGCATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGTGGGGTAGTAGTTCCAGCGATAACGCTGCaa                                                                                                           >  2:687238/1‑149 (MQ=255)
                             aGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCCGCTCCTACGCAATATTAAGCGCCGTTAACCCTCTCCGCACGGACAGGGTTAGTAGTTCCAGCGGTACCGCTGCaa                                                                                                           >  1:777044/1‑149 (MQ=255)
                              ggCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTAGGCAATTTTATGCGCCGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCaat                                                                                                          >  1:29361/1‑148 (MQ=255)
                               gCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGATCAAGCTCCTACGCAATATTAAGCGACGTTAACCCTATCCGCCCGGAGAGGGTTAGTAGTTACAGCCATAAAGCTGCAAGt                                                                                                         >  1:623791/1‑149 (MQ=255)
                                 cAGCGGTGGCGCAACGCCTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGGCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCg                                                                                                       >  1:914468/1‑149 (MQ=255)
                                  aGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGa                                                                                                      <  1:542049/149‑1 (MQ=255)
                                   gCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAt                                                                                                     <  2:870070/149‑1 (MQ=255)
                                    cGGTGGCGCAACGACTACCGTACCTTCGACATCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTAGGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGTTATCGCTGGaactactaa                                                                                                    <  1:707251/149‑8 (MQ=255)
                                         gcgcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGccc                                                                                               <  2:276017/149‑1 (MQ=255)
                                           gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                      >  2:392206/1‑108 (MQ=255)
                                           gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                      <  1:392206/108‑1 (MQ=255)
                                             aaCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTAAGCAATATTAAGCGACGTTAACCCTCAACGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCtt                                                                                           >  1:825416/1‑149 (MQ=255)
                                                    aCCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATa                                                                                    <  2:311596/149‑1 (MQ=255)
                                                     ccGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAc                                                                                   >  2:427576/1‑149 (MQ=255)
                                                      cGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTACGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCACGCTTCCCGATAcc                                                                                  >  2:17752/1‑149 (MQ=255)
                                                        tACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCAt                                                                                <  1:226455/149‑1 (MQ=255)
                                                          ccTTCGACTTCACCAACTCAGGTGAATCCGGCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGGCGATTTGCCCGCTTCCCGATACCATca                                                                              <  1:1122391/149‑1 (MQ=255)
                                                            ttCGACATCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcaca                                                                            <  2:262832/149‑1 (MQ=255)
                                                                  ttCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCt                                                                      >  1:569445/1‑149 (MQ=255)
                                                                    cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCTATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                                                                    >  1:892097/1‑149 (MQ=255)
                                                                    cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                                                                    <  1:682500/149‑1 (MQ=255)
                                                                    cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                                                                    <  1:1151642/149‑1 (MQ=255)
                                                                       cAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGACAg                                                                 >  1:1010522/1‑149 (MQ=255)
                                                                       cAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCACGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAg                                                                 >  1:436828/1‑149 (MQ=255)
                                                                        aaCTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTCAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGt                                                                <  2:165851/149‑1 (MQ=255)
                                                                           tCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACTCAATATTAAGCTACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAAt                                                             <  2:1083905/149‑1 (MQ=255)
                                                                            cAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAACATCACACCTTCTGCCAGTAATg                                                            <  1:1047352/149‑1 (MQ=255)
                                                                             aGGTGAATACGTCTTCGGCAGTTCCAGCTCCTACTCAATATTAAGCGACGTTAACCCTCTACGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGGCTATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGc                                                           <  2:524927/149‑1 (MQ=255)
                                                                                 gAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGt                                                       <  1:1030840/149‑1 (MQ=255)
                                                                                 gAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGt                                                       >  1:304103/1‑149 (MQ=255)
                                                                                         cTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACg                                               >  2:162149/1‑149 (MQ=255)
                                                                                          ttCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGc                                              >  1:437090/1‑149 (MQ=255)
                                                                                             ggCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTATCCGCCCGGAGAGGGTTAGTAGTTACAGCGATAACTCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCAGCCAGTAATGCCTGTCGCTGACGCTGt                                           >  2:910826/1‑149 (MQ=255)
                                                                                              gCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTa                                          >  1:122929/1‑149 (MQ=255)
                                                                                                aGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTaaa                                        >  1:39948/1‑149 (MQ=255)
                                                                                                 gTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAt                                       <  2:954399/149‑1 (MQ=255)
                                                                                                  ttCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTGAGTAGTTCCAGCGATACCGCTGCAAGTCGATTTGCCCGCTTCCCGATGCGAGCACACCTTCTGCACGTAATGCCAGTCGCTGACGCTGTACAtc                                      >  1:708590/1‑149 (MQ=255)
                                                                                                   tCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCCCACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAACtcg                                     >  1:1063143/1‑148 (MQ=255)
                                                                                                      aGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCTCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCTTGTCGCTGACGCTGTAAAtccgg                                   >  2:47661/1‑148 (MQ=255)
                                                                                                      aGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCTCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCTTGTCGCTGACGCTGTAAAtccgg                                   <  1:47661/148‑1 (MQ=255)
                                                                                                            tACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCTAGTCGATTTGCCCGCTTCACGATACCATCACAACTGCTTCCAGTAATGCCGGTCGATGACGATGTAAATCCGGGGCGGtt                            >  2:1183755/1‑149 (MQ=255)
                                                                                                             aCGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTa                           >  2:727961/1‑149 (MQ=255)
                                                                                                             aCGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTa                           >  1:752632/1‑149 (MQ=255)
                                                                                                                 aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGa                       >  2:854592/1‑149 (MQ=255)
                                                                                                                 aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGa                       >  2:1161180/1‑149 (MQ=255)
                                                                                                                 aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGAAACCATCACACCATCTGCCACTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGa                       >  2:185928/1‑149 (MQ=255)
                                                                                                                 aaTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGATAGTAGTTCCAGCGATAACGCTTCAAGTCGATTTGCCCGCTTCCAGAAACCATCACACCTTCTGCCAGAAATGCCTGTCTCTGACGCTGTAAATCCGGTCCGGTTAGCGa                       >  1:677499/1‑149 (MQ=255)
                                                                                                                   tatataGCGTCGTTAACCCTCTCCGCCCGGAGAGTGTTAGTTGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGaaa                     <  2:323956/144‑1 (MQ=255)
                                                                                                                     ttAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAtt                   >  1:170715/1‑149 (MQ=255)
                                                                                                                       aaGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgt                 >  2:449029/1‑149 (MQ=255)
                                                                                                                         gaactactAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTTAATCCGGTCCGGTGAGCGAAATTTTGc               >  2:361723/8‑149 (MQ=255)
                                                                                                                         gCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGGCCGGTTAGAGAAATTGTGc               >  1:1104818/1‑149 (MQ=255)
                                                                                                                          cGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGcc              <  2:378329/149‑1 (MQ=255)
                                                                                                                          cGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGcc              <  1:162149/149‑1 (MQ=255)
                                                                                                                           gACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCGGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATACCCCCTTCGGCCAGTAAGGCCTGTCGCTGACGCGTTAAATCCGGGCGGGTAACCGAAAGTGAGcag             >  2:988894/1‑147 (MQ=255)
                                                                                                                                      tctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAac  <  1:384725/149‑1 (MQ=255)
                                                                                                                                      tctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAac  <  2:90094/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
GCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAAC  >  NC_000913/475831‑476113

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