Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 475,971 C→G 5.0% intergenic (+20/+11) 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,9710CG 5.0% 210.4 / 4.2 80intergenic (+20/+11)ybaY/ybaZouter membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL
Reads supporting (aligned to +/- strand):  ref base C (34/42);  new base G (2/2);  total (36/44)
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.82e-01

ATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACC  >  NC_000913/475825‑476114
                                                                                                                                                  |                                                                                                                                               
aTTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCACGCCGccc                                                                                                                                               >  2:139024/1‑149 (MQ=255)
        cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                       <  2:623196/149‑1 (MQ=255)
        cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                                                                       <  2:265930/149‑1 (MQ=255)
         ggTACAGCACACCGCCGTGCCGGTTCAGGCCCGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCACCTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGt                                                                                                                                      <  1:177703/149‑1 (MQ=255)
         ggTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGt                                                                                                                                      <  1:832476/149‑1 (MQ=255)
         ggTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGt                                                                                                                                      >  2:899526/1‑149 (MQ=255)
          gTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGtt                                                                                                                                     <  2:588856/149‑1 (MQ=255)
             cAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAGGGTtaac                                                                                                                                  <  1:814979/149‑3 (MQ=255)
              aGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagtt                                                                                                                                 >  1:297346/1‑148 (MQ=255)
                cAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTagtagt                                                                                                                               >  1:27929/1‑149 (MQ=255)
                 aaaCCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGGGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGtt                                                                                                                              <  2:284501/149‑1 (MQ=255)
                  aaCCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTc                                                                                                                             <  2:381565/149‑1 (MQ=255)
                   aCCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTcc                                                                                                                            <  1:886489/149‑1 (MQ=255)
                      gccgTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGc                                                                                                                         <  2:868085/149‑1 (MQ=255)
                       ccgTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCg                                                                                                                        >  1:855797/1‑149 (MQ=255)
                         gttgcGGTTCAGGCCATCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCATCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAt                                                                                                                      <  2:55241/145‑1 (MQ=255)
                         gTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTTCCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCATCCCGGAGAGGGTTAGTAGTTCCAGCGAt                                                                                                                      <  2:615575/149‑1 (MQ=255)
                         gTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACATCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAt                                                                                                                      <  1:518473/149‑1 (MQ=255)
                          tGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATa                                                                                                                     >  2:160191/1‑149 (MQ=255)
                           gCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATaa                                                                                                                    <  1:173229/149‑1 (MQ=255)
                              ggTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATAATAAGCGACGTTAACCCTCTCCGCCGGGAGAGGGGTAGTTGTTCCAGCGATAACGc                                                                                                                 >  2:418384/1‑149 (MQ=255)
                                ttCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGCATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTg                                                                                                               >  1:768089/1‑149 (MQ=255)
                                ttCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAAAATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCAAGCGATAACGCTg                                                                                                               >  1:433421/1‑149 (MQ=255)
                                 tCAGGCCAGCGGTGGCGCAACGACTCCCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGc                                                                                                              <  2:720398/149‑1 (MQ=255)
                                   aGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCaa                                                                                                            >  2:72364/1‑149 (MQ=255)
                                   aGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTGCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGTTATCGCTGGaa                                                                                                            <  1:343683/149‑1 (MQ=255)
                                         gCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGAt                                                                                                      <  2:452558/149‑1 (MQ=255)
                                                   aaCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCtt                                                                                            >  1:548339/1‑149 (MQ=255)
                                                    aCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCAGTCTTCGGCAGATCCAGCTCCAACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTc                                                                                           <  1:562725/149‑1 (MQ=255)
                                                        cTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTAGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGa                                                                                       <  2:535009/149‑1 (MQ=255)
                                                           ctgTACCTTCGACTTACCCCACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTCAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGGTCCAGCGCTAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAc                                                                                    <  2:116702/147‑1 (MQ=255)
                                                           ccGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAGAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAc                                                                                    >  1:606890/1‑149 (MQ=255)
                                                           ccGTACCTTCGACTTCACAAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTGCTACGCAATATTAAGCGACGTTAACCCTCGCCGCCCAGAGAGGGTAAGTATTTCCAGCGATAACGCTGCACGTCGCTTAGTCCGCTTCCCGATaa                                                                                    >  1:585315/1‑148 (MQ=255)
                                                              tACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACTAt                                                                                 <  1:72364/149‑1 (MQ=255)
                                                              tACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCAt                                                                                 <  2:514937/149‑1 (MQ=255)
                                                                 cTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcac                                                                              >  1:837627/1‑149 (MQ=255)
                                                                  ttCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcaca                                                                             <  2:27929/149‑1 (MQ=255)
                                                                      aCTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCCAGTCGATTTGCCCGCTTCCCGCTACCATCACACCtt                                                                         >  1:7238/1‑149 (MQ=255)
                                                                      aCTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTACAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCtt                                                                         >  1:11055/1‑149 (MQ=255)
                                                                       cTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTc                                                                        <  1:686281/149‑1 (MQ=255)
                                                                          cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAAAATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGGAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                                                                     >  2:229416/1‑149 (MQ=255)
                                                                          cACCAACTCAGGTGAATCCGGCTTCGGCAGTTCCAGCTCCTACGCAATCTTAAGCGACGGTAACCCTCTCCTCCCGTAGAGGGTTAGTAGTTCCCGCGATCACGCTGCAAGTCGATTTGCCCGCTTCCCGCTACCATCACACCTTCTGc                                                                     <  2:588618/149‑1 (MQ=255)
                                                                                   aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTACCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGc                                                            <  1:719857/149‑1 (MQ=255)
                                                                                   aGGTGAATCCGTCTTCGGCAGTTCCAGCTCATACGCAATATTAAGCGACGTTAACCCTCTCCTCCCGGAGAGGGTTAGTAGTTCCATCGATAACGCTCCAAGTCGCTTTGCCCGCTTCCCGATACCATGACAACTTCTGCCAGGAATGc                                                            >  1:116948/1‑149 (MQ=255)
                                                                                    ggTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGcc                                                           <  1:178682/149‑1 (MQ=255)
                                                                                          tCCGTATTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAACTCGATTTGCCCGCTTCCCGCTACCATCACACCTTCTGCCAGTATTGCATGTCGc                                                     >  2:631340/1‑149 (MQ=255)
                                                                                           ccGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCt                                                    <  1:342261/149‑1 (MQ=255)
                                                                                           ccGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGATAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTTCCAGTAATGCCTGTCGCt                                                    >  1:16474/1‑149 (MQ=255)
                                                                                            cGTCTGCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTg                                                   <  1:914096/149‑1 (MQ=255)
                                                                                             gTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATTCCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                  >  1:333417/1‑149 (MQ=255)
                                                                                               cTTCGGCAGTTCCAGCTCCTACGCAATTTTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGTTAACGCTGCAAGTCGATTTGCCCGCGTCCCGATACGATCGCACCTTCTGCCAGTAATGCCTGTCGCTGACg                                                >  1:332441/1‑149 (MQ=255)
                                                                                                ttCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTATTAGTTCCCGCGATAACGCTGCAAGTCGTTTTGCCCGCTTCACGATCCCAACACACCTTCTGCCAGTATTGCTTGTCGATGACGc                                               >  2:921101/1‑149 (MQ=255)
                                                                                                ttCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCCCCCCTTCTGCCAGTAATGCCTGTCGCTGACGc                                               >  2:148135/1‑149 (MQ=255)
                                                                                                ttCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGc                                               >  1:835328/1‑149 (MQ=255)
                                                                                                 tCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                  <  2:835415/145‑1 (MQ=255)
                                                                                                 tCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                  >  1:835415/1‑145 (MQ=255)
                                                                                                  cGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCTCCCGGAGAGGGTTAGTAGGTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCATTAATGCCTGTCGCTGACGCTg                                             <  1:999102/149‑1 (MQ=255)
                                                                                                  cGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTg                                             <  1:160191/149‑1 (MQ=255)
                                                                                                  cGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGAAACTCGATTTGCCCGCTTCCCGATACCATCACACCTTATGCCAGTAATGCCTGTCGCTGACGCTg                                             >  1:938362/1‑149 (MQ=255)
                                                                                                   ggCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTTACCCTATCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGt                                            <  2:824591/149‑1 (MQ=255)
                                                                                                   ggCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGt                                            <  1:899526/149‑1 (MQ=255)
                                                                                                    gCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTa                                           >  2:469873/1‑149 (MQ=255)
                                                                                                          ccAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATCCCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccg                                     >  1:138660/1‑149 (MQ=255)
                                                                                                             gCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtc                                  <  2:736950/149‑1 (MQ=255)
                                                                                                             gCTCCTACGCAATATAAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAGAACGCTGCAAGTCGATTTGCCCGATTCCCGATACCCACCAACCTCATGCCCGTATTGCCTGTCCCTCAAGCTTTAAAtccgggc                                  >  1:321013/1‑147 (MQ=255)
                                                                                                               tCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtccg                                <  1:732414/149‑1 (MQ=255)
                                                                                                               tCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtccg                                <  1:634122/149‑1 (MQ=255)
                                                                                                                ccTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGGTAGTAGATCCAGCGATAACGCTGACAGTCGATTTGCCCGCTTCCCGATACCATCACAACTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtccgg                               >  1:569652/1‑149 (MQ=255)
                                                                                                                 cTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATccggtccggt                              <  2:855797/149‑1 (MQ=255)
                                                                                                                          atTAAGCGCCGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTGCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAt                     <  2:753884/149‑1 (MQ=255)
                                                                                                                           ttAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAAtt                    <  1:22781/149‑1 (MQ=255)
                                                                                                                              aGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgtg                 <  2:835328/149‑1 (MQ=255)
                                                                                                                              aGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgtg                 >  2:886725/1‑149 (MQ=255)
                                                                                                                               gCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAAAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGAAGCTGTAAATCCGGTCCGGTTTGCTCAATTGTGc                >  1:565300/1‑149 (MQ=255)
                                                                                                                                cGACGTTAACCCTCTCCGCCCGGAGAGGGTGAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCGCGCTACCCTCTCCGCGTCTGCGAGTAATGCCTGTCGCTGAGGGTGTATATCCGGTCCGGTTTGCTCAATTGTGcc               >  2:688887/1‑149 (MQ=255)
                                                                                                                                  ctactAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATCCCATCACACCTTCTGCCTGTAATGCCTGTCGCTGCCGCTGTAACTCCGGTCGGGTTAGCGAAATTGTGCCGt             >  1:195658/5‑149 (MQ=255)
                                                                                                                                     ttAACCCTCTCCGCCAGGAGAGGGTTAGTAGGTCGAGCGATAACGCGGCAAGTCGAGTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGc          <  1:676939/149‑1 (MQ=255)
                                                                                                                                        aCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAt       >  2:986342/1‑149 (MQ=255)
                                                                                                                                         cccTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCGAGTAATGCCTGTCGCTGACGGTGTAAATCCGGTCCGGGTAGCGATATTGTGCCGTGGCGAtt      >  2:112385/1‑149 (MQ=255)
                                                                                                                                           ctctCCGCCCGGAGAGGGTTATTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTaa    <  2:768089/149‑1 (MQ=255)
                                                                                                                                             ctcCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCGGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAacc  <  2:140436/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACC  >  NC_000913/475825‑476114

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: