Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,294,209 G→A 6.4% Y232Y (TAC→TAT)  ccmF ← heme lyase, CcmF subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,294,2090GA 6.4% 187.3 / 2.5 62Y232Y (TAC→TAT) ccmFheme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base A (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01

TGATTGCCAGCAGTAATGTCCACGCTTTGAAGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAA  >  NC_000913/2294061‑2294356
                                                                                                                                                    |                                                                                                                                                   
tgcttGCCAGCAGTAAGGACGTCGGTTTGAAGCTGGGGGGGTTTGCAGGGACCGCCCGTGAGTGCATCAGGGCAGGCCCCACCCGCGCCGGCATAAACGTGGCGTTTTCTACCGGATGCCAGAACCACCAGCCACCCCAGCTGAGTTCg                                                                                                                                                     <  2:105753/146‑1 (MQ=255)
               aTGTCCACGCTTTGAAGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGc                                                                                                                                       <  3:186067/148‑1 (MQ=255)
                  tCCACGCTTTGAAGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTCCCGGATCCCAGCACCCCCAGCCCCCCCAGCCGAGTTCGTAATAGGCCCATGCGGaa                                                                                                                                   >  1:54217/1‑149 (MQ=255)
                  tCCACGCGTTGAAGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATGCCAGAACCACCAGCCCCCCCAGCCGCGTTCGTAATAGGCCCCTGCggca                                                                                                                                   >  2:177081/1‑147 (MQ=255)
                     aCGCTTTGAAGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCCGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCg                                                                                                                                >  6:222377/1‑149 (MQ=255)
                     aCGCTTTGAAGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCg                                                                                                                                <  8:61402/149‑1 (MQ=255)
                     aCGCTTTGAAGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGACCCCACCAGCCACGGCTTAAACGAGGCGTTTTCTACCGGATCCCAGACCCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAAcc                                                                                                                                 >  1:252844/1‑148 (MQ=255)
                         tttGAAGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCa                                                                                                                            >  4:375181/1‑149 (MQ=255)
                           tGAAGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCAc                                                                                                                           <  6:295795/148‑1 (MQ=255)
                             aaGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGAt                                                                                                                        >  4:154013/1‑149 (MQ=255)
                                cTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGcc                                                                                                                     >  5:143065/1‑149 (MQ=255)
                                     gcgTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCg                                                                                                                >  5:339975/1‑149 (MQ=255)
                                       gttgttCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTc                                                                                                              <  3:285181/149‑1 (MQ=255)
                                        ttgttCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCa                                                                                                             <  3:320763/149‑1 (MQ=255)
                                           ttCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGa                                                                                                          <  1:372921/149‑1 (MQ=255)
                                                  aCCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCa                                                                                                   >  2:239253/1‑149 (MQ=255)
                                                  aCCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCa                                                                                                   >  4:253745/1‑149 (MQ=255)
                                                   ccGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCAc                                                                                                  <  3:327849/149‑1 (MQ=255)
                                                    cGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCa                                                                                                                                          <  7:183704/108‑1 (MQ=255)
                                                          tgagtgCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCATAATAGGCCCATGCGGaa                                                                                                                                   <  4:258123/109‑1 (MQ=255)
                                                          tgagtgCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCATAATAGGCCCATGCGGaa                                                                                                                                   <  4:253028/109‑1 (MQ=255)
                                                          tgagtgCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCATAATAGGCCCATGCGGaa                                                                                                                                   >  3:258123/1‑109 (MQ=255)
                                                          tgagtgCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCATAATAGGCCCATGCGGaa                                                                                                                                   >  3:253028/1‑109 (MQ=255)
                                                              tgCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGt                                                                                       <  3:253745/149‑1 (MQ=255)
                                                                cATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTca                                                                                     >  7:356208/1‑148 (MQ=255)
                                                                    aGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACg                                                                                  <  3:154013/148‑1 (MQ=255)
                                                                        cAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGCACCACCAGCCCCCCCAGCCGAGGTCGTAATAGGCCCATGCGGAACCGCGCACGATGCCGCGCGTCAGGAAGATCCACGCTGCCAGCGTCCACggcc                                                                               >  6:69547/1‑145 (MQ=255)
                                                                           tCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGt                                                                           >  4:365270/1‑148 (MQ=255)
                                                                             cccACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTaa                                                                         >  4:289101/1‑148 (MQ=255)
                                                                               caccaGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAAc                                                                      >  5:234377/1‑149 (MQ=255)
                                                                               caccaGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAAc                                                                      <  5:155426/149‑1 (MQ=255)
                                                                               caccaGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAAc                                                                      <  7:129106/149‑1 (MQ=255)
                                                                                        cGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAg                                                             >  7:255315/1‑149 (MQ=255)
                                                                                            aTAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGCGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCt                                                         <  2:252844/149‑1 (MQ=255)
                                                                                            aTAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCt                                                         >  1:30069/1‑149 (MQ=255)
                                                                                            aTAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCt                                                         <  8:255315/149‑1 (MQ=255)
                                                                                             tAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTg                                                        >  7:368262/1‑149 (MQ=255)
                                                                                               aaCGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCACATGCGGAACCGAGCACGATGCCGAGCTTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGt                                                       <  2:84056/148‑1 (MQ=255)
                                                                                                      cGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGGCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGc                                               <  5:69547/149‑1 (MQ=255)
                                                                                                        ttttCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACTCGCATAAGTGCTGTCCAGACGccc                                             <  6:143065/149‑1 (MQ=255)
                                                                                                              aCCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAg                                       >  1:36425/1‑149 (MQ=255)
                                                                                                               ccGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGTCGAGTAAAACGCGCATAAGTGCTGTCCAGCCGCCCGCTCAg                                       >  2:31370/1‑148 (MQ=255)
                                                                                                               ccGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCa                                        <  3:374215/147‑1 (MQ=255)
                                                                                                                 ggATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCaa                                    <  8:368262/149‑1 (MQ=255)
                                                                                                                  gATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCa                                     >  4:276274/1‑147 (MQ=255)
                                                                                                                      ccAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAaagaa                                <  3:276274/148‑1 (MQ=255)
                                                                                                                      ccAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAaagaa                                <  1:360673/148‑1 (MQ=255)
                                                                                                                      ccAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAaagaa                                <  3:289101/148‑1 (MQ=255)
                                                                                                                       cAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGAAGCCGAGCGTCAGGAAGATCCACGCTTCCAGCGTCCACCGACGAGTAAAACGCCCAAAAGTGCTGTCCAGACGCCCGCTCAGCAaagaag                               >  4:299215/1‑148 (MQ=255)
                                                                                                                        aGAACCACCAGCCACCCCGGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGc                              >  4:256462/1‑148 (MQ=255)
                                                                                                                            ccaccaGCCCCCCCAGCCGAGTTCGTAATAGGCCCATGCTGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGACTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATgg                         >  5:134055/1‑149 (MQ=255)
                                                                                                                             caccaGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATgg                         >  5:386965/1‑148 (MQ=255)
                                                                                                                             caccaGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATgg                         >  5:376976/1‑148 (MQ=255)
                                                                                                                             caccaGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATgg                         >  2:126469/1‑148 (MQ=255)
                                                                                                                             caccaGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCAACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATgg                         >  7:61456/1‑148 (MQ=255)
                                                                                                                                caGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGACCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCaa                      >  1:116564/1‑148 (MQ=255)
                                                                                                                                 aGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCaaaa                    <  2:30069/149‑1 (MQ=255)
                                                                                                                                 aGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCaaaa                    <  3:375181/149‑1 (MQ=255)
                                                                                                                                 aGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCaaaa                    <  1:239253/149‑1 (MQ=255)
                                                                                                                                       cccAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAAc              <  6:376976/149‑1 (MQ=255)
                                                                                                                                          aGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGcc           <  2:36425/149‑1 (MQ=255)
                                                                                                                                            ccGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCAc         <  6:234377/149‑1 (MQ=255)
                                                                                                                                                   cGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCCCGGACGAGTAAACCGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGaaa  >  1:70700/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGATTGCCAGCAGTAATGTCCACGCTTTGAAGCTGGCGCGTTGTTCAGTGACCGCCAGTGAGTGCATCAGCGCAGTCCCCACCAGCCACGGCATAAACGAGGCGTTTTCTACCGGATCCCAGAACCACCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAA  >  NC_000913/2294061‑2294356

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