Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,732,833 T→G E447D (GAA→GAC clpB ← protein disaggregation chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,732,8330TG100.0% 216.6 / NA 57E447D (GAA→GACclpBprotein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (28/29);  total (28/29)

CCAGGTCCCCCACACGGCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGG  >  NC_000913/2732688‑2732978
                                                                                                                                                 |                                                                                                                                                 
ccAGGTCCCCCACACGGCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGgtcttc                                                                                                                                              <  1:300137/151‑1 (MQ=255)
ccAGGACCTCCACACGGCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCTAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCAATCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGgtctt                                                                                                                                               >  1:567650/1‑150 (MQ=255)
   ggTCCCCCACACGGCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGtt                                                                                                                                           >  1:809541/1‑151 (MQ=255)
   ggTCCCCCACACGGCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGtt                                                                                                                                           >  2:268804/1‑151 (MQ=255)
   ggTCCCCCACACGGCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGtt                                                                                                                                           <  2:560325/151‑1 (MQ=255)
    gTCCCCCACACGGCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTg                                                                                                                                          >  1:3734/1‑151 (MQ=255)
        cccACACGGCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCa                                                                                                                                      <  2:732808/151‑1 (MQ=255)
         ccACACGGCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCat                                                                                                                                     <  1:529180/151‑1 (MQ=255)
                gCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGa                                                                                                                              <  1:536634/151‑1 (MQ=255)
                 cGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGAc                                                                                                                             <  2:624021/151‑1 (MQ=255)
                    gCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGtt                                                                                                                          <  1:46552/151‑1 (MQ=255)
                      cTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGtttt                                                                                                                        >  1:855424/1‑151 (MQ=255)
                          tCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTa                                                                                                                     <  2:906638/150‑1 (MQ=255)
                                     tttCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTtcatc                                                                                                          <  2:374857/150‑1 (MQ=255)
                                         gCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGa                                                                                                       >  1:599427/1‑149 (MQ=255)
                                         gCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGAct                                                                                                     <  1:111142/151‑1 (MQ=255)
                                               tCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTc                                                                                                >  1:625532/1‑150 (MQ=255)
                                               tCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTc                                                                                                <  1:666733/150‑1 (MQ=255)
                                                ccAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCAt                                                                                              <  2:694265/151‑1 (MQ=255)
                                                     tCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACg                                                                                         >  2:6794/1‑151 (MQ=255)
                                                      ccGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGc                                                                                        <  1:277941/151‑1 (MQ=255)
                                                          tttAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCtgt                                                                                    >  1:672265/1‑151 (MQ=255)
                                                           ttAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCtg                                                                                     <  1:12124/149‑1 (MQ=255)
                                                            tAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCtgtt                                                                                   >  2:136311/1‑150 (MQ=255)
                                                             aaTGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCtgttgt                                                                                 >  2:449581/1‑151 (MQ=255)
                                                               tGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTc                                                                               >  1:859705/1‑151 (MQ=255)
                                                                ggTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTc                                                                               <  2:584640/150‑1 (MQ=255)
                                                                              aaaGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGgat                                                                <  1:324952/151‑1 (MQ=255)
                                                                                 gCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGgatgat                                                             >  2:677470/1‑151 (MQ=255)
                                                                                 gCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGgatgat                                                             <  2:631175/151‑1 (MQ=255)
                                                                                 gCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGgatgat                                                             <  1:569015/151‑1 (MQ=255)
                                                                                  cGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATa                                                            <  1:409459/151‑1 (MQ=255)
                                                                                          tctcTGCTTTCCACTCTTCTTCTATCTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATAcgacga                                                      <  2:129487/149‑1 (MQ=255)
                                                                                          tctcTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGAt                                                     <  1:826699/150‑1 (MQ=255)
                                                                                             ctGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCg                                                   >  2:922302/1‑149 (MQ=255)
                                                                                                 tttCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAgtcgg                                             >  2:815336/1‑151 (MQ=255)
                                                                                                   tCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAgtcggtc                                           >  1:690018/1‑151 (MQ=255)
                                                                                                   tCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAgtcgg                                             >  2:534006/1‑149 (MQ=255)
                                                                                                    ccACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAgtcggtc                                           <  1:677261/150‑1 (MQ=255)
                                                                                                        tcttcttctAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACTATCGAGTCGGTCGAGt                                       >  2:44243/1‑150 (MQ=255)
                                                                                                         cttcttctAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGt                                       <  2:200798/149‑1 (MQ=255)
                                                                                                           tcttctAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGtt                                      <  2:822321/148‑1 (MQ=255)
                                                                                                            cttctAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGCTCGAGTCGGTCGAGttcttc                                  >  1:712568/1‑151 (MQ=255)
                                                                                                            cttctAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGttcttc                                  >  1:825238/1‑151 (MQ=255)
                                                                                                                tAACTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGt                              >  2:610308/1‑151 (MQ=255)
                                                                                                                 aaCTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGtt                             >  1:210400/1‑151 (MQ=255)
                                                                                                                 aaCTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGt                              <  1:194145/150‑1 (MQ=255)
                                                                                                                   cTCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGtttt                           >  1:460753/1‑151 (MQ=255)
                                                                                                                      ggAGTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTg                          >  2:586002/1‑149 (MQ=255)
                                                                                                                         gTACTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCa                     <  2:226685/151‑1 (MQ=255)
                                                                                                                           aCTGACGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCaa                    <  1:826438/150‑1 (MQ=255)
                                                                                                                                cGTTCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGc               >  2:106882/1‑150 (MQ=255)
                                                                                                                                   tCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGGCTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGTTCGAGTCGGTCGAGTTATTCTGGTTTTGAGTCAATCTGCATAc           >  1:760086/1‑151 (MQ=255)
                                                                                                                                   tCTTTGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATAc           >  2:454353/1‑151 (MQ=255)
                                                                                                                                       tGTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGaa        <  2:627152/150‑1 (MQ=255)
                                                                                                                                        gTCGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAAt       <  2:653226/150‑1 (MQ=255)
                                                                                                                                          cGCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCt    >  1:491188/1‑151 (MQ=255)
                                                                                                                                           gCTCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTg   <  2:869000/151‑1 (MQ=255)
                                                                                                                                             tCAGGTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTgg  >  2:361297/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCAGGTCCCCCACACGGCGAGCCTGTTCAATAGCGATTTTCGCCTGTTCCAGTTCCGCTTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGG  >  NC_000913/2732688‑2732978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: