Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 194,608 A→C 6.4% 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,6080AC 6.4% 278.0 / 6.9 89E38D (GAA→GAC) ‡PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (44/39);  new base C (2/4);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GGGAAGTACCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCCGCGCCTTCGGT  >  NC_002947/194461‑194754
                                                                                                                                                   |                                                                                                                                                  
gggAAGTACCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAg                                                                                                                                                   <  4:185647/149‑1 (MQ=255)
gggAAGTACCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAg                                                                                                                                                   >  6:53351/1‑149 (MQ=255)
   aaGTACCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCAGCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGc                                                                                                                                                >  5:134150/1‑149 (MQ=255)
   aaGTACCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGGTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGATGGTGACCGCCTGTTGGCCTGTGAGTAAGTGc                                                                                                                                                <  6:436021/149‑1 (MQ=255)
        ccccATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCTTTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGGTGGCCGGTGAGGAAGTGCTCAc                                                                                                                                            <  2:270973/148‑1 (MQ=255)
        ccccATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCAc                                                                                                                                            <  5:36917/148‑1 (MQ=255)
         cccATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAATTGCTCACCg                                                                                                                                          <  7:245635/149‑1 (MQ=255)
         cccATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGCGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCAGGTGAGGAAGTGCTCAcct                                                                                                                                          >  3:311596/1‑148 (MQ=255)
          ccATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGCACTTCCTCACCgg                                                                                                                                         <  7:234528/149‑1 (MQ=255)
          ccATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGCACTTCCTCACCgg                                                                                                                                         <  7:240841/149‑1 (MQ=255)
            aTCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCg                                                                                                                                          <  8:24901/146‑1 (MQ=255)
            aTCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCg                                                                                                                                          >  7:24901/1‑146 (MQ=255)
            aTCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTc                                                                                                                                       >  5:84560/1‑149 (MQ=255)
                ggTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCggg                                                                                                                                   >  4:307130/1‑149 (MQ=255)
                  tCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTg                                                                                                                                 <  6:224121/149‑1 (MQ=255)
                   cacaCGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTgg                                                                                                                                >  5:167332/1‑149 (MQ=255)
                    acacGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGgc                                                                                                                               >  1:122595/1‑149 (MQ=255)
                       cGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGgcg                                                                                                                              >  5:199416/1‑147 (MQ=255)
                        ggAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAg                                                                                                                           <  7:15342/149‑1 (MQ=255)
                               ccAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACttt                                                                                                                    >  2:14556/1‑149 (MQ=255)
                                cAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTg                                                                                                                   <  8:183520/149‑1 (MQ=255)
                                      gcagcGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAgct                                                                                                              <  7:37934/148‑1 (MQ=255)
                                       ccgcGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAgctgg                                                                                                            <  7:401819/147‑1 (MQ=255)
                                       cagcGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAgctgg                                                                                                            >  6:308532/1‑149 (MQ=255)
                                        agcGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAgctggc                                                                                                           <  1:230618/149‑1 (MQ=255)
                                          cGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAgctggctg                                                                                                         <  4:143903/149‑1 (MQ=255)
                                            ttGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGAt                                                                                                       >  5:368749/1‑149 (MQ=255)
                                              gTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATg                                                                                                      <  4:373773/148‑1 (MQ=255)
                                               tAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATgg                                                                                                     <  6:336304/148‑1 (MQ=255)
                                               tAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGc                                                                                                    >  3:115010/1‑149 (MQ=255)
                                                          aaaaGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGa                                                                                         <  6:368749/149‑1 (MQ=255)
                                                          aaaaGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGa                                                                                         <  6:199416/149‑1 (MQ=255)
                                                           aaaGCATTGTCGCCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGAt                                                                                        >  6:326141/1‑149 (MQ=255)
                                                             aGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTAc                                                                                                                       <  6:254148/116‑1 (MQ=255)
                                                             aGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTAc                                                                                                                       >  5:254148/1‑116 (MQ=255)
                                                              gCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCtt                                                                                     >  1:235689/1‑149 (MQ=255)
                                                               cATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGTCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTg                                                                                    <  7:201393/149‑1 (MQ=255)
                                                                aTTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTg                                                                                    >  1:48120/1‑148 (MQ=255)
                                                                aTTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACAGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTg                                                                                    <  4:268675/148‑1 (MQ=255)
                                                                 ttGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGa                                                                                  <  2:122595/149‑1 (MQ=255)
                                                                    tCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGt                                                                               >  6:71415/1‑149 (MQ=255)
                                                                    tCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGt                                                                               <  5:28553/149‑1 (MQ=255)
                                                                      ggCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGa                                                                             >  4:161964/1‑149 (MQ=255)
                                                                      ggCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGa                                                                             >  2:373396/1‑149 (MQ=255)
                                                                      ggCCAGGGTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTg                                                                              <  6:134150/148‑1 (MQ=255)
                                                                         cAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAg                                                                          >  1:296101/1‑149 (MQ=255)
                                                                         cAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAg                                                                          >  6:71129/1‑149 (MQ=255)
                                                                            gTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACGTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGcc                                                                        >  2:288463/1‑148 (MQ=255)
                                                                                   cGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAg                                                                 <  1:14556/148‑1 (MQ=255)
                                                                                       aTCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGcc                                                             >  4:242668/1‑148 (MQ=255)
                                                                                        tCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGGTGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGa                                                           >  6:296635/1‑149 (MQ=255)
                                                                                        tCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGa                                                           <  2:235689/149‑1 (MQ=255)
                                                                                        tCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGa                                                           >  3:100504/1‑149 (MQ=255)
                                                                                          cccAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGTTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGAt                                                          <  3:307130/148‑1 (MQ=255)
                                                                                           ccAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATgc                                                        >  5:367562/1‑149 (MQ=255)
                                                                                           ccAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATgc                                                        >  3:244862/1‑149 (MQ=255)
                                                                                             agagGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATgcgc                                                      <  3:161964/149‑1 (MQ=255)
                                                                                               agGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCt                                                    <  6:167332/149‑1 (MQ=255)
                                                                                                 ggCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGTGCCTGa                                                  >  4:360899/1‑149 (MQ=255)
                                                                                                 ggCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGa                                                  >  8:63760/1‑149 (MQ=255)
                                                                                                 ggCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGACCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGa                                                  >  8:66399/1‑149 (MQ=255)
                                                                                                  gCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGAc                                                 <  3:242668/149‑1 (MQ=255)
                                                                                                   cATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACa                                                >  6:129802/1‑149 (MQ=255)
                                                                                                        ggCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAg                                            >  6:406853/1‑148 (MQ=255)
                                                                                                           gtgtACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCAcc                                        >  7:293521/1‑149 (MQ=255)
                                                                                                              tACCTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGAc                                     <  1:288463/149‑1 (MQ=255)
                                                                                                                cTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGGTGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCt                                   <  5:296635/149‑1 (MQ=255)
                                                                                                                   aTCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGa                                 <  4:100504/148‑1 (MQ=255)
                                                                                                                    tCGAAGGTGACCGCCTGTTGGCCGGTGAGGACGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGc                               <  4:218260/149‑1 (MQ=255)
                                                                                                                      gAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGc                                           <  2:35958/135‑1 (MQ=255)
                                                                                                                      gAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGc                               >  3:128989/1‑147 (MQ=255)
                                                                                                                      gAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCAAAGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGc                                           >  1:35958/1‑135 (MQ=255)
                                                                                                                            gACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAggcgg                        >  6:203738/1‑148 (MQ=255)
                                                                                                                             aCCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAg                                               >  8:153805/1‑124 (MQ=255)
                                                                                                                             aCCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAg                                               <  7:153805/124‑1 (MQ=255)
                                                                                                                              ccGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAg                     <  7:63760/149‑1 (MQ=255)
                                                                                                                                     cggaccGGTGAGCACTTCCTCACCGGCCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGc              >  2:433684/5‑149 (MQ=255)
                                                                                                                                      tGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGACCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGcc             <  7:66399/149‑1 (MQ=255)
                                                                                                                                       ggCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCcg            >  5:77796/1‑149 (MQ=255)
                                                                                                                                       ggCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCcg            >  4:112325/1‑149 (MQ=255)
                                                                                                                                       ggCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCcg            >  5:23690/1‑149 (MQ=255)
                                                                                                                                        accGGTGAGCACTTCCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCcgc           >  7:22453/2‑149 (MQ=255)
                                                                                                                                         ccGGTGAGCACTTCCTCACCGGCCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCcgcg          <  1:433684/149‑1 (MQ=255)
                                                                                                                                          cGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCCgcgc         >  5:76630/1‑149 (MQ=255)
                                                                                                                                          cGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCCgcgc         <  6:367562/149‑1 (MQ=255)
                                                                                                                                             tGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGACAATGGAGCGCAGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCCGCGCCtt      >  2:407273/1‑149 (MQ=255)
                                                                                                                                                ggtaGTGCTCACCGGTCCGTGTGGGGCAGTTACTTTGGAGCTGGCTTATGGCCGACGGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGTCAGGCGGCAGCCCAGGCCGCGCCTTCgg   <  7:71368/146‑1 (MQ=255)
                                                                                                                                                ggAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCCGCGCCTTCgg   >  1:319260/1‑149 (MQ=255)
                                                                                                                                                ggAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCCGCGCCTTCgg   <  6:46418/149‑1 (MQ=255)
                                                                                                                                                 gAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCCGCGCCTTCGGt  <  4:115010/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGGAAGTACCCCATCGGGTCACACGGAGAGTCCAATGAGCAGCGTTGTAGCCATCGTCAAAAGCATTGTCGGCCAGGTTATTGCGGTATCCCCAGAGGGCATCCGGCGTGTACTTATCGAAGGTGACCGCCTGTTGGCCGGTGAGGAAGTGCTCACCGGTCCGGGTGGCGCAGTTACTTTGGAGCTGGCTGATGGCCGACTGCTCGATCTTGGACGTGACAGCCAATGGAGCGCCGATGCGCCTGACAGTAGCACCGACCTGAGCCAGGCGGCAGCCCAGGCCGCGCCTTCGGT  >  NC_002947/194461‑194754

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