Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 150,947 T→C T3A (ACT→GCT)  yadC ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913150,9470TC98.9% 298.4 / ‑5.7 88T3A (ACT→GCT) yadCputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (45/42);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AGGTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGATTTTT  >  NC_000913/150812‑151078
                                                                                                                                       |                                                                                                                                   
aGGTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTTAAGATAGCCtt                                                                                                                                  >  2:954403/1‑139 (MQ=255)
 ggTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGACTGTCtct                                                                                                                                 >  2:705304/1‑137 (MQ=255)
     aTCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAttt                                                                                                                             >  1:1291919/1‑139 (MQ=255)
      tCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAtttt                                                                                                                            >  2:1145040/1‑139 (MQ=255)
        tgtgCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTa                                                                                                                          <  1:1145040/139‑1 (MQ=255)
            cGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTcc                                                                                                                      <  2:601348/139‑1 (MQ=255)
            cGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTcc                                                                                                                      <  1:204982/139‑1 (MQ=255)
                aGTGTAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATg                                                                                                                  <  2:196457/139‑1 (MQ=255)
                  tGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGtt                                                                                                                >  2:1414628/1‑139 (MQ=255)
                    gAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTc                                                                                                              >  1:485585/1‑139 (MQ=255)
                    gAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTc                                                                                                              >  1:1218840/1‑139 (MQ=255)
                           cTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTACTg                                                                                                       >  1:973104/1‑139 (MQ=255)
                            tCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGt                                                                                                      >  1:838018/1‑139 (MQ=255)
                                     aTCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTg                                                                                             <  2:116494/139‑1 (MQ=255)
                                      tCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACCGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGTTAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgc                                                                                            <  1:954791/139‑1 (MQ=255)
                                       cAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCGTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                           <  1:617888/139‑1 (MQ=255)
                                               gCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGttttt                                                                                   >  1:355898/1‑139 (MQ=255)
                                                 attattTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                                 <  2:973104/139‑1 (MQ=255)
                                                     tttCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATg                                                                                                                  <  1:1048698/102‑1 (MQ=255)
                                                     tttCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATg                                                                                                                  >  2:1048698/1‑102 (MQ=255)
                                                       tCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                           >  1:313817/1‑123 (MQ=255)
                                                       tCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                           <  2:313817/123‑1 (MQ=255)
                                                       tCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTa                                                                           >  1:747799/1‑139 (MQ=255)
                                                          cacaATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTc                                                                        >  2:794713/1‑139 (MQ=255)
                                                             aaTAAAACTTGAATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTATTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgta                                                                     >  1:1429325/1‑139 (MQ=255)
                                                                    cTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                              <  2:485585/139‑1 (MQ=255)
                                                                      tGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATc                                                            <  1:1372910/139‑1 (MQ=255)
                                                                       gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCCTTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATca                                                           >  1:603921/1‑138 (MQ=255)
                                                                       gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                          <  2:381882/92‑1 (MQ=255)
                                                                       gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                          >  1:381882/1‑92 (MQ=255)
                                                                       gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGtt                                                                                      <  2:1079732/112‑1 (MQ=255)
                                                                       gTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGtt                                                                                      >  1:1079732/1‑112 (MQ=255)
                                                                               actgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgt                                                                      <  1:713522/120‑1 (MQ=255)
                                                                               actgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgt                                                                      >  2:713522/1‑120 (MQ=255)
                                                                                ctgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAAcaca                                                  <  2:355898/139‑1 (MQ=255)
                                                                                 tgactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACAt                                                 <  2:1218840/139‑1 (MQ=255)
                                                                                  gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTa                                                                                                                          <  1:897011/65‑1 (MQ=255)
                                                                                  gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                          >  2:563810/1‑81 (MQ=255)
                                                                                  gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                          <  1:563810/81‑1 (MQ=255)
                                                                                  gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATg                                                >  1:977099/1‑139 (MQ=255)
                                                                                  gactgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATg                                                >  1:1299697/1‑139 (MQ=255)
                                                                                  gactgTATTACAACAAGAATACAGAGCTAAAAAAAGAATGAACCTGAAGATAGCCTTCATTTTTa                                                                                                                          >  2:897011/1‑65 (MQ=255)
                                                                                   actgTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATgg                                               <  2:353562/139‑1 (MQ=255)
                                                                                      gTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTc                                            <  1:94534/139‑1 (MQ=255)
                                                                                       tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAttt                                                                                                                             >  1:873428/1‑57 (MQ=255)
                                                                                       tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAttt                                                                                                                             <  2:873428/57‑1 (MQ=255)
                                                                                       tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCt                                                                                         <  2:1438914/93‑1 (MQ=255)
                                                                                       tATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCt                                                                                         >  1:1438914/1‑93 (MQ=255)
                                                                                        aTTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGttt                                                                                                               <  2:638723/70‑1 (MQ=255)
                                                                                        aTTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGttt                                                                                                               >  1:638723/1‑70 (MQ=255)
                                                                                        aTTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                                 >  2:21837/1‑100 (MQ=255)
                                                                                        aTTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                                 <  1:21837/100‑1 (MQ=255)
                                                                                         ttACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATa                                         <  1:954403/139‑1 (MQ=255)
                                                                                              acaaGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                          >  2:273048/1‑69 (MQ=255)
                                                                                              acaaGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTa                                                                                                          <  1:273048/69‑1 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTc                                                                                                                       <  2:321747/53‑1 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTc                                                                                                                       >  1:321747/1‑53 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                                 >  1:800469/1‑91 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgt                                                                                 <  2:800469/91‑1 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                              >  2:1004896/1‑110 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                                 >  1:200436/1‑139 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                                 <  2:1299697/139‑1 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                                 >  1:717447/1‑139 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGAcc                                 >  2:1148428/1‑139 (MQ=255)
                                                                                                 aGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCAATTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAata                                                              <  1:1004896/110‑1 (MQ=255)
                                                                                                   aaTACAGTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTCCTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAg                               <  2:289106/139‑1 (MQ=255)
                                                                                                         gTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgt                                                                      >  1:319990/1‑94 (MQ=255)
                                                                                                         gTGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCgt                                                                      <  2:319990/94‑1 (MQ=255)
                                                                                                          tGCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAac                        >  1:910408/1‑139 (MQ=255)
                                                                                                           gCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                     >  2:98697/1‑61 (MQ=255)
                                                                                                           gCTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGtt                                                                                                     <  1:98697/61‑1 (MQ=255)
                                                                                                            cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAacac                      >  2:1460959/1‑139 (MQ=255)
                                                                                                            cTAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAacac                      >  2:26415/1‑139 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGtttt                                                                                    <  2:107921/76‑1 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGtttt                                                                                    >  1:107921/1‑76 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt                     >  2:1021119/1‑139 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt                     <  2:717447/139‑1 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt                     <  2:767209/139‑1 (MQ=255)
                                                                                                             tAAAAAAAGAATGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACt                     >  1:767209/1‑139 (MQ=255)
                                                                                                                      aaTGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                           <  1:1077749/60‑1 (MQ=255)
                                                                                                                      aaTGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTgcg                                                                                           >  2:1077749/1‑60 (MQ=255)
                                                                                                                        tGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGAt          <  1:1460959/139‑1 (MQ=255)
                                                                                                                        tGTACCTGAAGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGAt          <  1:1414628/139‑1 (MQ=255)
                                                                                                                                aaGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgtttg                                                                              >  1:919719/1‑63 (MQ=255)
                                                                                                                                aaGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTtttgtttg                                                                              <  2:919719/63‑1 (MQ=255)
                                                                                                                                aaGATAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGAttttt  >  2:621967/1‑139 (MQ=255)
                                                                                                                                    tAGCCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAaca                       >  1:942282/1‑114 (MQ=255)
                                                                                                                                    tAGCCTTCACTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAaca                       <  2:942282/114‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
AGGTCATCTGTGCGGCAGTGGAGGGGCCTCGATAGTTATCAACTCCTGCATTATTTCCCACAATAAAACTTGTATATGCACTGACTGTATTACAACAAGAATACAGTGCTAAAAAAAGAATGTACCTGAAGATAGTCTTCATTTTTACTCCAATGTTTCTTTATTGTTATGGTTTGCGCTGTTTTTGTTTGTTATTCGTAGTAAATATCCACTAACACATGGCTCATAAGTTGACCAGATGTAACACTCAATGCAGGATCGATTTTT  >  NC_000913/150812‑151078

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