Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,289,431 T→C G150G (GGT→GGC yraK → putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,289,4310TC98.1% 645.5 / ‑2.1 208G150G (GGT→GGCyraKputative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (111/93);  total (113/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

CGACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGTGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAT  >  NC_000913/3289341‑3289520
                                                                                          |                                                                                         
cGACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGa                                                                                         <  2:2962205/93‑1 (MQ=255)
cGACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGa                                                                                         <  1:1516476/93‑1 (MQ=255)
cGACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGa                                                                                         >  2:1091826/1‑93 (MQ=255)
 gACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTTCCATGGCGAt                                                                                        >  2:3609705/1‑93 (MQ=255)
  aCCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATa                                                                                       >  2:2562369/1‑93 (MQ=255)
  aCCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATa                                                                                       <  2:4611634/93‑1 (MQ=255)
   ccAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATAc                                                                                      >  2:1927197/1‑93 (MQ=255)
   ccAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATAc                                                                                      <  1:3920939/93‑1 (MQ=255)
       gTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtcc                                                                                  >  2:2936193/1‑92 (MQ=255)
         tGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCg                                                                                >  1:1377220/1‑93 (MQ=255)
          gTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGt                                                                               <  2:4429911/93‑1 (MQ=255)
           tctcTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGtt                                                                              <  1:4355159/93‑1 (MQ=255)
            ctctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTa                                                                             >  2:1992211/1‑93 (MQ=255)
            ctctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTa                                                                             <  1:2611551/93‑1 (MQ=255)
             tctGTTAACACCCGATTCACACAGATAAAAAACAACACTGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTAt                                                                            >  1:4743275/1‑93 (MQ=255)
             tctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTAt                                                                            <  2:3803564/93‑1 (MQ=255)
             tctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTAt                                                                            <  2:1471654/93‑1 (MQ=255)
              ctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATg                                                                           <  1:2828311/93‑1 (MQ=255)
              ctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATg                                                                           <  2:4513502/93‑1 (MQ=255)
              ctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATg                                                                           >  1:383852/1‑93 (MQ=255)
               tGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATgg                                                                          <  1:2699789/93‑1 (MQ=255)
               tGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATgg                                                                          <  2:1149540/93‑1 (MQ=255)
                gTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGa                                                                         >  1:3855017/1‑93 (MQ=255)
                  tAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGAcc                                                                       <  2:4101876/93‑1 (MQ=255)
                    acacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGTGATACGTCTCGTTATGGAccag                                                                     >  2:4041769/1‑91 (MQ=255)
                    acacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCg                                                                     >  2:3733921/1‑93 (MQ=255)
                    acacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCg                                                                     <  2:3928339/93‑1 (MQ=255)
                     cacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGt                                                                    >  1:3705060/1‑93 (MQ=255)
                      ccccGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGtt                                                                   <  2:3225820/92‑1 (MQ=255)
                      acCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGtt                                                                   >  2:815067/1‑93 (MQ=255)
                       cccGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTa                                                                  <  2:1245241/93‑1 (MQ=255)
                       cccGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTa                                                                  <  1:98572/93‑1 (MQ=255)
                       cccGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCAAGGGCGATACGTCTCGTTATGGACCCGTTa                                                                  >  2:3755355/1‑93 (MQ=255)
                        ccGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTGGCCAGGGCGATACGTCTCGTTATGGACCCGTTaa                                                                 <  1:185316/93‑1 (MQ=255)
                          gATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAata                                                               >  2:4087216/1‑93 (MQ=255)
                          gATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAata                                                               >  2:2214047/1‑93 (MQ=255)
                            ttCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGACAGGGCGATACGTCTCGTTATGGACCCGTTAATAtt                                                             >  1:412580/1‑93 (MQ=255)
                             tCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTg                                                            >  1:2062050/1‑93 (MQ=255)
                              cacacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTgg                                                           <  1:2422541/93‑1 (MQ=255)
                               acacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGt                                                          <  1:1143056/93‑1 (MQ=255)
                               acacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGGTAATATTGGt                                                          >  2:4356483/1‑93 (MQ=255)
                                 acaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTgg                                                        >  1:3134454/1‑92 (MQ=255)
                                 acaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGc                                                        >  2:2922191/1‑93 (MQ=255)
                                   aGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGa                                                      >  2:2435373/1‑93 (MQ=255)
                                   aGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGAGTTTGCCAGGGCGATACGTCTCGGTATGCACCCGTTAATATTGCTGCGa                                                      >  1:652059/1‑93 (MQ=255)
                                    gATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGaa                                                     >  2:3548708/1‑93 (MQ=255)
                                      tAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAaca                                                   >  2:3429523/1‑93 (MQ=255)
                                        aaaaaCAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAAcacc                                                 <  2:255661/93‑1 (MQ=255)
                                          aaaCAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAAcaccac                                               >  1:4641586/1‑93 (MQ=255)
                                           aacaacACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACaccacc                                              <  1:2477302/93‑1 (MQ=255)
                                           aacaacACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACaccacc                                              >  2:285947/1‑93 (MQ=255)
                                            acaacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCt                                             >  2:3284525/1‑93 (MQ=255)
                                             caacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCtt                                            <  2:3892975/93‑1 (MQ=255)
                                              aacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGTGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCttt                                           >  1:4103801/1‑93 (MQ=255)
                                              aacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCttt                                           <  2:1121720/93‑1 (MQ=255)
                                                 acGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACATCACCTTTAcc                                        >  2:558331/1‑93 (MQ=255)
                                                 acGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTAcc                                        <  1:1477651/93‑1 (MQ=255)
                                                  cGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTAccc                                       >  2:3794526/1‑93 (MQ=255)
                                                  cGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTAccc                                       <  1:519823/93‑1 (MQ=255)
                                                  cGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTAccc                                       <  1:3019577/93‑1 (MQ=255)
                                                    gTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTg                                     >  2:3675894/1‑93 (MQ=255)
                                                    gTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTg                                     >  1:2336433/1‑93 (MQ=255)
                                                      atatTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATGCGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTa                                   <  2:706087/93‑1 (MQ=255)
                                                      atatTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTa                                   <  2:2761359/93‑1 (MQ=255)
                                                       tatTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTAt                                  >  2:1588866/1‑93 (MQ=255)
                                                       tatTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTAt                                  <  1:38701/93‑1 (MQ=255)
                                                           ctgctgCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCa                              >  2:433733/1‑93 (MQ=255)
                                                           ctgctgCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCa                              >  1:4724572/1‑93 (MQ=255)
                                                            tgctgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCAc                             <  2:3172386/93‑1 (MQ=255)
                                                            tgctgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCAc                             >  2:711517/1‑93 (MQ=255)
                                                            tgctgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATCGACCCGTTAATATTGGTGCGAACACCACCTTTACCCCGTATGTCAc                             >  2:3566723/1‑93 (MQ=255)
                                                             gctgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCAcc                            >  2:125427/1‑93 (MQ=255)
                                                               tgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCac                          >  2:1692312/1‑92 (MQ=255)
                                                               tgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCaa                          >  2:1975190/1‑93 (MQ=255)
                                                                  aaCCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGcc                       <  2:2473218/93‑1 (MQ=255)
                                                                   aCCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCa                      >  2:1059414/1‑93 (MQ=255)
                                                                    cccAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAt                     >  2:2428428/1‑93 (MQ=255)
                                                                    cccAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAt                     >  1:2142306/1‑93 (MQ=255)
                                                                     ccAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAtt                    <  1:1448597/93‑1 (MQ=255)
                                                                      cAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCCCCTTTACCCTGTATGTCACCAAGCCAttt                   <  1:3527598/93‑1 (MQ=255)
                                                                      cAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAttt                   >  2:757912/1‑93 (MQ=255)
                                                                      cAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAttt                   >  2:3575746/1‑93 (MQ=255)
                                                                      cAAACGGGTGTTTGCCAGGGCGATACGTCTCGGTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAttt                   >  1:3748729/1‑93 (MQ=255)
                                                                        aaCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtct                 >  1:3633242/1‑93 (MQ=255)
                                                                             gtgtTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCt            <  1:3251946/93‑1 (MQ=255)
                                                                             gtgtTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCt            >  1:4664829/1‑93 (MQ=255)
                                                                              tgtTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTc           <  2:4471926/93‑1 (MQ=255)
                                                                              tgtTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTc           >  1:3198915/1‑93 (MQ=255)
                                                                              tgtTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTc           >  1:1311257/1‑93 (MQ=255)
                                                                              tgtTTGACGGGGCGATACGTCTCGTTATGGACCCGTTAATCTTGGTGCGAACCCCCCCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTc           <  2:4225693/93‑1 (MQ=255)
                                                                                 ttGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGAt        <  1:1775271/93‑1 (MQ=255)
                                                                                  tGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATg       >  2:2900750/1‑93 (MQ=255)
                                                                                    ccAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGAc     <  1:3881651/93‑1 (MQ=255)
                                                                                      aGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGGCCa   >  2:2526999/1‑93 (MQ=255)
                                                                                      aGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCa   >  1:4547314/1‑93 (MQ=255)
                                                                                      aGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCa   <  1:3007850/93‑1 (MQ=255)
                                                                                      aGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCa   >  2:1872636/1‑93 (MQ=255)
                                                                                      aGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCa   >  2:1951770/1‑93 (MQ=255)
                                                                                      aGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCa   <  2:262139/93‑1 (MQ=255)
                                                                                       gggCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAt  <  1:947178/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGTGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAT  >  NC_000913/3289341‑3289520

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