Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,549,932 T→A 5.0% intergenic (+222/+21) fimH → / ← gntP minor component of type 1 fimbriae/fructuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,9320TA 5.0% 554.3 / 13.2 159intergenic (+222/+21)fimH/gntPminor component of type 1 fimbriae/fructuronate transporter
Reads supporting (aligned to +/- strand):  ref base T (79/72);  new base A (4/4);  total (83/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 160 total aligned reads displayed.

TAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGA  >  NC_000913/4549784‑4550077
                                                                                                                                                    |                                                                                                                                                 
tAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAt                                                                                                                                                   <  4:485351/149‑1 (MQ=255)
  ggCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTg                                                                                                                                                 >  3:261774/1‑149 (MQ=255)
     cAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCa                                                                                                                                              >  3:25364/1‑149 (MQ=255)
              cgcgTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAGCCCACTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCAGAGCCTGCCAATGGCAGGctcagac                                                                                                                                     <  3:224991/149‑5 (MQ=39)
              cgcgTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                                                                     <  4:378766/149‑1 (MQ=255)
              cgcgTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCAGAGCCTGCCAATGGCAGGctcagac                                                                                                                                     <  3:175304/149‑5 (MQ=255)
                    aGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCtt                                                                                                                               <  4:765885/149‑1 (MQ=255)
                     gCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCttt                                                                                                                              <  2:585602/149‑1 (MQ=255)
                     gCGTCGCATCTGGCATAAACAAAGCGCAATTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTATCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCCTTGGCAGGCTCTGGTGTCCctt                                                                                                                              <  2:709309/149‑1 (MQ=255)
                      cGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTTTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCtttt                                                                                                                             <  2:275268/149‑1 (MQ=255)
                            aTCTGGCATAACCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCCGAGGGGCGAGGGGACTGTCCGGGCACCTTTTTATTCTTTGTCACAAGTCTGAGGCTGCAATTGGCAAGCTCTGCTGTCGTTTTATGCTa                                                                                                                       >  4:123445/1‑149 (MQ=25)
                             tCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTAc                                                                                                                      >  3:524009/1‑149 (MQ=255)
                              cTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTAcc                                                                                                                     <  3:206147/149‑1 (MQ=255)
                                  cATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGc                                                                                                                 <  3:719587/149‑1 (MQ=255)
                                   aTAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCt                                                                                                                >  4:654696/1‑149 (MQ=255)
                                    tAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCta                                                                                                               >  1:603501/1‑149 (MQ=255)
                                      aaCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGAAGGCTCTGGTGTCCTTTTACGCTACCATGCtaat                                                                                                             >  3:784351/1‑149 (MQ=255)
                                       aCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaata                                                                                                            >  2:310803/1‑149 (MQ=255)
                                       aCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGACCGGTCGGGGACAATTTTTGGCCTTTTCCTTAGTCTTGGGCTGGCATTGGCAAGCTCTGGTGTCCTTTTACGCCACCATTCtaataa                                                                                                           >  4:57755/1‑149 (MQ=19)
                                        cAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaataa                                                                                                           >  2:876845/1‑149 (MQ=255)
                                            gcgcACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAg                                                                                                       >  4:415620/1‑149 (MQ=255)
                                                 cTTTGCTGGTCTGTTCCCCTCACCCCAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaa                                                                                                  <  3:637060/149‑1 (MQ=255)
                                                  tttGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGca                                                                                                     >  4:626132/1‑145 (MQ=255)
                                                      cTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATc                                                                                             >  3:736135/1‑149 (MQ=255)
                                                       tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCa                                                                                            >  4:248998/1‑149 (MQ=255)
                                                           cTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGccc                                                                                        <  1:331410/149‑1 (MQ=255)
                                                            tGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTTCTTTTACGCTACCATGCTAATAACCAGCACAATAATCAGCCCa                                                                                       >  2:38062/1‑149 (MQ=255)
                                                            tGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCa                                                                                       >  4:366501/1‑149 (MQ=255)
                                                                 cccTCACCCTAACCCTCTCCCCGGAGGTGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCTGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCAc                                                                                  >  2:115406/1‑149 (MQ=255)
                                                                 cccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTgg                                                                                                                                                >  1:277028/1‑87 (MQ=255)
                                                                 cccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCAc                                                                                  >  2:658473/1‑149 (MQ=255)
                                                                  ccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACg                                                                                 <  2:603501/149‑1 (MQ=255)
                                                                  ccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACg                                                                                 >  2:382127/1‑149 (MQ=255)
                                                                  ccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACg                                                                                 >  2:791864/1‑149 (MQ=255)
                                                                   cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACgg                                                                                >  2:175452/1‑149 (MQ=255)
                                                                    tCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCTGCCCAACCACGGa                                                                               >  4:558538/1‑149 (MQ=255)
                                                                    tCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGa                                                                               <  3:290260/149‑1 (MQ=255)
                                                                    tCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGa                                                                               <  4:338129/149‑1 (MQ=255)
                                                                    gcACCCTAACCCTCTCCCCGTAGGAGCGAGGGGACTGTCCGGGCACATTTTTCGACTTTGTCATCAGTCTGCGCCTGCCATTGGGAGGCTCTGGTGTCCTTTTACGCTACCATGCTAAGAATCAGCACAATAATCAGCGCAACCAcaca                                                                               <  4:634137/148‑4 (MQ=38)
                                                                      aCCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGt                                                                             <  2:85859/149‑1 (MQ=255)
                                                                        ccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                                                                           >  2:509280/1‑149 (MQ=255)
                                                                        ccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                                                                           >  2:484020/1‑149 (MQ=255)
                                                                        ccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                                                                           >  1:301777/1‑149 (MQ=255)
                                                                        ccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                                                                           >  3:392017/1‑149 (MQ=255)
                                                                        ccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                                                                           >  2:822708/1‑149 (MQ=255)
                                                                         cTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGa                                                                          <  4:23184/149‑1 (MQ=255)
                                                                          tAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGa                                                                               <  2:589526/143‑1 (MQ=255)
                                                                          tAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAc                                                                         <  1:494867/149‑1 (MQ=255)
                                                                               ctctCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCt                                                                    >  2:812104/1‑149 (MQ=255)
                                                                                tctcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGc                                                                                                                 <  4:44596/103‑1 (MQ=255)
                                                                                    cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGcc                                                                                         >  3:476615/1‑123 (MQ=255)
                                                                                    cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGcc                                                                                         <  4:476615/123‑1 (MQ=255)
                                                                                    cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagc                                                               >  1:369473/1‑149 (MQ=255)
                                                                                     ccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagca                                                              >  4:87612/1‑149 (MQ=255)
                                                                                       ggAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATTAGCCCAACCACGGAGTTGACCAGCTCCAGCAGa                                                            <  3:558538/149‑1 (MQ=255)
                                                                                          ggggCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAccc                                                         <  3:346657/149‑1 (MQ=255)
                                                                                             tcGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTGTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAg                                                      <  3:79255/148‑1 (MQ=255)
                                                                                               gAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGt                                                    <  4:806662/149‑1 (MQ=255)
                                                                                                 ggggACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGttt                                                  >  4:714512/1‑149 (MQ=255)
                                                                                                  gggACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGtttt                                                 >  1:18957/1‑149 (MQ=255)
                                                                                                   ggACTGTCCGGGCACATTTTTATACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTc                                                <  1:418028/149‑1 (MQ=255)
                                                                                                     aCTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCaa                                              <  3:87612/149‑1 (MQ=255)
                                                                                                     aCTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCaa                                              <  2:210787/149‑1 (MQ=255)
                                                                                                      cTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAAc                                             <  3:603649/149‑1 (MQ=255)
                                                                                                        gTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACgt                                           <  3:57755/149‑1 (MQ=255)
                                                                                                         tCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACgtg                                          >  1:716629/1‑149 (MQ=255)
                                                                                                         tCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACgtg                                          >  1:725621/1‑149 (MQ=255)
                                                                                                               cacaTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCtttt                                    <  2:735119/149‑1 (MQ=255)
                                                                                                                    ttttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGa                               >  1:322702/1‑149 (MQ=255)
                                                                                                                     tttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGAc                              >  1:58885/1‑149 (MQ=255)
                                                                                                                      ttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAg                                                      >  1:101979/1‑124 (MQ=255)
                                                                                                                      ttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAg                                                      <  2:101979/124‑1 (MQ=255)
                                                                                                                       tAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAg                            <  1:641028/149‑1 (MQ=255)
                                                                                                                         gACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACa                             >  3:149168/1‑146 (MQ=255)
                                                                                                                         gACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACa                             <  4:149168/146‑1 (MQ=255)
                                                                                                                         gACTTTGTCAGCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTa                                   >  1:816950/1‑140 (MQ=255)
                                                                                                                            tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                       >  2:167209/1‑149 (MQ=255)
                                                                                                                            tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCAAAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                       <  1:577751/149‑1 (MQ=255)
                                                                                                                             ttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaaa                      <  3:123445/149‑1 (MQ=255)
                                                                                                                             ttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaaa                      <  2:725621/149‑1 (MQ=255)
                                                                                                                             ttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaaa                      <  3:398529/149‑1 (MQ=255)
                                                                                                                             ttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaaa                      <  2:218752/149‑1 (MQ=255)
                                                                                                                                  atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAc                 >  1:641317/1‑149 (MQ=255)
                                                                                                                                  atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAc                 <  1:157197/149‑1 (MQ=255)
                                                                                                                                  atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAc                 >  2:250706/1‑149 (MQ=255)
                                                                                                                                  atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAc                 >  1:374697/1‑149 (MQ=255)
                                                                                                                                  atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAc                 >  2:235599/1‑149 (MQ=255)
                                                                                                                                      accagAGCCTGCCAATGGCAGGCTCAGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTaa                  <  1:470882/140‑1 (MQ=255)
                                                                                                                                      accagAGCCTGCCAATGGCAGGCTCAGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTaa                  >  2:470882/5‑144 (MQ=255)
                                                                                                                                       tCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTg            >  3:872948/1‑149 (MQ=255)
                                                                                                                                        cagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTg            <  4:824403/146‑1 (MQ=255)
                                                                                                                                        cagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTCTGa           >  4:572444/3‑149 (MQ=255)
                                                                                                                                             ccTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaacaac      <  1:569115/149‑1 (MQ=255)
                                                                                                                                             ccTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaacaac      <  1:658473/149‑1 (MQ=255)
                                                                                                                                               tGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCa    >  1:795626/1‑149 (MQ=255)
                                                                                                                                               tGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCa    >  1:606333/1‑149 (MQ=255)
                                                                                                                                               tGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCa    >  1:751412/1‑149 (MQ=255)
                                                                                                                                                gCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGCTCAAAGTAACCTTTGAACAACCAg   >  4:269202/1‑149 (MQ=255)
                                                                                                                                                   aTTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAcccc                                                        >  3:358620/1‑93 (MQ=255)
                                                                                                                                                   aTTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAcccc                                                        <  4:358620/93‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGA  >  NC_000913/4549784‑4550077

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