Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,314,800 A→C 5.6% Q438P (CAG→CCG) ‡ rcsD → phosphotransfer intermediate protein in two‑component regulatory system with RcsBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,314,8000AC 5.6% 261.7 / 7.1 71Q438P (CAG→CCG) ‡rcsDphosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base C (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

CCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACT  >  NC_000913/2314655‑2314947
                                                                                                                                                 |                                                                                                                                                   
ccGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAgg                                                                                                                                                   <  3:836137/148‑1 (MQ=255)
  gCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGAt                                                                                                                                                 <  4:197565/148‑1 (MQ=255)
       tGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGa                                                                                                                                            <  1:797866/148‑1 (MQ=255)
               cAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGAt                                                                                                                                                 <  2:835018/135‑1 (MQ=255)
               cAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGAt                                                                                                                                                 >  1:835018/1‑135 (MQ=255)
                   aCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGAt                                                                                                                                                 <  2:179355/131‑1 (MQ=255)
                   aCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGAt                                                                                                                                                 >  1:179355/1‑131 (MQ=255)
                    cATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAAc                                                                                                                               >  3:646440/1‑148 (MQ=255)
                    cATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAAc                                                                                                                               >  4:410114/1‑148 (MQ=255)
                     aTCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGAt                                                                                                                                                 >  4:402244/1‑129 (MQ=255)
                     aTCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGAt                                                                                                                                                 <  3:402244/129‑1 (MQ=255)
                     aTCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACa                                                                                                                              >  4:1089674/1‑148 (MQ=255)
                         ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGAt                                                                                                                                                 <  2:77918/125‑1 (MQ=255)
                         ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGAt                                                                                                                                                 >  1:77918/1‑125 (MQ=255)
                         ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaa                                                                                                                          <  3:410114/148‑1 (MQ=255)
                                  cGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCa                                                                                                                 >  1:1090793/1‑148 (MQ=255)
                                                gggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACa                                                                                                                              >  3:658312/1‑121 (MQ=255)
                                                gggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACa                                                                                                                              <  4:658312/121‑1 (MQ=255)
                                                 ggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCt                                                                                                     <  1:473792/145‑1 (MQ=255)
                                                 ggATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCt                                                                                                     >  2:473792/1‑145 (MQ=255)
                                                     tattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaa                                                                                                                          <  4:41356/120‑1 (MQ=255)
                                                     tattCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaa                                                                                                                          >  3:41356/1‑120 (MQ=255)
                                                        tCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaaa                                                                                           >  2:209621/1‑148 (MQ=255)
                                                          aGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaaaaa                                                                                         <  2:1090793/148‑1 (MQ=255)
                                                          aGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaaaaa                                                                                         <  1:209621/148‑1 (MQ=255)
                                                             cgacgaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACca                                                                                      <  4:646440/148‑1 (MQ=255)
                                                             cgacgaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACca                                                                                      >  4:1081857/1‑148 (MQ=255)
                                                                               tGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGa                                                                                                                                            >  2:106384/1‑76 (MQ=255)
                                                                               tGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGa                                                                                                                                            <  1:106384/76‑1 (MQ=255)
                                                                                 tATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACa                                                                                                                              <  2:407243/88‑1 (MQ=255)
                                                                                 tATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACa                                                                                                                              >  1:407243/1‑88 (MQ=255)
                                                                                 tATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                                  <  3:1089674/148‑1 (MQ=255)
                                                                                 tATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                                  <  3:1081857/148‑1 (MQ=255)
                                                                                  aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGc                                                                                                                                       >  1:1020380/1‑78 (MQ=255)
                                                                                  aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGc                                                                                                                                       <  2:1020380/78‑1 (MQ=255)
                                                                                  aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCa                                                                                                                     <  2:171201/96‑1 (MQ=255)
                                                                                  aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCa                                                                                                                     >  4:1090573/1‑96 (MQ=255)
                                                                                  aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCa                                                                                                                     <  3:1090573/96‑1 (MQ=255)
                                                                                  aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCa                                                                                                                     >  1:171201/1‑96 (MQ=255)
                                                                                    gagaTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGc                                                                                                         <  4:1126593/106‑1 (MQ=255)
                                                                                    gagaTCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGc                                                                                                         >  3:1126593/1‑106 (MQ=255)
                                                                                        tCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAgg                                                                                                               >  3:518850/1‑96 (MQ=255)
                                                                                        tCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAgg                                                                                                               <  4:518850/96‑1 (MQ=255)
                                                                                          cGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACa                                                         >  3:859819/1‑148 (MQ=255)
                                                                                            cATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACaa                                                                                                                             >  3:913804/1‑78 (MQ=255)
                                                                                            cATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACaa                                                                                                                             <  4:913804/78‑1 (MQ=255)
                                                                                            cATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAcc                                                                                       <  3:394899/116‑1 (MQ=255)
                                                                                            cATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAcc                                                                                       >  4:394899/1‑116 (MQ=255)
                                                                                            cATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGga                                                                           <  1:710154/128‑1 (MQ=255)
                                                                                            cATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGga                                                                           >  2:710154/1‑128 (MQ=255)
                                                                                               gTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                                  <  4:802762/134‑1 (MQ=255)
                                                                                               gTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                                  >  3:802762/1‑134 (MQ=255)
                                                                                                         cAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTaa                                                                                                                                 <  1:1107132/61‑1 (MQ=255)
                                                                                                         cAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTaa                                                                                                                                 >  2:1107132/1‑61 (MQ=255)
                                                                                                          aGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGaa                                         <  4:859819/148‑1 (MQ=255)
                                                                                                           ggTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATg                                                <  4:1117983/140‑1 (MQ=255)
                                                                                                           ggTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATg                                                >  3:1117983/1‑140 (MQ=255)
                                                                                                             tCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACAtt                                                       >  4:312405/1‑131 (MQ=255)
                                                                                                             tCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACAtt                                                       <  3:312405/131‑1 (MQ=255)
                                                                                                             tCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATg                                                >  2:629643/1‑138 (MQ=255)
                                                                                                             tCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATg                                                <  1:629643/138‑1 (MQ=255)
                                                                                                             tCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGa                                      >  2:1099438/1‑148 (MQ=255)
                                                                                                              cgcgcggcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                                      >  2:204268/1‑115 (MQ=255)
                                                                                                              cgcgcggcgcACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                                      <  1:204268/115‑1 (MQ=255)
                                                                                                                               ttttCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaaaa                                                                                          <  2:852235/78‑1 (MQ=255)
                                                                                                                               ttttCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGaaaa                                                                                          >  1:852235/1‑78 (MQ=255)
                                                                                                                                   cATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGgag                <  1:1099438/148‑1 (MQ=255)
                                                                                                                                         tCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATg                                                                                                <  3:48581/62‑1 (MQ=38)
                                                                                                                                         tCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATg                                                                                                >  4:48581/1‑62 (MQ=38)
                                                                                                                                         tCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg          >  2:96338/1‑148 (MQ=255)
                                                                                                                                         tCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGcgcg          >  1:1057665/1‑148 (MQ=255)
                                                                                                                                                 aGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACt  <  2:1057665/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCGCATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGAGCGCGGCTAAACT  >  NC_000913/2314655‑2314947

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