Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 666,104 G→T 6.7% M990I (ATG→ATT) ‡ yfhM → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947666,1040GT 6.7% 188.2 / 4.0 60M990I (ATG→ATT) ‡yfhMlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base T (2/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

ACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTCGAAATCACCAC  >  NC_002947/665959‑666250
                                                                                                                                                 |                                                                                                                                                  
acacTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGcc                                                                                                                                                 <  6:336463/149‑1 (MQ=255)
   cTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGaa                                                                                                                                              <  2:312014/149‑1 (MQ=255)
    tACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAg                                                                                                                                             >  3:36650/1‑149 (MQ=255)
      cGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCCATGCCGAAgtg                                                                                                                                           <  2:275357/149‑3 (MQ=255)
       ggTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGc                                                                                                                                           >  3:201489/1‑148 (MQ=255)
       ggTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCa                                                                                                                                          <  1:130526/149‑1 (MQ=255)
              aTAGTCGCCCAGCAGCGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGt                                                                                                                                   >  1:219252/1‑149 (MQ=255)
              aTAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGt                                                                                                                                   >  2:108426/1‑149 (MQ=255)
                 gTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTgg                                                                                                                                 <  8:224859/148‑1 (MQ=255)
                    gCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTcgc                                                                                                                             >  4:194013/1‑149 (MQ=255)
                     cccAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTcgcc                                                                                                                            >  7:141942/1‑149 (MQ=255)
                      ccAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTcgccg                                                                                                                           >  7:276730/1‑149 (MQ=255)
                          cagAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTcgccgcgcc                                                                                                                       >  5:293756/1‑149 (MQ=255)
                          cagAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTcgccgcgcc                                                                                                                       <  3:194013/149‑1 (MQ=255)
                             agCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGtt                                                                                                                    <  2:206152/149‑1 (MQ=255)
                                   ccGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGAtt                                                                                                                <  8:141942/147‑1 (MQ=255)
                                                         gTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTc                                                                                        <  3:325919/149‑1 (MQ=255)
                                                          tGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTcc                                                                                       >  2:25210/1‑149 (MQ=255)
                                                              ggCCAGGCTACGGTCGACATCCCCGACTTCAAAGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACATCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCAGGc                                                                                   >  7:211617/1‑149 (MQ=255)
                                                               gCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGcc                                                                                  <  4:201489/149‑1 (MQ=255)
                                                               gCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGACGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGcc                                                                                  >  7:9755/1‑149 (MQ=255)
                                                                ccAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGc                                                                                                         <  7:249254/125‑1 (MQ=255)
                                                                ccAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGc                                                                                                         >  8:249254/1‑125 (MQ=255)
                                                                ccAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCg                                                                                 >  1:83442/1‑149 (MQ=255)
                                                                    gCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTgg                                                                             <  8:64722/149‑1 (MQ=255)
                                                                             gACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTAc                                                                    <  3:19464/149‑1 (MQ=255)
                                                                             gACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTAc                                                                    >  7:90804/1‑149 (MQ=255)
                                                                             gACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTAc                                                                    <  2:278484/149‑1 (MQ=255)
                                                                                   cccGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCt                                                              >  3:116692/1‑149 (MQ=255)
                                                                                    ccGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTg                                                             >  1:52555/1‑149 (MQ=255)
                                                                                    ccGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTg                                                             >  3:20695/1‑149 (MQ=255)
                                                                                    ccGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTg                                                             >  1:88014/1‑149 (MQ=255)
                                                                                    ccGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTg                                                             <  2:83442/149‑1 (MQ=255)
                                                                                     cGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCTGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTgg                                                            <  2:219252/149‑1 (MQ=255)
                                                                                         ttCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTg                                                                              <  6:86945/127‑1 (MQ=255)
                                                                                           cAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTg                                                       <  4:36650/148‑1 (MQ=255)
                                                                                               ggCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCt                                                  >  2:94593/1‑149 (MQ=255)
                                                                                               ggCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCt                                                  >  7:77667/1‑149 (MQ=255)
                                                                                                 cGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTgg                                                <  5:192289/149‑1 (MQ=255)
                                                                                                       gCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAAcc                                          >  7:213370/1‑149 (MQ=255)
                                                                                                        cGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCt                                         <  8:276730/149‑1 (MQ=255)
                                                                                                             gATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCg                                    >  4:186890/1‑149 (MQ=255)
                                                                                                                 gcgcAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCt                                         >  4:306545/1‑140 (MQ=255)
                                                                                                                 gcgcAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCt                                         <  3:306545/140‑1 (MQ=255)
                                                                                                                  cgcAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCcgc                               <  4:116692/149‑1 (MQ=255)
                                                                                                                         cTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAg                         >  8:270626/1‑148 (MQ=255)
                                                                                                                          catgccGAAGTGCTCTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGc                        >  2:140165/5‑148 (MQ=255)
                                                                                                                              ccGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACt                    >  3:173185/1‑148 (MQ=255)
                                                                                                                              ccGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACt                    >  2:67042/1‑148 (MQ=255)
                                                                                                                              ccGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACt                    >  1:83230/1‑148 (MQ=255)
                                                                                                                               cGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAAc                                           <  6:338994/124‑1 (MQ=255)
                                                                                                                               cGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAAc                                           >  5:338994/1‑124 (MQ=255)
                                                                                                                                 aaGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGc                <  4:173185/149‑1 (MQ=255)
                                                                                                                                    ttgccTTCGGCCATGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTc             >  2:203702/5‑149 (MQ=255)
                                                                                                                                    agCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTc             <  8:77667/149‑1 (MQ=255)
                                                                                                                                     tgccTTCGGCCATGCCGAAGTCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGt              <  4:231829/144‑1 (MQ=255)
                                                                                                                                     tgccTTCGGCCATGCCGAAGTCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGt              >  3:231829/4‑147 (MQ=255)
                                                                                                                                           cGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGGAACTGAGCGTCTAAAt        <  1:67042/147‑1 (MQ=255)
                                                                                                                                              cATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTCGAAATcacca   >  3:204996/1‑149 (MQ=255)
                                                                                                                                               aTGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTCGAAATcaccac  >  6:68109/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTCGAAATCACCAC  >  NC_002947/665959‑666250

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: