Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 194,608 A→C 12.5% E38D (GAA→GAC) ‡ PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947194,6080AC12.5% 169.5 / 16.6 64E38D (GAA→GAC) ‡PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base C (3/5);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

AAGTACCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCCG  >  NC_002947/194464‑194744
                                                                                                                                                |                                                                                                                                        
aaGTACCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCAGTTGGCCGGTGAGGAAGTGc                                                                                                                                      >  5:284728/1‑149 (MQ=255)
  gTACCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTc                                                                                                                                    <  3:229702/149‑1 (MQ=255)
    aCCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCAc                                                                                                                                  >  2:350574/1‑149 (MQ=255)
     ccccATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCAc                                                                                                                                  <  5:91345/148‑1 (MQ=255)
      cccATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCCCCg                                                                                                                                >  6:6305/1‑149 (MQ=255)
       ccATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGCACTTCCTCACCgg                                                                                                                               <  1:212783/149‑1 (MQ=255)
        cATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGCACTTCCTCACCggc                                                                                                                              <  7:53543/149‑2 (MQ=255)
             ggTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCggg                                                                                                                         <  3:123844/149‑1 (MQ=255)
              gTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGt                                                                                                                        <  2:6601/149‑1 (MQ=255)
              gTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGt                                                                                                                        <  4:162272/149‑1 (MQ=255)
              gTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGt                                                                                                                        <  1:83698/149‑1 (MQ=255)
                     ggAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAg                                                                                                                 <  5:256627/149‑1 (MQ=255)
                             cAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTg                                                                                                         >  5:49000/1‑149 (MQ=255)
                               aTGAGCAGCGTTGTAGCCATCGTCAAAAGGATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATGGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGa                                                                                                       <  3:172374/149‑1 (MQ=255)
                               aTGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTg                                                                                                         >  7:67765/1‑147 (MQ=255)
                                tGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAg                                                                                                      >  8:127000/1‑149 (MQ=255)
                                tGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAg                                                                                                      <  1:234166/149‑1 (MQ=255)
                                   gcagcGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAgctg                                                                                                   >  3:264713/1‑149 (MQ=255)
                                      gcGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAgctggct                                                                                                <  6:157817/149‑1 (MQ=255)
                                      gcGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAgctggct                                                                                                <  4:264713/149‑1 (MQ=255)
                                         ttGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGAt                                                                                             <  4:114371/149‑1 (MQ=255)
                                         ttGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGAt                                                                                             >  4:382370/1‑149 (MQ=255)
                                                cATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGAc                                                                                      >  7:312291/1‑149 (MQ=255)
                                                     tCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTc                                                                                 >  5:289904/1‑149 (MQ=255)
                                                       aaaaGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGTCCGACTGCTCGa                                                                               >  3:210630/1‑149 (MQ=255)
                                                       aaaaGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGa                                                                               >  8:165377/1‑149 (MQ=255)
                                                            cATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTg                                                                          <  5:160180/149‑1 (MQ=255)
                                                                   ggCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGa                                                                   <  1:350574/149‑1 (MQ=255)
                                                                    gCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGGTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGa                                                                   >  2:77407/1‑148 (MQ=255)
                                                                      cAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAg                                                                >  5:202301/1‑149 (MQ=255)
                                                                         gTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGcc                                                              >  4:183038/1‑148 (MQ=255)
                                                                             ttGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATg                                                          <  8:312291/148‑1 (MQ=255)
                                                                                     tCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGa                                                 <  6:49000/149‑1 (MQ=255)
                                                                                      ccccAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCGCCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGAt                                                >  5:237118/1‑149 (MQ=255)
                                                                                       cccAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATgg                                                         <  3:202698/139‑1 (MQ=255)
                                                                                       cccAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCCGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATgg                                                         >  4:202698/1‑139 (MQ=255)
                                                                                        ccAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATgc                                              >  2:33165/1‑149 (MQ=255)
                                                                                        ccAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATgc                                              >  7:334270/1‑149 (MQ=255)
                                                                                          agagGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATgcgc                                            <  7:127000/149‑1 (MQ=255)
                                                                                          agagGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATgcgc                                            <  6:202301/149‑1 (MQ=255)
                                                                                              ggCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGa                                        <  3:379554/149‑1 (MQ=255)
                                                                                                   ccGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAg                                     >  3:177703/1‑147 (MQ=255)
                                                                                                    cGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAg                                  <  6:289904/149‑1 (MQ=255)
                                                                                                      gCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCa                                >  4:295016/1‑149 (MQ=255)
                                                                                                        gtgtACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCAcc                              <  3:382370/149‑1 (MQ=255)
                                                                                                           tACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGAc                           >  3:365798/1‑149 (MQ=255)
                                                                                                                aTCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGa                       <  3:183038/148‑1 (MQ=255)
                                                                                                                aTCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAg                      >  4:290270/1‑149 (MQ=255)
                                                                                                                  cGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGcc                    >  4:254043/1‑149 (MQ=255)
                                                                                                                   gAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGc                     >  7:165285/1‑147 (MQ=255)
                                                                                                                    aaGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAg                  >  3:18139/1‑149 (MQ=255)
                                                                                                                     aGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAgg                 >  8:22708/1‑149 (MQ=255)
                                                                                                                     aGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAgg                 >  7:151602/1‑149 (MQ=255)
                                                                                                                       gTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAg                                                                <  1:192170/100‑1 (MQ=255)
                                                                                                                       gTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAg                                                                >  2:192170/1‑100 (MQ=255)
                                                                                                                          aCCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCa            <  8:334270/149‑1 (MQ=255)
                                                                                                                                 gTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAgg     >  7:360277/1‑149 (MQ=255)
                                                                                                                                   ggaccGGTGAGCACTTCCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAg           <  7:304534/138‑1 (MQ=255)
                                                                                                                                   ggaccGGTGAGCACTTCCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAg           >  8:304534/4‑141 (MQ=255)
                                                                                                                                     gCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCTCCGACCTGAGCCAGGCGGCAGCCCAGGCcg  <  1:33165/148‑1 (MQ=255)
                                                                                                                                      ccGGTGAGCACTTCCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAgc                                                      <  2:299660/95‑1 (MQ=255)
                                                                                                                                      ccGGTGAGCACTTCCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAgc                                                      >  1:299660/1‑95 (MQ=255)
                                                                                                                                      ccGGTGAGCACTTCCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGa                       <  6:560/126‑1 (MQ=255)
                                                                                                                                      ccGGTGAGCACTTCCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGa                       >  5:560/1‑126 (MQ=255)
                                                                                                                                                |                                                                                                                                        
AAGTACCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCCG  >  NC_002947/194464‑194744

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: ‑