Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,135,131 C→A 5.1% K540N (AAG→AAT fadE ← medium‑long chain acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,135,1310CA 5.1% 254.8 / 3.1 78K540N (AAG→AATfadEmedium‑long chain acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (33/41);  new base A (2/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

GCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGTGATGAAGATCGGCG  >  NC_002947/2134987‑2135273
                                                                                                                                                |                                                                                                                                              
gCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCtt                                                                                                                                              <  7:211569/147‑1 (MQ=255)
gCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCtt                                                                                                                                              <  6:361953/147‑1 (MQ=255)
gCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGa                                                                                                                                            <  3:402456/149‑1 (MQ=255)
 cGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAg                                                                                                                                           >  5:46187/1‑149 (MQ=255)
      cAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGcc                                                                                                                                       <  5:367397/148‑1 (MQ=255)
        gggCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGa                                                                                                                                    >  5:185075/1‑149 (MQ=255)
          gcgcGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATc                                                                                                                                  <  7:316834/149‑1 (MQ=255)
              ggAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTgg                                                                                                                              <  5:200933/149‑1 (MQ=255)
                   tAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTc                                                                                                                         <  8:197285/149‑1 (MQ=255)
                         tGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGccc                                                                                                                   >  1:49282/1‑149 (MQ=255)
                             tCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGc                                                                                                                >  1:99271/1‑148 (MQ=255)
                                gggCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAAt                                                                                                            <  4:29825/149‑1 (MQ=255)
                                  gCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGACGTAGACCAAGGTACTGGCCGCAGTACCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGAACCTGGTCTTCACGCCCGGCCAATGc                                                                                                          >  5:287441/1‑149 (MQ=255)
                                   cATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGcc                                                                                                         <  3:260824/149‑1 (MQ=255)
                                    aTCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCa                                                                                                        <  6:327279/149‑1 (MQ=255)
                                     tCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCAt                                                                                                       <  5:295670/149‑1 (MQ=255)
                                       gCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCAtt                                                                                                      <  6:185075/148‑1 (MQ=255)
                                        cccGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTc                                                                                                    >  2:183770/1‑149 (MQ=255)
                                        cccGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTc                                                                                                    <  5:212756/149‑1 (MQ=255)
                                         ccGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTcc                                                                                                   <  6:46187/149‑1 (MQ=255)
                                             cACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCTAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGa                                                                                               <  1:45614/149‑1 (MQ=255)
                                                ccGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCGGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCa                                                                                            >  5:78543/1‑149 (MQ=255)
                                                ccGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCa                                                                                            >  6:194530/1‑149 (MQ=255)
                                                    tCAAAGCGACCCAGCCCCGGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATAAGCAGGTCATCGAATTACTTGAGTGCCTGATCGTGGTCTTCACGACCGGCCAATGCCATTTCCTTGAGCACGaa                                                                                        >  5:22322/1‑149 (MQ=255)
                                                     cAAAGCGACCCAGCCCCAGGCCGTAGACCANGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCATTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAAc                                                                                       >  1:32978/1‑149 (MQ=255)
                                                     cAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCATTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAAc                                                                                       >  1:38981/1‑149 (MQ=255)
                                                        aGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCTGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGa                                                                                    <  6:78543/149‑1 (MQ=255)
                                                        aGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGa                                                                                    <  2:288844/149‑1 (MQ=255)
                                                         gCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGAt                                                                                   >  4:381000/1‑149 (MQ=255)
                                                         gCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGAt                                                                                   <  7:79288/149‑1 (MQ=255)
                                                          cGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATg                                                                                  <  1:114557/149‑1 (MQ=255)
                                                            aCCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTTATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGc                                                                                <  2:49282/149‑1 (MQ=255)
                                                             cccAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCa                                                                               <  3:381000/149‑1 (MQ=255)
                                                              ccAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAg                                                                              <  2:12200/149‑1 (MQ=255)
                                                                aGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGc                                                                             >  6:185011/1‑148 (MQ=255)
                                                                   cccAGGCCGTAGACCAAGGTACTGGCCGGATTGCCGTGGTCGAAAATGATGTGTTTCATCAGCAGGTCATCGTATTCCTTGTGTGCCTGATCGTTGTATTGACGCCCGGGCAATGCCATTTCCTTGAGCAGGAACGGATGGCAGCGAAt                                                                         <  3:182899/149‑1 (MQ=255)
                                                                    ccAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCTGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATg                                                                        <  5:285141/149‑1 (MQ=255)
                                                                       ggCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGgcg                                                                     <  2:99271/149‑1 (MQ=255)
                                                                        gCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGgcgc                                                                    >  3:155242/1‑149 (MQ=255)
                                                                          cGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGccc                                                                  >  6:311937/1‑149 (MQ=255)
                                                                           gTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCt                                                                 <  5:194530/149‑1 (MQ=255)
                                                                            tAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTg                                                                <  4:155242/149‑1 (MQ=255)
                                                                               aCCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGcc                                                             >  4:300578/1‑149 (MQ=255)
                                                                                  aaGGTACTGTCCTCATTTCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGGCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGaa                                                          <  3:300369/149‑1 (MQ=255)
                                                                                         tGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGatcat                                                    >  2:40894/1‑148 (MQ=255)
                                                                                          ggCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGTCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGatcat                                                    >  1:41052/1‑147 (MQ=255)
                                                                                           gCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAg                                                 >  2:294859/1‑149 (MQ=255)
                                                                                           gCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAg                                                 <  3:160054/149‑1 (MQ=255)
                                                                                            ccGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAgg                                                >  2:118957/1‑149 (MQ=255)
                                                                                                   gCCGGCGGCTAACAGGATGTTTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGGCCGGCCAATGCCATTTCCTTGAGGACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTTCGt                                          <  6:22322/148‑1 (MQ=255)
                                                                                                   gCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCATTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGt                                          <  2:38981/148‑1 (MQ=255)
                                                                                                   gCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCATTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGt                                          <  2:32978/148‑1 (MQ=255)
                                                                                                             aaCAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGATGCGTGAAAGGATgat                               >  6:195585/1‑147 (MQ=255)
                                                                                                              aCAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTg                              <  1:118957/149‑1 (MQ=255)
                                                                                                              aCAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTg                              >  3:198183/1‑149 (MQ=255)
                                                                                                                 ggATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGa                                                        <  1:136571/120‑1 (MQ=255)
                                                                                                                 ggATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGa                                                        >  2:136571/1‑120 (MQ=255)
                                                                                                                 ggATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCa                           <  1:40894/149‑1 (MQ=255)
                                                                                                                  gATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCAc                          >  8:187741/1‑149 (MQ=255)
                                                                                                                  gATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCAc                          >  7:278232/1‑149 (MQ=255)
                                                                                                                    tgtgTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGatc                                                      <  2:48822/119‑1 (MQ=255)
                                                                                                                    tgtgTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGatc                                                      >  1:48822/1‑119 (MQ=255)
                                                                                                                        tttCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCAc                          >  5:387327/1‑143 (MQ=255)
                                                                                                                        tttCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCAc                          <  6:387327/143‑1 (MQ=255)
                                                                                                                        tttCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGa                    >  3:189628/1‑149 (MQ=255)
                                                                                                                            atcaGCAGGTCATCGAATTCCTTGAGAGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGt                >  5:6608/1‑149 (MQ=255)
                                                                                                                            atcaGCAGGTAAGCTAATTCCTTGAGTGCCTGATCGTGCTCTTCACGGCCTGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGt                <  6:287441/149‑1 (MQ=255)
                                                                                                                                gcagGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGTCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGtgatg            <  2:41052/149‑1 (MQ=255)
                                                                                                                                  agGTAATCTAATTCCTTGAGTGCCTTATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGTTTGGCAGCTAATGGCGCCCTGGCCTAAAATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGTGATGaa          <  1:294859/149‑1 (MQ=255)
                                                                                                                                       aTCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGtgatg            >  3:161802/1‑142 (MQ=255)
                                                                                                                                       aTCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGtgatg            >  4:380778/1‑142 (MQ=255)
                                                                                                                                       aTCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGtgatg            >  5:54053/1‑142 (MQ=255)
                                                                                                                                       aTCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGtgat             <  3:142615/141‑1 (MQ=255)
                                                                                                                                       aTCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGtgat             >  4:142615/1‑141 (MQ=255)
                                                                                                                                        tCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGTGATGAAGATCgg    <  8:278232/149‑1 (MQ=255)
                                                                                                                                           aaTTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGaaa                                      >  1:151920/1‑112 (MQ=255)
                                                                                                                                           aaTTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGaaa                                      <  2:151920/112‑1 (MQ=255)
                                                                                                                                           aaTTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGTGATGAAGATCGgcg  >  7:21611/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAGGTACTGGCCGCATTGCCGGCGGCGAACAGGATGTGTTTCATCAGCAGGTCATCGAATTCCTTGAGTGCCTGATCGTGGTCTTCACGCCCGGCCAATGCCATTTCCTTGAGCACGAACGGATGGCAGCGAATGGCGCCCTGGCCGAAGATCATCAGGTTGCGTGAAAGGATGTTGGCACCTTCGACCGTGATGAAGATCGGCG  >  NC_002947/2134987‑2135273

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: