Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 894,627 A→G 5.2% M282V (ATG→GTG)  potF → putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913894,6270AG 5.2% 240.5 / 2.7 77M282V (ATG→GTG) potFputrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base G (2/2);  total (35/42)
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.13e-01

CCTGGCAAACGGCGATATTTGCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGA  >  NC_000913/894479‑894775
                                                                                                                                                    |                                                                                                                                                    
ccTGGCAAACGGCGATATTTGCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCa                                                                                                                                                      <  6:179760/149‑1 (MQ=255)
ccTGGCAAACGGCGATATTTGCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCa                                                                                                                                                      <  2:272132/149‑1 (MQ=255)
 cTGGCAAACGGCGATATTTGCGTCGCTATCGGCTGGGCGGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCAt                                                                                                                                                     <  8:379918/149‑1 (MQ=255)
 cTGGCAAACGGCGATATTTGCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCAt                                                                                                                                                     >  1:169942/1‑149 (MQ=255)
    gCAAACGGCGATATTTGCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGcc                                                                                                                                                  >  1:118292/1‑149 (MQ=255)
       aaCGGCGATATTTGCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGc                                                                                                                                               <  7:371258/149‑1 (MQ=255)
         cGGCGATATTTGCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCAt                                                                                                                                                     >  6:194916/1‑141 (MQ=255)
         cGGCGATATTTGCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCAt                                                                                                                                                     <  5:194916/141‑1 (MQ=255)
                atTTGCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCa                                                                                                                                       >  3:178026/1‑148 (MQ=255)
                 attGCGGCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGGCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCaa                                                                                                                                      <  7:162594/147‑1 (MQ=255)
                    gCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCaaaaa                                                                                                                                   <  1:156896/148‑1 (MQ=255)
                      gTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACa                                                                                                                                 >  7:246470/1‑148 (MQ=255)
                       tCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCa                                                                                                                                                      >  8:58693/1‑126 (MQ=255)
                       tCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCa                                                                                                                                                      <  7:58693/126‑1 (MQ=255)
                       tCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACaaa                                                                                                                               <  2:169942/149‑1 (MQ=255)
                         gCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGa                                                                                                                             <  5:129508/149‑1 (MQ=255)
                         gCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCAGATGCCAAAAACAAAGa                                                                                                                             <  2:277371/149‑1 (MQ=255)
                             tCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCTCGAATGAATCGAAGAATTGAGTGAATGTCTCGATCTCGATTCCAAAATAAGGGGCGATGGAGTTCTTTGATGTATTAGCCATGCCTGCGGATGCCAAAAACAAAGACGaa                                                                                                                         <  6:291343/149‑1 (MQ=255)
                             tCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGa                                                                                                                          <  6:25385/148‑1 (MQ=255)
                              cGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGc                                                                                                                                               <  1:102479/126‑1 (MQ=255)
                              cGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGc                                                                                                                                               >  2:102479/1‑126 (MQ=255)
                               ggCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGc                                                                                                                       <  1:345724/149‑1 (MQ=255)
                                    ggCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTAt                                                                                                                   <  8:246470/148‑1 (MQ=255)
                                         gTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGTCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTc                                                                                                             <  1:234660/149‑1 (MQ=255)
                                         gTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTc                                                                                                             >  4:83651/1‑149 (MQ=255)
                                         gTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGGTCTTTGATGTATTCGACATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTc                                                                                                             >  3:305005/1‑149 (MQ=255)
                                          tGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTcc                                                                                                            >  1:149692/1‑149 (MQ=255)
                                                 tGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGAGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTATTTGATGTATTCGCCATGCCTGCCGATGCAAAAAACAAATACGAATCCTATCATTTCCTGAAtt                                                                                                      >  4:24649/1‑148 (MQ=255)
                                                   gCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTAcc                                                                                                   >  5:146130/1‑149 (MQ=255)
                                                       gCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACctgct                                                                                               >  8:250436/1‑149 (MQ=255)
                                                          tCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTgcg                                                                                            >  3:186249/1‑149 (MQ=255)
                                                               ccGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCGATCAGTTCCTGAATTACCTGCTGCGCCCg                                                                                        >  1:294611/1‑148 (MQ=255)
                                                                        ggaagCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGtagtag                                                                              <  4:281429/149‑1 (MQ=255)
                                                                         gaagCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAgc                                                                             <  4:132122/149‑1 (MQ=255)
                                                                              gaagaaTGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCata                                                                        <  7:250436/149‑1 (MQ=255)
                                                                              gaagaaTGGCGTGAATGTCTCGTTATCGATTCAAAAAGAAGTGGCGATGGCGTTCTTTGATGTATTCGCCATGCATGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTACTGAATTACCTGATGCGCCAGGATGTAGTAGCGCata                                                                        >  6:109142/1‑149 (MQ=255)
                                                                                agaaTGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCGTGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATAtt                                                                      >  3:223078/1‑149 (MQ=255)
                                                                                agaaTGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCGTGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATAtt                                                                      >  3:227864/1‑149 (MQ=255)
                                                                                  aaTGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTc                                                                    <  2:149692/149‑1 (MQ=255)
                                                                                      tcgggAATGTCTCGTTCTCGATTACAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCGTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGa                                                                 <  6:321712/144‑1 (MQ=255)
                                                                                      gCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGa                                                                 <  6:157423/148‑1 (MQ=255)
                                                                                        gTGAAGGTCTCTTTCTCTATTACAAAAGCAGGTGCGCCGCCTTTAATTGATTTATTCGCAATTCTTTCGTATTCGAAAATCATAGTCTATGCCTATAACTTTCTGATTTTCAAGATGCGCCCGGATGTAGTAGCGCATATTTCCGACCa                                                              <  2:78882/149‑1 (MQ=255)
                                                                                         tGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCAt                                                             <  4:178026/149‑1 (MQ=255)
                                                                                           aaTGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTAc                                                                                                    <  3:142827/108‑1 (MQ=255)
                                                                                           aaTGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTAc                                                                                                    >  4:142827/1‑108 (MQ=255)
                                                                                           aaTGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGGTCTTTGATGTATTCGCCGTGCCGGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCAtgt                                                           <  4:227864/149‑1 (MQ=255)
                                                                                            aTGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCGTGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCAtgt                                                           <  4:223078/148‑1 (MQ=255)
                                                                                            aTGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCAtgtg                                                          <  3:83651/149‑1 (MQ=255)
                                                                                             tGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCAGTCCTGCGGATGCCAAAAACAAAGCCGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGAcccttggt                                                         >  5:100649/1‑143 (MQ=255)
                                                                                               tctcGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGTAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGTATGTAGTAGCGCATATTTCCGACCATGTGTTc                                                       <  6:316977/149‑1 (MQ=255)
                                                                                               tctcGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGtt                                                        <  4:186249/148‑1 (MQ=255)
                                                                                                   gtgctcGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGACGCCTATCAGTTCCTGAATTACCTGCGGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTGCTa                                                     <  1:236975/144‑1 (MQ=255)
                                                                                                      ctcGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGcc                                                 >  4:205447/1‑148 (MQ=255)
                                                                                                       tcGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGtt                                                        >  4:262622/1‑140 (MQ=255)
                                                                                                       tcGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGtt                                                        <  3:262622/140‑1 (MQ=255)
                                                                                                         gATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATgtagta                                                                               >  4:58472/1‑115 (MQ=255)
                                                                                                         gATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATgtagta                                                                               <  3:58472/115‑1 (MQ=255)
                                                                                                         gATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGAcc                                                               >  7:213087/1‑131 (MQ=255)
                                                                                                         gATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGAcc                                                               <  8:213087/131‑1 (MQ=255)
                                                                                                         gATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACg                                             >  5:198459/1‑149 (MQ=255)
                                                                                                              cAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCaa                                         <  1:348036/148‑1 (MQ=255)
                                                                                                              cAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCaa                                         >  1:68979/1‑148 (MQ=255)
                                                                                                              cAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAAt                                        <  2:118292/149‑1 (MQ=255)
                                                                                                               aaaaGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGACGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAAt                                        >  8:189015/1‑148 (MQ=255)
                                                                                                               aaaaGAAGGGGCGATGGCGTTCTTTGATGAATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACAATGTGTTCTATGCCAACGCCAAt                                        >  4:103535/1‑148 (MQ=255)
                                                                                                                aaaGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGtt                                                        <  1:38256/131‑1 (MQ=255)
                                                                                                                aaaGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGtt                                                        >  2:38256/1‑131 (MQ=255)
                                                                                                                    aattggCTATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGTATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAagca                                  <  3:24649/145‑1 (MQ=255)
                                                                                                                    aagGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAagca                                  <  6:146130/149‑1 (MQ=255)
                                                                                                                       gggCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAATACCAAAGACGAAGCGGATCAGTTCCTGAATTACCTGCAGAGCCCCGAGGTAGTAGAGCATATTTCCGAACATGTGTTCTAGGCACAAGCCAATAACTCAGcc                               >  6:189767/1‑149 (MQ=255)
                                                                                                                            aTGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATTCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTATATGCCAACGCCAATAAAGCGGCCAcgcc                          >  4:187645/1‑149 (MQ=255)
                                                                                                                            aTGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCAcgcc                          >  3:266848/1‑149 (MQ=255)
                                                                                                                               tctttCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATTACAACGCCAATAAAGCAGCCACGCCGCt                       <  5:109142/146‑1 (MQ=255)
                                                                                                                                 gTTCTTTGATGTATTCGCCATGCCTGATGATGCCAAAAACAAAGACGAAGCATATCAGTTACTGACTTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCAGGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTgg                     <  6:100649/149‑1 (MQ=255)
                                                                                                                                     tttGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATgtgt                                                         <  1:81044/109‑1 (MQ=255)
                                                                                                                                     tttGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATgtgt                                                         >  2:81044/1‑109 (MQ=255)
                                                                                                                                             aTTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCAAACTCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAg         >  6:134680/1‑149 (MQ=255)
                                                                                                                                              ttCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCg             >  5:182279/1‑144 (MQ=255)
                                                                                                                                              ttCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCg             <  6:182279/144‑1 (MQ=255)
                                                                                                                                              ttCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGt        >  2:321378/1‑149 (MQ=255)
                                                                                                                                                 gCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTcc      <  6:198459/148‑1 (MQ=255)
                                                                                                                                                    aTGCCTGCGGATGCCAAAAACAAAGACTAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTga  <  3:205447/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTGGCAAACGGCGATATTTGCGTCGCTATCGGCTGGGCAGGTGATGTCTGGCAGGCGTCAAACCGCGCGAAGGAAGCGAAGAATGGCGTGAATGTCTCGTTCTCGATTCCAAAAGAAGGGGCGATGGCGTTCTTTGATGTATTCGCCATGCCTGCGGATGCCAAAAACAAAGACGAAGCCTATCAGTTCCTGAATTACCTGCTGCGCCCGGATGTAGTAGCGCATATTTCCGACCATGTGTTCTATGCCAACGCCAATAAAGCAGCCACGCCGCTGGTGAGTGCGGAAGTCCGTGA  >  NC_000913/894479‑894775

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