Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,323,124 A→C 5.3% L43R (CTC→CGC)  mrcB ← transglycosylase/transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,323,1240AC 5.3% 269.5 / 2.2 84L43R (CTC→CGC) mrcBtransglycosylase/transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (38/41);  new base C (3/2);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGTAGCCGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATTATACGC  >  NC_002947/5322982‑5323268
                                                                                                                                              |                                                                                                                                                 
aCGGTAGCCGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTaaa                                                                                                                                             <  8:363118/149‑1 (MQ=255)
       ccGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGaa                                                                                                                                      <  5:330400/149‑1 (MQ=255)
        ccccattaTCTAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACcgcg                                                                                                                                        <  1:348508/139‑1 (MQ=255)
               aTCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCaa                                                                                                                              >  4:41431/1‑149 (MQ=255)
               aTCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCaa                                                                                                                              >  5:162631/1‑149 (MQ=255)
               aTCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCaa                                                                                                                              >  5:412471/1‑149 (MQ=255)
                  gAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCTTCGCGGTTTACCTCGATGCCgtcgtcca                                                                                                                           <  7:220017/149‑8 (MQ=255)
                    gCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACcgcg                                                                                                                                        >  2:348508/1‑134 (MQ=255)
                              gAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCaa                                                                                                               >  6:338737/1‑149 (MQ=255)
                              gAAGTCGTGCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACTGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCaa                                                                                                               <  7:129252/149‑1 (MQ=255)
                               aaGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAAc                                                                                                              <  4:62988/149‑1 (MQ=255)
                                aGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAAc                                                                                                              >  6:111740/1‑148 (MQ=255)
                                aGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAAc                                                                                                              >  6:154030/1‑148 (MQ=255)
                                    gTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGc                                                                                                         >  3:56789/1‑149 (MQ=255)
                                       cTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTg                                                                                                       <  7:283766/148‑1 (MQ=255)
                                           cTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                                  >  8:258122/1‑149 (MQ=255)
                                           cTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                                  >  8:170452/1‑149 (MQ=255)
                                            tCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                                  >  8:332094/1‑148 (MQ=255)
                                            tCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                                  >  1:16129/1‑148 (MQ=255)
                                            tCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTc                                                                                                 >  4:272176/1‑149 (MQ=255)
                                            tCAACTTCTGGCCGGTGAACAGCTCCAGCGATCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                                  >  1:243751/1‑148 (MQ=255)
                                            tCAAATTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTGCAGCGGTTGCCGGAGAACTTTTCCGGGACGACAGCTTCGCGTTTAACCGCGCACCCGGCGACGATCACCAGGCCAGCCAGGCTgcgcgt                                                                                                  >  6:327568/1‑146 (MQ=255)
                                              aaCTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCAc                                                                                                                       <  8:89840/125‑1 (MQ=255)
                                              aaCTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCAc                                                                                                                       >  7:89840/1‑125 (MQ=255)
                                                   cTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagcc                                                                                           <  5:317850/148‑1 (MQ=255)
                                                      gCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagcccagc                                                                                       >  4:209194/1‑149 (MQ=255)
                                                      gCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagccca                                                                                         <  6:41951/147‑1 (MQ=255)
                                                        cGGTGCACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCACCGATCACCAGGCCAACCAGCCAGAGCTTCAAcgcccagcc                                                                                      >  2:86583/1‑148 (MQ=255)
                                                        cGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagcccagcc                                                                                      >  3:181976/1‑148 (MQ=255)
                                                         ggTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAgcccagccca                                                                                    <  5:405875/149‑1 (MQ=255)
                                                          gTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGcccagcccag                                                                                   >  3:147204/1‑149 (MQ=255)
                                                                  cTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGgcgc                                                                          <  4:147204/149‑1 (MQ=255)
                                                                     cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGcgcg                                                                         <  3:199471/147‑1 (MQ=255)
                                                                     cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCgg                                                                        <  5:61292/148‑1 (MQ=255)
                                                                     cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCgg                                                                        <  4:181976/148‑1 (MQ=255)
                                                                     cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCgg                                                                        <  6:412471/148‑1 (MQ=255)
                                                                          gTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCgg                                                                  <  2:238902/149‑1 (MQ=255)
                                                                           tCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGc                                                                 <  3:272176/149‑1 (MQ=255)
                                                                           tCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGc                                                                 >  2:364542/1‑149 (MQ=255)
                                                                           tCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGc                                                                 <  6:162631/149‑1 (MQ=255)
                                                                           tCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGc                                                                 <  7:170452/149‑1 (MQ=255)
                                                                             gggCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCg                                                               >  8:125376/1‑149 (MQ=255)
                                                                                 aTACACTTTGGCCGGGATGGTCCAGCGCTTGTCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCg                                                           >  3:165360/1‑149 (MQ=255)
                                                                                 aTACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCg                                                           >  4:301207/1‑149 (MQ=255)
                                                                                  tACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCgg                                                          <  5:154030/149‑1 (MQ=255)
                                                                                  tACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCgg                                                          <  2:243751/149‑1 (MQ=255)
                                                                                   acacTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCggg                                                         >  8:144211/1‑149 (MQ=255)
                                                                                     acTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCg                                                       >  7:361852/1‑149 (MQ=255)
                                                                                         tGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGttt                                                    <  3:209194/148‑1 (MQ=255)
                                                                                              gggATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAgg                                              >  5:51661/1‑149 (MQ=255)
                                                                                               ggATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAggg                                             <  3:301207/149‑1 (MQ=255)
                                                                                               ggATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAggg                                             >  8:125232/1‑149 (MQ=255)
                                                                                                 aTGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGt                                           <  3:41431/149‑1 (MQ=255)
                                                                                                 aTGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGAGAGCGGCCGGTCGGGCGTTTCTGAGGGGt                                           <  3:181187/149‑1 (MQ=255)
                                                                                                   ggTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGc                                         >  2:329150/1‑149 (MQ=255)
                                                                                                       cAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGc                                                                                                                  <  3:110480/73‑1 (MQ=255)
                                                                                                               gCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgg                              <  8:381855/148‑1 (MQ=255)
                                                                                                               gCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgg                              <  5:111740/148‑1 (MQ=255)
                                                                                                               gCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcg                               <  7:332094/147‑1 (MQ=255)
                                                                                                               gCCGGAGAACTGTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgg                              <  8:269054/148‑1 (MQ=255)
                                                                                                                ccGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcggg                             <  1:161810/148‑1 (MQ=255)
                                                                                                                 cGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCATCAAAGCCCAGCCCAGCAAAGCGGCGGGACCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggt                            >  1:131354/1‑149 (MQ=255)
                                                                                                                    agaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGGGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcgg                        <  5:338737/149‑1 (MQ=255)
                                                                                                                     gaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcggg                       >  1:190714/1‑149 (MQ=255)
                                                                                                                     gaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcggg                       >  7:9595/1‑149 (MQ=255)
                                                                                                                     gaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcgg                        <  1:364542/148‑1 (MQ=255)
                                                                                                                       aCTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCTAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTc                     <  7:258122/149‑1 (MQ=255)
                                                                                                                       aCTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTc                     <  7:144211/149‑1 (MQ=255)
                                                                                                                         ttttCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCAt                   >  5:186327/1‑149 (MQ=255)
                                                                                                                         ttttCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCAt                   >  2:119147/1‑149 (MQ=255)
                                                                                                                          tttCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATg                  >  4:312380/1‑149 (MQ=255)
                                                                                                                          tttCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATg                  >  7:107970/1‑149 (MQ=255)
                                                                                                                              cTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTGAGGTTCGGGTCATggcgg              >  1:138990/1‑149 (MQ=255)
                                                                                                                                ggACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCggg                                      >  5:186998/1‑123 (MQ=255)
                                                                                                                                ggACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCggg                                      <  6:186998/123‑1 (MQ=255)
                                                                                                                                 gacgacGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGgcggcg            <  2:354000/148‑1 (MQ=255)
                                                                                                                                 gacgacGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGcggcgg           <  3:312380/149‑1 (MQ=255)
                                                                                                                                  acgacgGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGa          <  7:125232/149‑1 (MQ=255)
                                                                                                                                    gacgGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGAtt        <  1:119147/149‑1 (MQ=255)
                                                                                                                                     ccgGCTTCGCGGTTTACCTCGATGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcggg                       <  5:217062/132‑1 (MQ=255)
                                                                                                                                     ccgGCTTCGCGGTTTACCTCGATGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcggg                       >  6:217062/2‑133 (MQ=255)
                                                                                                                                     ccgGCTTCGCGGTTTACCTCGATGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGAGGGGAGTTCGGGTTCGGGTCATGGCGGCGGATta       >  8:351027/2‑149 (MQ=255)
                                                                                                                                         cATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATTATACg   <  4:165360/149‑1 (MQ=255)
                                                                                                                                          aTCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATTATACGc  >  7:68115/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
ACGGTAGCCGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATTATACGC  >  NC_002947/5322982‑5323268

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