Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,149,241 A→C 5.1% K150Q (AAG→CAG)  PP_2765 → sulfonate monooxygenase MsuD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,149,2410AC 5.1% 254.4 / 4.1 79K150Q (AAG→CAG) PP_2765sulfonate monooxygenase MsuD
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base C (2/2);  total (40/39)
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.01e-01

TTGGCAAGCATCGACCAGTTCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCTCGCGGGCCGCGCGCGACATGGCCTC  >  NC_002947/3149094‑3149386
                                                                                                                                                    |                                                                                                                                                 
ttGGCAAGCATCGACCAGTTCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGaa                                                                                                                                                  <  1:72934/149‑1 (MQ=255)
  ggCAAGCATCGACCAGTTCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAg                                                                                                                                                 <  8:14266/148‑1 (MQ=255)
     aaGCATCGACCAGTTCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCa                                                                                                                                             >  2:13322/1‑149 (MQ=255)
     aaGCATCGACCAGTTCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCa                                                                                                                                             >  7:302077/1‑149 (MQ=255)
      aGCATCGACCAGTTCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCAt                                                                                                                                            >  4:47924/1‑149 (MQ=255)
        cATCGACCAGTTCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCt                                                                                                                                          <  4:177893/149‑1 (MQ=255)
            gACCAGTTCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGa                                                                                                                                       <  4:131359/148‑1 (MQ=255)
                 gTTCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAg                                                                                                                                  >  5:261491/1‑148 (MQ=255)
                  ttCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGa                                                                                                                                <  6:261491/149‑1 (MQ=255)
                     aCCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACaa                                                                                                                             >  1:370625/1‑149 (MQ=255)
                        gcgGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAAttt                                                                                                                          <  7:345570/149‑1 (MQ=255)
                           ggCCGCATCTCCATCAACGTGGTGTGCGGCTGGTTCAAGGGCGCGTTCCGCGCCATCTGAGAGACATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTTCCCTGAAGGGCATCTGGACCCAGGACAATTTCAg                                                                                                                       <  6:313385/149‑1 (MQ=255)
                                 aTCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTcaa                                                                                                                 >  5:73184/1‑147 (MQ=255)
                                 aTCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCa                                                                                                                 >  5:256523/1‑149 (MQ=255)
                                    gCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACgg                                                                                                              <  5:335995/149‑1 (MQ=255)
                                         cAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGAtt                                                                                                         >  4:35953/1‑149 (MQ=255)
                                            cGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCt                                                                                                      <  2:370625/149‑1 (MQ=255)
                                             gtggtgTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTa                                                                                                     <  1:13322/149‑1 (MQ=255)
                                             gtggtgTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTa                                                                                                     <  2:160485/149‑1 (MQ=255)
                                               ggtgTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTccc                                                                                                   >  2:47066/1‑149 (MQ=255)
                                                gtgtCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCg                                                                                                  >  8:219808/1‑149 (MQ=255)
                                                   tCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGttt                                                                                               <  7:219808/149‑1 (MQ=255)
                                                     cGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTc                                                                                              <  5:320575/148‑1 (MQ=255)
                                                      ggCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCaa                                                                                            >  6:150233/1‑149 (MQ=255)
                                                        cTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGCCGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAAc                                                                                           <  1:303512/149‑1 (MQ=255)
                                                        cTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCa                                                                                             >  2:207721/1‑146 (MQ=255)
                                                           gTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGAc                                                                                        >  1:301528/1‑148 (MQ=255)
                                                             tCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGACTGGACCATGATGAGCGCTATCGCCCATCCGAGGAGTTCATTCGTGCACTGAAGAGAATCTGGACCCAGGACAATTTCAGCTTCAAAGGCGATTTCTACCGTTTCAACGAGTa                                                                                      >  1:248840/1‑148 (MQ=255)
                                                              cAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTaca                                                                                    >  8:210868/1‑149 (MQ=255)
                                                                 gggCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACAcc                                                                                  >  4:338262/1‑148 (MQ=255)
                                                                       gTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTcaagc                                                                           >  1:44437/1‑149 (MQ=255)
                                                                        ttCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTACCACGGCGATTTCTACCGTTTCAACGACTACACCCTcaagcc                                                                          >  8:101763/1‑149 (MQ=255)
                                                                                cATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGAt                                                                                                          >  2:44331/1‑109 (MQ=255)
                                                                                cATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGAt                                                                                                          <  1:44331/109‑1 (MQ=255)
                                                                                cATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTcaagcc                                                                          <  8:50810/141‑1 (MQ=255)
                                                                                cATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTcaagcc                                                                          >  7:50810/1‑141 (MQ=255)
                                                                                  tCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCa                                                                                                                        >  1:220773/1‑93 (MQ=255)
                                                                                  tCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCa                                                                                                                        <  2:220773/93‑1 (MQ=255)
                                                                                    ggCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTgc                                                               <  1:172883/148‑1 (MQ=255)
                                                                                     gCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTgcag                                                             <  3:338262/149‑1 (MQ=255)
                                                                                       gAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCaa                                                          >  3:148016/1‑149 (MQ=255)
                                                                                            aTGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCg                                                       >  7:142621/1‑148 (MQ=255)
                                                                                             tGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCTATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCa                                                     <  5:150233/149‑1 (MQ=255)
                                                                                              ggCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTgca                                                              >  2:107335/1‑139 (MQ=255)
                                                                                              ggCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTgca                                                              <  1:107335/139‑1 (MQ=255)
                                                                                               gCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCAc                                                    <  1:252242/148‑1 (MQ=255)
                                                                                                    aCCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGa                                               >  6:219569/1‑148 (MQ=255)
                                                                                                    aCCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGa                                               >  1:2740/1‑148 (MQ=255)
                                                                                                          atgaGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCtt                                         <  7:210868/148‑1 (MQ=255)
                                                                                                           tgaGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCCACCGCACCCGGAAATcgtg                                        >  4:118092/1‑145 (MQ=255)
                                                                                                            gaGCGCTATCGCCGATCCGAGGAGTACAATCGTGCCCTGAAGCGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGAATTCTACCGATTCAACGACTACACCCACAAGCCCAAGCCGCTGCAGTAACCGCACCCGGAAATCTTc                                        >  3:244884/1‑147 (MQ=255)
                                                                                                                    tCGCCGATCCGAGGAGTTCATTCGTGCCCTTTAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCgg                                <  2:248840/147‑1 (MQ=255)
                                                                                                                    tCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACGGCACCCGGAAATCTTCCAGGGCgg                                <  1:207721/147‑1 (MQ=255)
                                                                                                                    tCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCgg                                <  6:256523/147‑1 (MQ=255)
                                                                                                                    tCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGa                              <  8:142621/149‑1 (MQ=255)
                                                                                                                      gCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAg                             >  2:366133/1‑148 (MQ=255)
                                                                                                                      gCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGc                            >  5:193383/1‑149 (MQ=255)
                                                                                                                        cGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCTc                          <  2:2740/149‑1 (MQ=255)
                                                                                                                              gaggagTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAggg                                   <  2:319732/134‑1 (MQ=255)
                                                                                                                              gaggagTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCTCGCGGGc                    >  7:229473/1‑149 (MQ=255)
                                                                                                                               aggagTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCTCGCGGGcc                   <  3:47924/149‑1 (MQ=255)
                                                                                                                                ggagTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCACGCGGGCcg                  >  2:262924/1‑149 (MQ=255)
                                                                                                                                   gTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCTCGCGGGCcgcgc               <  5:219569/149‑1 (MQ=255)
                                                                                                                                    ttCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCgg                                <  5:375929/131‑1 (MQ=255)
                                                                                                                                    ttCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCgg                                >  6:375929/1‑131 (MQ=255)
                                                                                                                                      cATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCTCGCGGGCcgcgcgc             <  1:47066/148‑1 (MQ=255)
                                                                                                                                         cagatGCCCTTCAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGa                              <  5:390192/125‑1 (MQ=255)
                                                                                                                                         cagatGCCCTTCAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGa                              >  6:390192/5‑129 (MQ=255)
                                                                                                                                          cGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCTCGCGGGCCGCGCGCGACa         <  3:35953/148‑1 (MQ=255)
                                                                                                                                           gatGCCCTTCAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGa                              <  2:247999/125‑1 (MQ=255)
                                                                                                                                           gatGCCCTTCAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGa                              >  1:247999/3‑127 (MQ=255)
                                                                                                                                           gTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCt                                          <  6:170411/114‑1 (MQ=255)
                                                                                                                                           gTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCt                                          >  5:170411/1‑114 (MQ=255)
                                                                                                                                           gTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCTCGCGGGCCGCGCGCGACATg       >  2:370048/1‑149 (MQ=255)
                                                                                                                                           cctGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTGCAACGACTACACCCCCACGCCCAAGCCGCTGCAGCACCCGCACCCGGACATCTTCCAGGGCGGGAGCTCGCGGGCCGCGCGCGACa         <  5:359627/146‑1 (MQ=255)
                                                                                                                                               cccTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCt                                          >  4:368553/1‑111 (MQ=255)
                                                                                                                                               cccTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCt                                          <  3:368553/111‑1 (MQ=255)
                                                                                                                                                ccTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAATCCCAAGCCGCt                                                                 <  6:171158/87‑1 (MQ=255)
                                                                                                                                                ccTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCt                                                                 >  5:171158/1‑87 (MQ=255)
                                                                                                                                                 cTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCTCGGGGGCCGCGCGCGACATGGCCTc  >  5:340378/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTGGCAAGCATCGACCAGTTCACCGCGGGCCGCATCGCCATCAACGTGGTGTCCGGCTGGTTCAAGGGCGAGTTCCGCGCCATCGGCGAGCCATGGCTGGACCATGATGAGCGCTATCGCCGATCCGAGGAGTTCATTCGTGCCCTGAAGGGCATCTGGACCCAGGACAATTTCAGCTTCCACGGCGATTTCTACCGTTTCAACGACTACACCCTCAAGCCCAAGCCGCTGCAGCAACCGCACCCGGAAATCTTCCAGGGCGGGAGCTCGCGGGCCGCGCGCGACATGGCCTC  >  NC_002947/3149094‑3149386

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: