Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 391,944 G→C 5.2% V740L (GTA→CTA)  soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947391,9440GC 5.2% 247.5 / 3.4 77V740L (GTA→CTA) soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base C (2/2);  total (41/36)
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.57e-01

ACCAACGCCTGGACCAAGCTCGACGTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTGAGCGAGCTGACCGACATCGACATGG  >  NC_002947/391797‑392091
                                                                                                                                                   |                                                                                                                                                   
acCAACGCCTGGACCAAGCTCGACGTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGt                                                                                                                                                    <  2:272454/149‑1 (MQ=255)
 ccAACGCCTGGACCAAGCTCGACGTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTa                                                                                                                                                   >  3:13268/1‑149 (MQ=255)
   aaCGCCTGGACCAAGCTCGACGTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTAtt                                                                                                                                                 >  6:396426/1‑149 (MQ=255)
        cTGTACCATGCGCGACGTGGTCATGGACGGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAgc                                                                                                                                            <  8:139753/149‑2 (MQ=255)
        cTGGACCAAGCTCGACGTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGt                                                                                                                                            >  6:147252/1‑149 (MQ=255)
         tGGACCAAGCTCGACGTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGtg                                                                                                                                           >  6:93356/1‑149 (MQ=255)
              cAAGCTCGACGTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGCTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACATCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCTGtggatg                                                                                                                                       >  2:372687/1‑148 (MQ=255)
                aGCTCGACGTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAg                                                                                                                                    >  1:222392/1‑149 (MQ=255)
                  cTCGACGTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCt                                                                                                                                  <  5:396426/149‑1 (MQ=255)
                        gTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTATATCa                                                                                                                            <  1:291218/149‑1 (MQ=255)
                           ggCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGCCCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCAc                                                                                                                           >  6:235993/1‑147 (MQ=255)
                           ggCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCa                                                                                                                         >  4:375874/1‑149 (MQ=255)
                           ggCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCa                                                                                                                         <  4:55060/149‑1 (MQ=255)
                            gCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCAc                                                                                                                           >  1:327625/1‑146 (MQ=255)
                                ggCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGAcc                                                                                                                     >  2:212105/1‑148 (MQ=255)
                                   ccGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAAt                                                                                                                 >  6:202998/1‑149 (MQ=255)
                                    cGCTACGGCGTGTTGTGCAAGGAAGACGGCTTGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACGACCACCAGCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATg                                                                                                                <  2:165023/149‑1 (MQ=255)
                                      cTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGc                                                                                                              <  5:147252/149‑1 (MQ=255)
                                           gCCTTATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCagag                                                                                                         <  2:363990/149‑1 (MQ=255)
                                           gCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCagag                                                                                                         <  3:353653/149‑1 (MQ=255)
                                            ccTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGc                                                                                                              >  6:231164/1‑143 (MQ=255)
                                            ccTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGc                                                                                                              <  5:231164/143‑1 (MQ=255)
                                             cTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCt                                                                                                       >  3:395418/1‑149 (MQ=255)
                                                aTGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGaa                                                                                                    <  5:166798/149‑1 (MQ=255)
                                                 tgtgCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGaa                                                                                                    <  4:302925/148‑1 (MQ=255)
                                                    gCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGgtg                                                                                                >  5:107014/1‑149 (MQ=255)
                                                     cAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGgtg                                                                                                <  2:222392/148‑1 (MQ=255)
                                                     cAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCAACGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGgtg                                                                                                >  4:50711/1‑148 (MQ=255)
                                                       aGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGACTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTac                                                                                             <  7:301286/149‑1 (MQ=255)
                                                         gAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTactt                                                                                           <  2:288505/149‑1 (MQ=255)
                                                            gACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCACACCGAATGGCCAGAGCTGAAGGTGTacttaac                                                                                        >  6:310316/1‑149 (MQ=255)
                                                             aCGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTacttcact                                                                                       >  6:183019/1‑149 (MQ=255)
                                                              cGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTacttcactt                                                                                      <  1:212105/149‑1 (MQ=255)
                                                                gCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGc                                                                                                              <  6:245582/123‑1 (MQ=255)
                                                                gCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGc                                                                                                              >  5:245582/1‑123 (MQ=255)
                                                                 cATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCgg                                                                                   >  6:39581/1‑149 (MQ=255)
                                                                     gTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCa                                                                                <  3:342403/148‑1 (MQ=255)
                                                                     gTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCAc                                                                               <  7:219632/149‑1 (MQ=255)
                                                                      tCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCAcc                                                                              >  8:68513/1‑149 (MQ=255)
                                                                      tCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCAcc                                                                              >  5:9195/1‑149 (MQ=255)
                                                                      tCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCAcc                                                                              <  2:356250/149‑1 (MQ=255)
                                                                       cTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCg                                                                             >  3:55035/1‑149 (MQ=255)
                                                                              gacgGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTg                                                                      <  1:178294/149‑1 (MQ=255)
                                                                              gacgGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTg                                                                      <  5:2896/149‑1 (MQ=255)
                                                                              gacgGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTg                                                                      >  3:245552/1‑149 (MQ=255)
                                                                              gacgGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTg                                                                      >  2:396165/1‑149 (MQ=255)
                                                                                         cGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGAGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACAGACTGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGa                                                           <  5:310316/149‑1 (MQ=255)
                                                                                          gCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGAc                                                          >  8:116157/1‑149 (MQ=255)
                                                                                           ccTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTac                                                                                             <  4:302580/113‑1 (MQ=255)
                                                                                           ccTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGAcc                                                         <  2:327625/149‑1 (MQ=255)
                                                                                            cTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGAccc                                                        >  6:128296/1‑149 (MQ=255)
                                                                                                gTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTc                                                    <  3:196659/149‑1 (MQ=255)
                                                                                                   ggCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCgg                                                 <  6:9195/149‑1 (MQ=255)
                                                                                                    gCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCGGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGc                                                <  8:32042/149‑1 (MQ=255)
                                                                                                      gACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGccc                                              >  6:238448/1‑149 (MQ=255)
                                                                                                      gACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGccc                                              >  5:322760/1‑149 (MQ=255)
                                                                                                         aaCCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTaa                                           >  8:405901/1‑149 (MQ=255)
                                                                                                           ccACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTaa                                           >  5:199569/1‑147 (MQ=255)
                                                                                                            cACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAg                                        >  3:34368/1‑149 (MQ=255)
                                                                                                                       gACCACCACCACCGGCGGCGCTGCCCGCCTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCa                                   <  5:270540/143‑1 (MQ=255)
                                                                                                                       gACCACCACCACCGGCGGCGCTGCCCGCCTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCa                                   >  6:276823/1‑143 (MQ=255)
                                                                                                                       gACCACCACCACCGGCGGCGCTGCCCGCCTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCa                                   <  5:276823/143‑1 (MQ=255)
                                                                                                                       gACCACCACCACCGGCGGCGCTGCCCGCCTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCa                                   >  6:270540/1‑143 (MQ=255)
                                                                                                                         ccaccaccaccGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCAGCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAgctgct                            >  1:112735/1‑148 (MQ=255)
                                                                                                                         ccaccaccaccGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAgctgct                            >  6:38595/1‑148 (MQ=255)
                                                                                                                         ccaccaccaccGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAgctgct                            >  8:253659/1‑148 (MQ=255)
                                                                                                                          caccaccaccGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTga                          <  4:13268/149‑1 (MQ=255)
                                                                                                                            ccaccaccGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTgagc                        >  7:252052/1‑149 (MQ=255)
                                                                                                                            ccaccaccGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTgagc                        <  5:202998/149‑1 (MQ=255)
                                                                                                                             caccaccGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTgagc                        >  4:103723/1‑148 (MQ=255)
                                                                                                                             caccaccGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTgagc                        >  1:274274/1‑148 (MQ=255)
                                                                                                                               ccaccGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTgagcgag                     <  3:103723/149‑1 (MQ=255)
                                                                                                                                 accGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTGAGCGAGCt                   >  4:299156/1‑149 (MQ=255)
                                                                                                                                        gcgcTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTGAGCGAGCTGACcgac            <  5:93356/149‑1 (MQ=255)
                                                                                                                                          gcTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTGAGCGAGCTGACcgacat          <  5:38595/149‑1 (MQ=255)
                                                                                                                                            tGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGcc                                               >  4:317469/1‑110 (MQ=255)
                                                                                                                                            tGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGcc                                               <  3:317469/110‑1 (MQ=255)
                                                                                                                                                  cgTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTGAGCGAGCTGACCGACATCGACATgg  <  7:116157/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCAACGCCTGGACCAAGCTCGACGTGGGCAAGGCCCGCTACGGCCTGATGTGCAAGGAAGACGGCATGGTCTTCGACGACGGCGTAACCGCCTGCGTCGGCGACAACCACTTCATCATGACCACCACCACCGGCGGCGCTGCCCGCGTATTGCAGTGGATGGAGCTGTATCACCAGACCGAATGGCCAGAGCTGAAGGTGTACTTCACTTCGGTCACCGACCACTGGGCCACCATGACCCTGTCCGGCCCTAACAGCCGCAAGCTGCTGAGCGAGCTGACCGACATCGACATGG  >  NC_002947/391797‑392091

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