Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 931,487 C→T 7.5% R4712H (CGC→CAC)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947931,4870CT 7.5% 194.3 / 4.9 67R4712H (CGC→CAC) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base T (3/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCT  >  NC_002947/931341‑931634
                                                                                                                                                   |                                                                                                                                                   
cATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCgg                                                                                                                                                   >  5:173546/1‑149 (MQ=255)
cATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCgg                                                                                                                                                   >  2:108842/1‑149 (MQ=255)
 aTCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGc                                                                                                                                                  >  6:115685/1‑149 (MQ=255)
     gCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGt                                                                                                                                              <  6:173546/149‑1 (MQ=255)
       cTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGa                                                                                                                                            <  1:2045/149‑1 (MQ=255)
       cTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGa                                                                                                                                            <  7:276065/149‑1 (MQ=255)
          ggTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAggg                                                                                                                                         >  6:118510/1‑149 (MQ=255)
           gTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGt                                                                                                                                        <  4:181578/149‑1 (MQ=255)
            tGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTTAAGTGTCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTCTCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGt                                                                                                                                        <  1:96541/149‑1 (MQ=255)
            tGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCAATGGCATCGACGGTGGCGCTGCCGATTTCCTAGCCTGCCGGGTTGGTCAAGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGt                                                                                                                                        >  8:25705/1‑148 (MQ=255)
             gAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGATGGTgt                                                                                                                                      >  3:128868/1‑148 (MQ=255)
                gCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGGCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGa                                                                                                                                            >  1:281979/1‑140 (MQ=255)
                gCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGa                                                                                                                                            <  2:281979/140‑1 (MQ=255)
                 cTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGacac                                                                                                                                   >  1:69395/1‑148 (MQ=255)
                 cTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACAcc                                                                                                                                  <  3:324540/149‑1 (MQ=255)
                     gcagcTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGATGGTGTCACCCTCa                                                                                                                              <  4:128868/149‑1 (MQ=255)
                       agcTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCAcc                                                                                                                            <  3:13699/149‑1 (MQ=255)
                          tCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGTGGCCGGTGAGGGTGACACCATCACCgg                                                                                                                          >  6:268415/1‑148 (MQ=255)
                               gTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACgg                                                                                                                    <  5:118510/149‑1 (MQ=255)
                                ttGCCCCTTGCCGGGGTCAGCGTGTTGGTTATGTGGCCACTTGCATCGACGGTTGCGCTGCCGATTTCCTCGCCTGCCGGTTTTGTCACTGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGTTGTTGACACCATCACCGGAAACGGc                                                                                                                   <  7:232850/149‑1 (MQ=255)
                                ttGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGc                                                                                                                   <  3:107178/149‑1 (MQ=255)
                                    ccAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGTTTCCGGTGATGGTGTCACCCTCACCggccgcgg                                                                                                                    >  4:58278/1‑138 (MQ=38)
                                    ccAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGTTTCCGGTGATGGTGTCACCCTCACCggccgcgg                                                                                                                    <  3:58278/144‑7 (MQ=38)
                                      aGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAggg                                                                                                                                         <  8:59464/121‑1 (MQ=255)
                                      aGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAggg                                                                                                                                         >  7:59464/1‑121 (MQ=255)
                                           cGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAAc                                                                                                                      >  6:204703/1‑135 (MQ=255)
                                           cGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAAc                                                                                                                      <  5:204703/135‑1 (MQ=255)
                                                tCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGTGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGcc                                                                                                   >  4:223305/1‑149 (MQ=255)
                                                    cGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATg                                                                                               >  3:289795/1‑149 (MQ=255)
                                                     gTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATgg                                                                                              <  4:289795/149‑1 (MQ=255)
                                                     gTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATgg                                                                                              <  1:91144/149‑1 (MQ=255)
                                                      ttaTGTTGAAGTGGCCACTGGCATCGACGGTGGCGATGACGATTTCATCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGAGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATgg                                                                                              <  8:243517/146‑1 (MQ=255)
                                                        aTGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCgg                                                                                           >  5:106662/1‑149 (MQ=255)
                                                           gTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATc                                                                                                                               >  6:135648/1‑110 (MQ=255)
                                                           gTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATc                                                                                                                               <  5:135648/110‑1 (MQ=255)
                                                           gTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGgtgt                                                                                        >  1:144889/1‑149 (MQ=255)
                                                             gAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTgg                                                                                      >  5:101491/1‑149 (MQ=255)
                                                             gAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTgg                                                                                      >  8:169156/1‑149 (MQ=255)
                                                               aGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTa                                                                                    <  6:201658/149‑1 (MQ=255)
                                                                 tGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATc                                                                                  >  4:141476/1‑149 (MQ=255)
                                                                 tGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATc                                                                                  >  4:142099/1‑149 (MQ=255)
                                                                           cATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTg                                                                        <  7:252930/149‑1 (MQ=255)
                                                                                 cGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCAGCGGCCGGTGATGGTGACACCATCACCGGAAACGGCCAGTTCGCCTGGTGCCTATGGCGGTGTGGTATCGGGTTCGGTGAt                                                                      <  4:349190/145‑1 (MQ=255)
                                                                                      gcgcTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggc                                                             <  6:101491/149‑1 (MQ=255)
                                                                                      gcgcTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTgg                                                              <  1:204860/148‑1 (MQ=255)
                                                                                      gcgcTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGCGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggc                                                             <  4:140382/149‑1 (MQ=255)
                                                                                        gcTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggcgg                                                           <  1:108842/149‑1 (MQ=255)
                                                                                         cTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggcgg                                                           >  6:278175/1‑148 (MQ=255)
                                                                                          tGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGac                                                                                                                                     >  5:103481/1‑73 (MQ=255)
                                                                                          tGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGac                                                                                                                                     <  6:103481/73‑1 (MQ=255)
                                                                                          tGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggcgg                                                           >  6:376585/1‑147 (MQ=255)
                                                                                            ccGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGTGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCa                                                         <  3:223305/147‑1 (MQ=255)
                                                                                            ccGATGTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGATGGCGGCAGa                                                       >  7:336735/1‑149 (MQ=255)
                                                                                                    cTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGtt                                                 >  8:23021/1‑147 (MQ=255)
                                                                                                    cTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGtt                                                 <  7:23021/147‑1 (MQ=255)
                                                                                                    cTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTtgc                                               >  6:374677/1‑149 (MQ=255)
                                                                                                     tCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGAc                                                      >  4:16305/1‑141 (MQ=255)
                                                                                                     tCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGAc                                                      <  3:16305/141‑1 (MQ=255)
                                                                                                              cGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCAt                                      <  3:142099/149‑1 (MQ=255)
                                                                                                                      tCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAg                               >  5:5419/1‑148 (MQ=255)
                                                                                                                      tCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGc                              >  5:208771/1‑149 (MQ=255)
                                                                                                                         cGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGcacc                           <  5:376585/149‑1 (MQ=255)
                                                                                                                                   tGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCt                 <  7:25705/149‑1 (MQ=255)
                                                                                                                                    gcgcCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTg                >  7:175577/1‑149 (MQ=255)
                                                                                                                                     cgcCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTg                <  6:208771/148‑1 (MQ=255)
                                                                                                                                      gcCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgtt              <  6:5419/149‑1 (MQ=255)
                                                                                                                                          aCTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgttggcc          <  2:69395/149‑1 (MQ=255)
                                                                                                                                              cgccgcGGCCGGTGAGTGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgttggccgtt       <  4:261218/148‑1 (MQ=255)
                                                                                                                                               gccgcGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgttggccgttg      <  6:106662/148‑1 (MQ=255)
                                                                                                                                                  gcgGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCt  <  5:374677/149‑1 (MQ=255)
                                                                                                                                                  gcgGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCt  >  8:256004/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCT  >  NC_002947/931341‑931634

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: