Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,586,583 T→C D94G (GAT→GGT)  b1504 ← putative fimbrial‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,586,5830TC99.5% 693.2 / ‑4.9 199D94G (GAT→GGT) b1504putative fimbrial‑like protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (105/93);  total (105/94)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 204 total aligned reads displayed.

GGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTG  >  NC_000913/1586493‑1586672
                                                                                          |                                                                                         
 gACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCt                                                                                        >  2:5694736/1‑93 (MQ=255)
 gACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCt                                                                                        >  2:354833/1‑93 (MQ=255)
 gACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCt                                                                                        <  1:5051931/93‑1 (MQ=255)
 gACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTACTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCt                                                                                        >  2:6416960/1‑93 (MQ=255)
   cGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCa                                                                                      >  2:4894638/1‑93 (MQ=255)
   cGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCa                                                                                      >  2:1793843/1‑93 (MQ=255)
   cGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCa                                                                                      >  1:5967815/1‑93 (MQ=255)
    gCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCaa                                                                                     >  1:1546227/1‑93 (MQ=255)
      atatTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATa                                                                                   >  1:3386245/1‑93 (MQ=255)
        atTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATAcc                                                                                 >  1:5400136/1‑93 (MQ=255)
         tttGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCg                                                                                >  2:130442/1‑93 (MQ=255)
          ttGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGt                                                                               <  1:6167931/93‑1 (MQ=255)
          ttGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGt                                                                               <  2:395424/93‑1 (MQ=255)
           tGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACCAAGAATTATTTTCTGCACCCTCAATACCGtt                                                                              >  1:4803191/1‑93 (MQ=255)
           tGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGtt                                                                              >  1:589495/1‑93 (MQ=255)
              gTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTaaa                                                                           >  2:4777172/1‑93 (MQ=255)
               tCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAAt                                                                          <  2:5831106/93‑1 (MQ=255)
               tCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAAt                                                                          <  1:3096390/93‑1 (MQ=255)
                 cAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGc                                                                        >  1:2315669/1‑93 (MQ=255)
                  aaTATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCa                                                                       <  2:3187734/93‑1 (MQ=255)
                  aaTATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCa                                                                       >  1:1929751/1‑93 (MQ=255)
                     atTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCaaac                                                                    >  1:2827175/1‑91 (MQ=255)
                       ttCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCg                                                                  <  2:1287103/93‑1 (MQ=255)
                        tCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGa                                                                 >  1:3798991/1‑93 (MQ=255)
                        tCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGa                                                                 <  1:879347/93‑1 (MQ=255)
                         cTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGaa                                                                <  1:368812/93‑1 (MQ=255)
                         cTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGaa                                                                <  2:4976145/93‑1 (MQ=255)
                         cTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGaa                                                                >  2:2486031/1‑93 (MQ=255)
                          tataCCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAAc                                                               >  1:1032788/1‑93 (MQ=255)
                           ataCCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAAcc                                                              >  1:4326585/1‑93 (MQ=255)
                             aCCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAAcccc                                                            >  1:6101578/1‑93 (MQ=255)
                              cccAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGATccccc                                                           <  2:6275892/93‑1 (MQ=255)
                               ccAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCg                                                          <  2:2730502/93‑1 (MQ=255)
                                cagaCCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGt                                                         <  2:4009806/90‑1 (MQ=255)
                                cAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGt                                                         <  2:1144481/93‑1 (MQ=255)
                                cAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGt                                                         <  2:1116085/93‑1 (MQ=255)
                                 aaaCCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGtt                                                        <  1:812983/93‑1 (MQ=255)
                                 aaaCCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGtt                                                        <  1:5548056/93‑1 (MQ=255)
                                 aaaCCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGtt                                                        >  1:5809141/1‑93 (MQ=255)
                                  aagcGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTg                                                       <  2:2049336/90‑1 (MQ=255)
                                  aaCCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTg                                                       >  1:235966/1‑93 (MQ=255)
                                  aaCCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTg                                                       >  1:3583579/1‑93 (MQ=255)
                                   aCCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGt                                                      >  2:1484824/1‑93 (MQ=255)
                                    ccGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTc                                                     >  1:4541216/1‑93 (MQ=255)
                                    ccGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTc                                                     >  2:1563728/1‑93 (MQ=255)
                                     cGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTcc                                                    <  2:4117072/93‑1 (MQ=255)
                                      ggAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTccc                                                   >  1:4213720/1‑93 (MQ=255)
                                      ggAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTccc                                                   >  1:5210549/1‑93 (MQ=255)
                                       gAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCt                                                  >  2:3288240/1‑93 (MQ=255)
                                       gAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCt                                                  <  1:511131/93‑1 (MQ=255)
                                        aGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCtt                                                 >  1:593947/1‑93 (MQ=255)
                                        aGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCtt                                                 >  1:1727350/1‑93 (MQ=255)
                                         ggCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCttt                                                >  1:3362062/1‑93 (MQ=255)
                                           ccGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgc                                              <  1:3178538/93‑1 (MQ=255)
                                           ccGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgc                                              <  2:4297353/93‑1 (MQ=255)
                                           ccGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgc                                              <  1:4637568/93‑1 (MQ=255)
                                            cGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgct                                             >  1:2775574/1‑93 (MQ=255)
                                            cGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgct                                             >  1:1573139/1‑93 (MQ=255)
                                               aTTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCa                                          <  2:4484516/93‑1 (MQ=255)
                                                ttGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCAt                                         <  1:441014/93‑1 (MQ=255)
                                                  gCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTc                                       <  2:6120212/93‑1 (MQ=255)
                                                   cttccttcATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCa                                      >  1:1260147/1‑93 (MQ=255)
                                                    ttccttcATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCAc                                     <  1:2149790/93‑1 (MQ=255)
                                                      ccttcATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACtt                                   <  1:970215/93‑1 (MQ=255)
                                                      ccttcATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACtt                                   >  1:4530472/1‑93 (MQ=255)
                                                       cttcATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTa                                  <  1:134927/93‑1 (MQ=255)
                                                         tcATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAAc                                <  1:2496770/93‑1 (MQ=255)
                                                         tcATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAAc                                <  1:490113/93‑1 (MQ=255)
                                                         tcATCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAAc                                <  1:3633667/93‑1 (MQ=255)
                                                           aTCCAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGt                              >  2:5354395/1‑93 (MQ=255)
                                                             ccAGTTTCTACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTaa                            >  1:6406255/1‑93 (MQ=255)
                                                              cAGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAAt                           <  2:5993184/93‑1 (MQ=255)
                                                               aGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATc                          <  2:1189586/93‑1 (MQ=255)
                                                               aGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATc                          <  2:5754485/93‑1 (MQ=255)
                                                               aGTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATc                          >  2:121963/1‑93 (MQ=255)
                                                                gTTTCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCt                         <  1:6059885/93‑1 (MQ=255)
                                                                  ttCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGa                       >  2:3917257/1‑93 (MQ=255)
                                                                   tCAACAAAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGaa                      <  2:4668535/93‑1 (MQ=255)
                                                                    caacaaAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGaaa                     <  2:5413308/93‑1 (MQ=255)
                                                                     aacaaAGTATTATTTTCTGCCCCCTCAATACCGTTAAATGCAACCCGACCCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGaaaa                    <  1:1066249/93‑1 (MQ=255)
                                                                     aacaaAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGaaaa                    >  1:5682674/1‑93 (MQ=255)
                                                                      acaaAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAg                   >  1:6255948/1‑93 (MQ=255)
                                                                      acaaAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAg                   >  1:4000282/1‑93 (MQ=255)
                                                                      acaaAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAg                   >  1:5314800/1‑93 (MQ=255)
                                                                      acaaAGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAg                   >  2:2262996/1‑93 (MQ=255)
                                                                        aaaGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAggg                 >  2:845379/1‑93 (MQ=255)
                                                                        aaaGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAggg                 <  1:1674071/93‑1 (MQ=255)
                                                                        aaaGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAggg                 >  2:3109894/1‑93 (MQ=255)
                                                                         aaGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGa                <  2:373477/93‑1 (MQ=255)
                                                                          aGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGAc               <  2:5131377/93‑1 (MQ=255)
                                                                          aGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGAc               <  2:4815102/93‑1 (MQ=255)
                                                                          aGTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGAc               <  2:4249030/93‑1 (MQ=255)
                                                                           gTATTATTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACg              <  1:4133920/93‑1 (MQ=255)
                                                                            tattatTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACgg             <  1:5239079/93‑1 (MQ=255)
                                                                            tattatTTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACgg             <  1:1359089/93‑1 (MQ=255)
                                                                             attattTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGc            >  1:630513/1‑93 (MQ=255)
                                                                             attattTTCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGc            >  2:293099/1‑93 (MQ=255)
                                                                                  tttCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGAc       >  2:5578269/1‑93 (MQ=255)
                                                                                   ttCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACt      >  1:613977/1‑93 (MQ=255)
                                                                                    tCTGCACCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACtt     >  2:6078810/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCTGAAAAGGGACGGCTGGACTTGTG  >  NC_000913/1586493‑1586672

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