Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,549,932 T→A 6.5% 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 6.5% 424.7 / 16.2 123intergenic (+222/+21)fimH/gntPminor component of type 1 fimbriae/fructuronate transporter
Reads supporting (aligned to +/- strand):  ref base T (60/55);  new base A (5/3);  total (65/58)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 126 total aligned reads displayed.

TAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAA  >  NC_000913/4549784‑4550079
                                                                                                                                                    |                                                                                                                                                   
tAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAt                                                                                                                                                     >  3:418710/1‑149 (MQ=255)
    ccAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGc                                                                                                                                                 >  1:424320/1‑149 (MQ=255)
          aaGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctc                                                                                                                                           <  3:642431/149‑1 (MQ=255)
            gACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTg                                                                                                                                         >  4:87531/1‑149 (MQ=255)
            gACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTg                                                                                                                                         <  1:683526/149‑1 (MQ=255)
                 gTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTc                                                                                                                                    <  2:772009/149‑1 (MQ=255)
                 gTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTc                                                                                                                                    >  1:535814/1‑149 (MQ=255)
                  tCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTcc                                                                                                                                   >  1:736237/1‑149 (MQ=255)
                   cAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCt                                                                                                                                  <  3:311472/149‑1 (MQ=255)
                       gTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTa                                                                                                                              <  3:598234/149‑1 (MQ=255)
                       gTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTa                                                                                                                              >  2:528749/1‑149 (MQ=255)
                             tCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTGCCCCTCACCCTCACCGTCTCGCCGTAGGGGCGAGGGGACTGTCCGGGCACATTGTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTAc                                                                                                                        <  2:608612/149‑1 (MQ=39)
                              cTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctca                                                                                                                                          <  1:730917/130‑2 (MQ=255)
                               tGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCa                                                                                                                      >  4:556328/1‑149 (MQ=255)
                                   aTAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCt                                                                                                                  >  1:122560/1‑149 (MQ=255)
                                   aTAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCt                                                                                                                  >  2:169682/1‑149 (MQ=255)
                                      aaCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTgg                                                                                                                                        <  2:672936/124‑1 (MQ=255)
                                      aaCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTgg                                                                                                                                        >  1:672936/1‑124 (MQ=255)
                                      aaCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaat                                                                                                               >  4:301751/1‑149 (MQ=255)
                                      aaCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaat                                                                                                               >  4:389783/1‑149 (MQ=255)
                                        cAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaataa                                                                                                             >  3:603955/1‑149 (MQ=255)
                                        cAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaataa                                                                                                             >  4:713249/1‑149 (MQ=255)
                                        cAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaataa                                                                                                             >  2:616281/1‑149 (MQ=255)
                                         aaaGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTaataat                                                                                                            <  1:573742/149‑1 (MQ=255)
                                          aaGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTc                                                                                                                                    <  3:601112/124‑1 (MQ=255)
                                          aaGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTc                                                                                                                                    >  4:601112/1‑124 (MQ=255)
                                          aaGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATc                                                                                                           <  3:135858/149‑1 (MQ=255)
                                           aGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCa                                                                                                          <  4:326555/149‑1 (MQ=255)
                                            gcgcACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                                                                       >  1:526003/1‑119 (MQ=255)
                                            gcgcACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                                                                       <  2:526003/119‑1 (MQ=255)
                                             cgcACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGc                                                                                                        <  1:528749/149‑1 (MQ=255)
                                                 cTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTg                                                                                                                                         >  2:729565/1‑112 (MQ=255)
                                                 cTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTg                                                                                                                                         <  1:729565/112‑1 (MQ=255)
                                                 cTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaa                                                                                                    <  1:280520/149‑1 (MQ=255)
                                                  tttGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaat                                                                                                   >  3:165707/1‑149 (MQ=255)
                                                   ttGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaata                                                                                                  <  3:630243/149‑1 (MQ=255)
                                                       tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAtt                                                                                                                                                    >  1:88555/1‑95 (MQ=255)
                                                       tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAtt                                                                                                                                                    <  2:88555/95‑1 (MQ=255)
                                                        ggTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAg                                                                                             <  4:165707/149‑1 (MQ=255)
                                                        ggTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAg                                                                                             <  2:535814/149‑1 (MQ=255)
                                                          tCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGcc                                                                                           <  1:137367/149‑1 (MQ=255)
                                                                   cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCta                                                                                                                 >  2:285163/1‑118 (MQ=255)
                                                                   cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCta                                                                                                                 <  1:285163/118‑1 (MQ=255)
                                                                   cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCta                                                                                                                 <  1:431545/118‑1 (MQ=255)
                                                                   cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCta                                                                                                                 >  2:431545/1‑118 (MQ=255)
                                                                    tCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAGTAATCAGCCCAACCACGGa                                                                                 >  4:157160/1‑149 (MQ=255)
                                                                    tCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGa                                                                                 <  2:444443/149‑1 (MQ=255)
                                                                        ccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                                                                             >  3:151410/1‑149 (MQ=255)
                                                                         cTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGa                                                                            <  1:616281/149‑1 (MQ=255)
                                                                          tAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATc                                                                                                           <  2:629772/117‑1 (MQ=255)
                                                                          tAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAc                                                                           <  3:389783/149‑1 (MQ=255)
                                                                               ctctCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTcc                                                                                                                                   <  2:512259/88‑1 (MQ=255)
                                                                               ctctCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTcc                                                                                                                                   >  1:512259/1‑88 (MQ=255)
                                                                               ctctCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCt                                                                      >  4:506985/1‑149 (MQ=255)
                                                                                  tcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCca                                                                   >  1:733333/1‑149 (MQ=255)
                                                                                  tcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCca                                                                   >  4:443527/1‑149 (MQ=255)
                                                                                    cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagc                                                                 >  3:110047/1‑149 (MQ=255)
                                                                                      cGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagcag                                                               <  3:96555/149‑1 (MQ=255)
                                                                                              cGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAgg                                                       <  4:418710/149‑1 (MQ=255)
                                                                                                  gggACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCGTGCTAATAATCAGCACAAGAATCAGCCCAACCACGGAGTTGACCAGATCCAGCAGACCCAAGGtttt                                                   >  1:595757/1‑149 (MQ=255)
                                                                                                  gggACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGtttt                                                   >  1:398344/1‑149 (MQ=255)
                                                                                                  gggACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGtttt                                                   <  2:309840/149‑1 (MQ=255)
                                                                                                        gTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACgt                                             <  1:326401/149‑1 (MQ=255)
                                                                                                         tCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACgtg                                            >  1:103159/1‑149 (MQ=255)
                                                                                                            gggCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCt                                         >  1:137010/1‑149 (MQ=255)
                                                                                                            gggCACATTTTTAGACCTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCt                                         <  3:646986/149‑1 (MQ=255)
                                                                                                               cacaTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCtttt                                      >  3:342407/1‑149 (MQ=255)
                                                                                                                acaTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACTTAGTTGACCAGCTCCAGCAGCCCCCAGGTTTTCAACGTGTCTTTTa                                     <  4:624723/149‑1 (MQ=255)
                                                                                                                 caTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTAc                                    <  3:431267/149‑1 (MQ=255)
                                                                                                                   tttttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCaa                                                >  3:717015/1‑135 (MQ=255)
                                                                                                                    ttttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGa                                 <  4:603955/149‑1 (MQ=255)
                                                                                                                     tttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGAc                                <  3:179338/149‑1 (MQ=255)
                                                                                                                      ttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACa                               <  3:292709/149‑1 (MQ=255)
                                                                                                                            tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                         >  3:411791/1‑149 (MQ=255)
                                                                                                                            tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                         >  2:740190/1‑149 (MQ=255)
                                                                                                                            tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                         >  3:87974/1‑149 (MQ=255)
                                                                                                                            tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                         >  1:453468/1‑149 (MQ=255)
                                                                                                                            tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                         >  3:551996/1‑149 (MQ=255)
                                                                                                                             ttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaaa                        <  3:277472/149‑1 (MQ=255)
                                                                                                                             ttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaaa                        <  3:587864/149‑1 (MQ=255)
                                                                                                                             ttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaaa                        <  4:214567/149‑1 (MQ=255)
                                                                                                                                 catcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTaa                    <  1:169682/149‑1 (MQ=255)
                                                                                                                                  atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAc                   >  1:543977/1‑149 (MQ=255)
                                                                                                                                  atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAc                   >  1:169278/1‑149 (MQ=255)
                                                                                                                                   tcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAcc                  >  2:709145/1‑149 (MQ=255)
                                                                                                                                   tcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAcc                  >  2:50996/1‑149 (MQ=255)
                                                                                                                                      gTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCttt               <  1:379086/149‑1 (MQ=255)
                                                                                                                                       tCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTg              >  1:585005/1‑149 (MQ=255)
                                                                                                                                       tCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTg              >  4:721046/1‑149 (MQ=255)
                                                                                                                                        cagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGa             >  2:3786/3‑149 (MQ=255)
                                                                                                                                        cagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGa             >  4:746079/3‑149 (MQ=255)
                                                                                                                                        cagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGa             >  3:279257/3‑149 (MQ=255)
                                                                                                                                         tGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaa            <  4:391352/149‑1 (MQ=255)
                                                                                                                                          gAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaac           <  3:489365/149‑1 (MQ=255)
                                                                                                                                             ccTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaacaac        <  3:406045/149‑1 (MQ=255)
                                                                                                                                             ccTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaacaac        >  4:360351/1‑149 (MQ=255)
                                                                                                                                             ccTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaacaac        >  4:104294/1‑149 (MQ=255)
                                                                                                                                               tGCCATTGGAATGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGTACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCa      <  2:595757/149‑1 (MQ=255)
                                                                                                                                               tGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCa      >  1:401196/1‑149 (MQ=255)
                                                                                                                                                  cATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACGCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGGAACCTTTGAACAACCAGaa   >  2:546265/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAA  >  NC_000913/4549784‑4550079

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