Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,293,433 T→C 5.6% H491R (CAC→CGC)  ccmF ← heme lyase, CcmF subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,293,4330TC 5.6% 204.8 / 2.3 71H491R (CAC→CGC) ccmFheme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base C (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGGCCAGTCACCTCTTTGACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCC  >  NC_000913/2293289‑2293580
                                                                                                                                                |                                                                                                                                                   
gggCCAGTCACCTCTTTGACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAgc                                                                                                                                                 >  6:331286/1‑149 (MQ=255)
 ggCCAGTCACCTCTTTGACATCACGGAAGGTGAAGCGATATTCATGCATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAgcc                                                                                                                                                <  5:199420/149‑1 (MQ=255)
   ccAGTCACCTCTTTGACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACATTTGTCACTGCCAGCCCAAGGTGAgccgc                                                                                                                                              <  8:300601/149‑1 (MQ=255)
   ccAGTCACCTCTTTGACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAgccgc                                                                                                                                              <  5:170519/149‑1 (MQ=255)
      gTCACCTCTTTGACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGccac                                                                                                                                           >  1:281770/1‑149 (MQ=255)
       tCACCTCTTTGACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCCACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCCAGGTGAGCCGccac                                                                                                                                           >  6:358699/1‑148 (MQ=255)
         aCCTCTTTGACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGCCTTCATGCGCCCAGCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGATCCTCCACCAt                                                                                                                                        >  6:300669/1‑149 (MQ=255)
          ccTCTTTGACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATc                                                                                                                                       <  4:48199/149‑1 (MQ=255)
              tttGACATCACTGAAGGTGAAGCTATAATCATGAATATCGACGCTATCGCCTGACTTCATGCTCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAAAAATTGGCACAGACAGCCCAAGGTGAGCCGCCACCATccccc                                                                                                                                   <  7:171959/149‑1 (MQ=255)
                 gACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGcca                                                                                                                                            >  2:13130/1‑137 (MQ=255)
                 gACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGcca                                                                                                                                            <  1:13130/137‑1 (MQ=255)
                  aCATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAAt                                                                                                                                >  1:73/1‑148 (MQ=255)
                  aCATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAAt                                                                                                                                <  1:264535/148‑1 (MQ=255)
                  aCATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAAt                                                                                                                                <  5:132763/148‑1 (MQ=255)
                   cATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATaa                                                                                                                              <  6:12876/149‑1 (MQ=255)
                    aTCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCCTGCGCACATCACGCTCACCGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAAc                                                                                                                             >  8:272251/1‑149 (MQ=255)
                      cACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCAAAGGTGAGCCGCCACCATCCCCCAAAAACt                                                                                                                            >  6:71915/1‑148 (MQ=255)
                        cGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTg                                                                                                                           >  7:376444/1‑147 (MQ=255)
                         ggAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAgcc                                                                                                                                                <  8:25583/125‑1 (MQ=255)
                         ggAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAgcc                                                                                                                                                >  7:25583/1‑125 (MQ=255)
                               tGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACTATTGTCACTGCCAGCCCAAGGCGAGCCGCCACCATCCCCCAATAACTGAAGGTGGtt                                                                                                                  >  1:344381/1‑149 (MQ=255)
                               tGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGCGAGCCGCCACCATCCCCCAATAACTGAAGGTGGtt                                                                                                                  >  1:338888/1‑149 (MQ=255)
                                   gCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCg                                                                                                              <  3:187011/149‑1 (MQ=255)
                                           cATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGt                                                                                                      >  7:165353/1‑149 (MQ=255)
                                            aTGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTg                                                                                                     <  7:216439/149‑1 (MQ=255)
                                            aTGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCAAAGGTGAGCCGCCACCATCCCCAAATAACTGAAGGTGGTTTTCGTGCCGCGt                                                                                                      >  3:327887/1‑148 (MQ=255)
                                                atatCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCCCCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGaaa                                                                                                  >  3:164692/1‑148 (MQ=255)
                                                     gACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACg                                                                                              >  1:179497/1‑147 (MQ=255)
                                                     gACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGt                                                                                             <  3:151795/148‑1 (MQ=255)
                                                        gCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGc                                                                                                           >  6:370735/1‑131 (MQ=255)
                                                        gCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGc                                                                                                           <  5:370735/131‑1 (MQ=255)
                                                            tCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGc                                                                                     <  8:28727/149‑1 (MQ=255)
                                                              gCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCGGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCtt                                                                                   <  5:180337/149‑1 (MQ=255)
                                                               ccGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTc                                                                                  >  7:364464/1‑149 (MQ=255)
                                                               ccGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTc                                                                                  <  8:9184/149‑1 (MQ=255)
                                                                 ggACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCg                                                                                <  2:179497/149‑1 (MQ=255)
                                                                       cATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTg                                                                          >  7:232145/1‑149 (MQ=255)
                                                                              acaTCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCAACCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCAc                                                                   >  6:170131/1‑149 (MQ=255)
                                                                                     gCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATc                                                            >  5:185575/1‑149 (MQ=255)
                                                                                       tCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATc                                                            >  5:349848/1‑147 (MQ=255)
                                                                                          aCGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCcagca                                                       <  5:83537/149‑1 (MQ=255)
                                                                                              tataGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTcc                                                                                 <  5:21123/119‑1 (MQ=255)
                                                                                              tataGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTcc                                                                                 >  6:21123/1‑119 (MQ=255)
                                                                                              tataGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGc                                                    >  1:29165/1‑148 (MQ=255)
                                                                                              tataGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTCAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGc                                                    >  1:240178/1‑148 (MQ=255)
                                                                                                taGTTCTGGCTAAAGGCAATCCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCAt                                                 >  4:127606/1‑149 (MQ=255)
                                                                                                        gCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGc                                         >  5:50615/1‑149 (MQ=255)
                                                                                                         cTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGGGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGcc                                        >  7:314128/1‑149 (MQ=255)
                                                                                                          tAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCAtt                                                >  4:108957/1‑140 (MQ=255)
                                                                                                          tAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCAtt                                                <  3:108957/140‑1 (MQ=255)
                                                                                                           acaGGCACAGCCAACAATTGTCACTGCCAGCCCAAGGCGAGCCGCCACCATCCCCCAATAACAGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGa                                      <  2:344381/147‑1 (MQ=255)
                                                                                                            aaGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGCGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGa                                      <  2:338888/148‑1 (MQ=255)
                                                                                                                 aaTGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCg                                >  6:118258/1‑149 (MQ=255)
                                                                                                                  aTGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGt                               >  4:300695/1‑149 (MQ=255)
                                                                                                                      caacaaTTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATc                           <  7:272251/149‑1 (MQ=255)
                                                                                                                      caacaaTTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATc                           >  4:352203/1‑149 (MQ=255)
                                                                                                                       aacaaTTGTCACTGCCAGCCCCAGGTGAGCCGCCACCATCCCCCAATAACTGACTGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCg                          <  5:300669/149‑1 (MQ=255)
                                                                                                                       aacaaTTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCAcc                                 >  6:67471/1‑142 (MQ=255)
                                                                                                                       aacaaTTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCAcc                                 <  5:67471/142‑1 (MQ=255)
                                                                                                                         caaTTGTCACTGCCAGCCCAAGGTGCGCCGCCACCCTCCCCCCATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCACTTGCCAGCAACGCAATCCAGAAGGCCATTGCCAGGCCGAGCACCGTCATCGc                         >  6:149881/1‑148 (MQ=255)
                                                                                                                              gTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAAc                   >  6:184703/1‑149 (MQ=255)
                                                                                                                                     ccAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCtt            >  2:31038/1‑149 (MQ=255)
                                                                                                                                     ccAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCtt            >  1:75959/1‑149 (MQ=255)
                                                                                                                                     ccAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGGTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGc              <  5:184703/147‑1 (MQ=255)
                                                                                                                                      catcccAAGGTGAGCCGCCACCATCCCCCACTAACTGAAGGTGGTTTTCGTGCCGCGTGAAATAAGTAGCGCAGCTTCCGCAATTGCCAGCACCTCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCttt           <  2:281770/146‑1 (MQ=255)
                                                                                                                                      cAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCttt           <  8:165353/149‑1 (MQ=255)
                                                                                                                                       aGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTc          <  2:29165/149‑1 (MQ=255)
                                                                                                                                       aGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTc          <  8:232145/149‑1 (MQ=255)
                                                                                                                                          ccAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGaa       <  6:50615/149‑1 (MQ=255)
                                                                                                                                            aaGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACa     <  5:331286/149‑1 (MQ=255)
                                                                                                                                               gTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGCGCACCGTCATCGCCCCAACTTTGCTTTCGACCAGcc  >  1:331225/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGGCCAGTCACCTCTTTGACATCACGGAAGGTGAAGCGATATTCATGAATATCGACGCTATCGCCGGACTTCATGCGCACATCACGCTCAACGCTATAGTTCTGGCTAAAGGCAATGCCAACAATTGTCACTGCCAGCCCAAGGTGAGCCGCCACCATCCCCCAATAACTGAAGGTGGTTTTCGTGCCGCGTGAAATACGTAGCGCAGCTTCCGCAATTGCCAGCACCGCAATCCAGCAGGCCATTGCCAGGCCGAGCACCGTCATCGCCACAACTTTGCTTTCGAACAGCC  >  NC_000913/2293289‑2293580

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: