Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,291,373 G→A L11L (TTG→TTA yraK → putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,291,3730GA97.6% 305.6 / ‑1.7 84L11L (TTG→TTAyraKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (32/50);  total (32/52)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACCTCTGCAA  >  NC_000913/3291226‑3291518
                                                                                                                                                   |                                                                                                                                                 
gTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTttattg                                                                                                                                                <  1:18355/151‑1 (MQ=255)
gTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTttatt                                                                                                                                                 <  1:350574/150‑1 (MQ=255)
gTACGCTGGGGAAAGCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATATCAGGTCTTttatt                                                                                                                                                 <  1:1091656/150‑1 (MQ=255)
  aCGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGa                                                                                                                                               <  1:1126078/150‑1 (MQ=255)
   cGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGata                                                                                                                                             <  2:1106244/151‑1 (MQ=255)
       ggggAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCca                                                                                                                                         <  1:13700/151‑1 (MQ=255)
        gggAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCcac                                                                                                                                        <  2:19704/151‑1 (MQ=255)
         ggAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCcaca                                                                                                                                       <  2:1208800/151‑1 (MQ=255)
          gAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACAt                                                                                                                                      <  2:954151/151‑1 (MQ=255)
           aaaTCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATc                                                                                                                                     <  1:140690/151‑1 (MQ=255)
             aTCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCt                                                                                                                                   <  2:549671/151‑1 (MQ=255)
               cAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTgc                                                                                                                                 <  2:1112670/151‑1 (MQ=255)
                aGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTgcg                                                                                                                                >  2:188167/1‑151 (MQ=255)
                  gCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCt                                                                                                                              <  1:1073941/151‑1 (MQ=255)
                      aTCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATg                                                                                                                          >  1:475439/1‑151 (MQ=255)
                               gTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCag                                                                                                                 <  1:1034936/151‑1 (MQ=255)
                               gTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCag                                                                                                                 <  1:781573/151‑1 (MQ=255)
                               gTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCa                                                                                                                  <  2:491534/150‑1 (MQ=255)
                               gTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCa                                                                                                                  <  2:280579/150‑1 (MQ=255)
                                 cAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCagga                                                                                                               <  2:244334/151‑1 (MQ=255)
                                   gTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCaggagg                                                                                                             <  1:862892/151‑1 (MQ=255)
                                          ttCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGtgtg                                                                                                        <  1:727292/149‑1 (MQ=255)
                                                      gATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACt                                                                                          <  1:697617/151‑1 (MQ=255)
                                                      gATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACt                                                                                          <  2:25227/151‑1 (MQ=255)
                                                        ttCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACt                                                                                          >  1:1064896/1‑149 (MQ=255)
                                                            gAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGt                                                                                    <  1:495312/151‑1 (MQ=255)
                                                             aaCCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGt                                                                                    <  1:552226/150‑1 (MQ=255)
                                                             aaCCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTg                                                                                   <  1:854822/151‑1 (MQ=255)
                                                              aCCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGc                                                                                  <  1:443324/151‑1 (MQ=255)
                                                              aCCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGc                                                                                  >  2:585407/1‑151 (MQ=255)
                                                               ccTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCg                                                                                 <  2:974718/151‑1 (MQ=255)
                                                                cTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCg                                                                                 >  2:319514/1‑150 (MQ=255)
                                                                 tACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGa                                                                                >  2:26029/1‑150 (MQ=255)
                                                                  aCAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCTCATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGAc                                                                               <  1:401052/150‑1 (MQ=255)
                                                                  aCAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGAc                                                                               <  2:1017374/150‑1 (MQ=255)
                                                                  aCAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGAc                                                                               <  2:972255/150‑1 (MQ=255)
                                                                    aaCCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACaaa                                                                            <  2:968930/151‑1 (MQ=255)
                                                                    aaCCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACaa                                                                             <  2:514841/150‑1 (MQ=255)
                                                                    aaCCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACaa                                                                             >  2:1022248/1‑150 (MQ=255)
                                                                     aCCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAAt                                                                           >  2:1251251/1‑151 (MQ=255)
                                                                           gTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTAca                                                                       <  2:402235/149‑1 (MQ=255)
                                                                           gTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTAca                                                                       <  2:568934/149‑1 (MQ=255)
                                                                           gTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTAca                                                                       <  1:523982/149‑1 (MQ=255)
                                                                            tATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGTGACAAATTAcagc                                                                     >  1:1162451/1‑150 (MQ=255)
                                                                             aTTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTAcagca                                                                    <  2:617853/150‑1 (MQ=255)
                                                                                aCAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAgtg                                                                 >  2:261900/1‑150 (MQ=255)
                                                                                    ggCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGtt                                                            >  2:24690/1‑151 (MQ=255)
                                                                                       acaGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTATAGCAGTGTGGTtgat                                                         >  1:1030761/1‑151 (MQ=255)
                                                                                        caGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTtgatg                                                        >  1:839163/1‑151 (MQ=255)
                                                                                          gcgcAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTgatgat                                                      <  1:98215/151‑1 (MQ=255)
                                                                                            gcAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGt                                                    <  2:845081/151‑1 (MQ=255)
                                                                                            gcAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGt                                                    >  1:480620/1‑151 (MQ=255)
                                                                                              aGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTa                                                  <  2:30918/151‑1 (MQ=255)
                                                                                               gTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTAc                                                 <  1:109102/151‑1 (MQ=255)
                                                                                                tGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACg                                                >  2:525813/1‑151 (MQ=255)
                                                                                                 gCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACgg                                               <  2:105821/151‑1 (MQ=255)
                                                                                                 gCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACg                                                <  2:132075/150‑1 (MQ=255)
                                                                                                  ccATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGt                                              >  2:808209/1‑151 (MQ=255)
                                                                                                     ttAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGa                                            >  2:1099244/1‑150 (MQ=255)
                                                                                                       aaGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAAc                                          <  1:857588/150‑1 (MQ=255)
                                                                                                             aaaaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagacag                                   <  2:1196126/151‑1 (MQ=255)
                                                                                                             aaaaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagacag                                   <  2:845898/151‑1 (MQ=255)
                                                                                                             aaaaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagaca                                    <  2:769138/150‑1 (MQ=255)
                                                                                                             aaaaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagac                                     >  2:258868/1‑149 (MQ=255)
                                                                                                              aaaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagacaga                                  <  2:850403/151‑1 (MQ=255)
                                                                                                              aaaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagac                                     >  1:261585/1‑148 (MQ=255)
                                                                                                               aaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagac                                     >  1:408797/1‑147 (MQ=255)
                                                                                                               aaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGAt                                 <  1:740246/151‑1 (MQ=255)
                                                                                                               aaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGAt                                 <  2:1040260/151‑1 (MQ=255)
                                                                                                               aaaaaTGAAAAGAGCGCCTCTTATAACAGGACTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGCAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagac                                     >  1:618761/1‑147 (MQ=255)
                                                                                                                aaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGAt                                 >  1:219645/1‑150 (MQ=255)
                                                                                                                   aTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACg                             <  2:1005014/151‑1 (MQ=255)
                                                                                                                     gAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTg                           <  1:20002/151‑1 (MQ=255)
                                                                                                                       aaaGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGAc                         >  1:503979/1‑151 (MQ=255)
                                                                                                                        aaGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGt                            >  2:794866/1‑147 (MQ=255)
                                                                                                                        aaGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGAc                         >  1:191754/1‑150 (MQ=255)
                                                                                                                         agagCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAg                       <  1:1065161/151‑1 (MQ=255)
                                                                                                                          gagCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAgg                      <  2:414876/151‑1 (MQ=255)
                                                                                                                           agCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAgg                      >  1:229582/1‑150 (MQ=255)
                                                                                                                              gcCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACgg                  >  1:1005848/1‑151 (MQ=255)
                                                                                                                                  ctTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAg               >  1:328251/1‑150 (MQ=255)
                                                                                                                                  ctTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGCGGTTGATGATGTTACGGTGAACCAGACAGATAATGTGACAGGACGGGa                >  2:943865/1‑149 (MQ=255)
                                                                                                                                    tataACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGttt            >  1:1113752/1‑151 (MQ=255)
                                                                                                                                     ataACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGttt            <  2:1201947/150‑1 (MQ=255)
                                                                                                                                      taACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGGGAACCAGACAGATAACGTGACAGGACGGGAGTTTAc          >  1:900154/1‑151 (MQ=255)
                                                                                                                                             tCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACCTCTGc    >  1:517139/1‑150 (MQ=255)
                                                                                                                                                tttATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACCTCTGCaa  <  2:146197/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACCTCTGCAA  >  NC_000913/3291226‑3291518

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