Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,525,756 A→G 6.5% N211S (AAT→AGT)  yfeN → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,525,7560AG 6.5% 176.4 / 2.6 61N211S (AAT→AGT) yfeNputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base G (2/2);  total (29/32)
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.00e-01

CTGGCGACAACGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGG  >  NC_000913/2525609‑2525903
                                                                                                                                                   |                                                                                                                                                   
cTGGCGACAACGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAAt                                                                                                                                                    <  1:863311/149‑1 (MQ=255)
 tGGCGACAACGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATg                                                                                                                                                   >  1:511426/1‑149 (MQ=255)
   gCGACAACGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGgc                                                                                                                                                 >  1:1259711/1‑149 (MQ=255)
        atcGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGt                                                                                                                                            <  2:1231626/147‑1 (MQ=255)
         aCGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGGTGGGGCGAGGAAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGtt                                                                                                                                           <  2:1027684/149‑1 (MQ=255)
         aCGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGtt                                                                                                                                           <  1:1008878/149‑1 (MQ=255)
         aCGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGtt                                                                                                                                           <  1:501551/149‑1 (MQ=255)
         aCGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGtt                                                                                                                                           <  1:501345/149‑1 (MQ=255)
          cGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTg                                                                                                                                          <  2:149236/149‑1 (MQ=255)
          cGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTGTATGCTGGGCAGGAAGTAGTTCAGCCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTg                                                                                                                                          <  1:278436/149‑1 (MQ=255)
             cTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAAAGAGTACGAGTTTGACCGTGACGCTACCTGTGCAGCGGGGAATGGCGGTAACGAAGGGCTGAATGTCGCAGTTGCAc                                                                                                                                       >  1:513054/1‑149 (MQ=255)
             cTATGTCGCCGGCTGGGTTGCAGGCTACAACTTTATGCTGGGCAGCGAGAAATTCACTGTTACCAACTGGACAGGGTACGAGATTGACCGTGACGCTCCCTAGCCAGCGGGGAATGGCTGTAAAGATGGGCTGAATGGCGCAGCTGccc                                                                                                                                       >  2:328778/1‑147 (MQ=255)
                  tCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAGTGGCGGGAAAGAAGGGCTGAATGGCGCAGTTGAACTgtgg                                                                                                                                  >  1:358675/1‑149 (MQ=255)
                      cGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGa                                                                                                                              <  2:1303044/149‑1 (MQ=255)
                       ggCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGCTGGaa                                                                                                                             >  2:1085352/1‑149 (MQ=255)
                            ggTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCaa                                                                                                                        <  2:1427206/149‑1 (MQ=255)
                              ttGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACa                                                                                                                      <  1:911068/149‑1 (MQ=255)
                               tGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTGCGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGTCGCAGTTGCACTGTGGTGGAATGCAACAt                                                                                                                     >  2:1072592/1‑149 (MQ=255)
                                 cAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCAGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATca                                                                                                                   <  2:971759/149‑1 (MQ=255)
                                   ggTTACAACTTTATGCTGGGCAGCAAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATcaca                                                                                                                 <  2:261770/149‑1 (MQ=255)
                                    gTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAGTGGCGCAGTTGCACTGTGGTGGAATGCAACATcacac                                                                                                                >  2:1426354/1‑149 (MQ=255)
                                    gTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAGTGGCGCAGTTGCACTGTGGTGGAATGCAACATcacac                                                                                                                >  2:770087/1‑149 (MQ=255)
                                         aaCTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGTTGGAATGCAACATCACACATtac                                                                                                           >  1:389735/1‑149 (MQ=255)
                                          aCTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATtact                                                                                                          <  1:1413148/149‑1 (MQ=255)
                                           cTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATtacta                                                                                                         >  1:623340/1‑149 (MQ=255)
                                            tttATGCCGGGGAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATtactac                                                                                                        <  1:660584/149‑1 (MQ=255)
                                                 gCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGa                                                                                                   <  2:839683/149‑1 (MQ=255)
                                                  cTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGAt                                                                                                  >  2:1114111/1‑149 (MQ=255)
                                                            agaAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATAGa                                                                                        >  1:468697/1‑149 (MQ=255)
                                                             gaAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGat                                                                                       <  1:19747/149‑1 (MQ=255)
                                                               aaTTCATTCTTCGCACCTGGAGCGAGGACGAGTTGGCTCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGGACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGatat                                                                                     <  2:1304397/149‑1 (MQ=255)
                                                                aTTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCGGTCTCGATCTg                                                                                    >  2:1138257/1‑149 (MQ=255)
                                                                 ttCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTCAAGAAGGGCTGAATGGCGCAGTGGCACTGTGGGGGAATGCAACCTCACACATTACCACGGGGATTCAGTAGCGCTATGc                                                                                   >  1:1244335/1‑149 (MQ=255)
                                                                 ttCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAAGGCTGAATGGCGCAGTTGCAGTTTGATGGAATGCAGCATCACCCATTCCTACGGGGATTCATTATCGatatag                                                                                   >  2:520195/1‑147 (MQ=255)
                                                                    aCTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGa                                                                                >  1:1250668/1‑149 (MQ=255)
                                                                        ttACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATGGAGATGCGGATGAc                                                                            >  1:2748/1‑149 (MQ=255)
                                                                         tACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACa                                                                           <  2:560594/149‑1 (MQ=255)
                                                                                 ggAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACACCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGc                                                                   <  1:1138257/149‑1 (MQ=255)
                                                                                        gtacgaggTTGACCGTGACGTTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGAtt                                                            <  1:610549/149‑1 (MQ=255)
                                                                                         tacgagtTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGAttt                                                           >  1:597081/1‑149 (MQ=255)
                                                                                         taagagtTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGAttt                                                           >  2:1115452/4‑149 (MQ=255)
                                                                                            gagtTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACCTCACACATTACTACGGGGATTCAGTATCGATCTGCGGATGACAAACTCGGCGAAGATTTTTa                                                        >  1:47616/1‑149 (MQ=255)
                                                                                                 tGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAgg                                                   <  1:640289/149‑1 (MQ=255)
                                                                                                  gACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGa                                                  <  2:706775/149‑1 (MQ=255)
                                                                                                       tGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGa                                             >  1:777852/1‑149 (MQ=255)
                                                                                                       tGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGAGGATGACAAACTGGGCGAAGATTTTTACCAGGCTGCGa                                             >  1:1264068/1‑149 (MQ=255)
                                                                                                        gACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGAt                                            <  1:371120/149‑1 (MQ=255)
                                                                                                            cTACCTAGGCAGCGGGTAATGGCGCTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATc                                        <  1:1207240/149‑1 (MQ=255)
                                                                                                             tACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGCTAATCt                                       >  1:964094/1‑149 (MQ=255)
                                                                                                                  aTGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGGTGCGATAATATACTcc                                  >  2:655932/1‑149 (MQ=255)
                                                                                                                    gCAGCGGGTAATGGCGGTAAAGAAGGGCTGAGTGGCGCAGTTGCACTGTGGTGGACTGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCCACTCCAt                                <  1:1426354/149‑1 (MQ=255)
                                                                                                                    gCAGCGGGTAATGGCGGTAAAGAAGGGCTGAGTGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCCACTCCAt                                <  1:770087/149‑1 (MQ=255)
                                                                                                                           gTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATATACTCCATCAAAttt                         <  2:1260421/149‑1 (MQ=255)
                                                                                                                             aaTGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTaa                       <  1:1114111/149‑1 (MQ=255)
                                                                                                                                   ggTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCTTCAACTTTAACTTCtg                 >  1:368435/1‑149 (MQ=255)
                                                                                                                                   ggTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGAGGCGCTAAACTACTCCAGCAAAATTAACATCtg                 >  2:343241/1‑149 (MQ=255)
                                                                                                                                      aaaGGAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCtgatt              <  2:504279/149‑1 (MQ=255)
                                                                                                                                        agaagGGGTGAATGGCGCAGTTGCGCTGTGGTGGAACGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCtgattga            <  2:511426/149‑1 (MQ=255)
                                                                                                                                        agaagGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACCCATTACTACGGGGAGTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCtgattga            <  2:506380/149‑1 (MQ=255)
                                                                                                                                          aagGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACCTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTaattgatt          >  1:939633/1‑149 (MQ=255)
                                                                                                                                          aagGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTgattgatt          >  1:1308912/1‑149 (MQ=255)
                                                                                                                                                  aaTGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCgg  <  1:128711/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGGCGACAACGGCTATGTCGCCGGCTGGGTTGCAGGTTACAACTTTATGCTGGGCAGCGAGAAATTCACTCTTACCAACTGGAACGAGTACGAGTTTGACCGTGACGCTACCTATGCAGCGGGTAATGGCGGTAAAGAAGGGCTGAATGGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGG  >  NC_000913/2525609‑2525903

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