Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,706 G→C 7.8% L39V (CTC→GTC)  parA ← chromosome partition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,7060GC 7.8% 193.1 / 6.4 64L39V (CTC→GTC) parAchromosome partition protein
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base C (3/2);  total (34/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 = 5.08e-01

TTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTGC  >  NC_002947/1558‑1849
                                                                                                                                                     |                                                                                                                                               
ttGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAg                                                                                                                                                 <  3:96477/149‑1 (MQ=255)
   gCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGt                                                                                                                                               <  4:96605/148‑1 (MQ=255)
    ccGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCg                                                                                                                                             <  1:350237/149‑1 (MQ=255)
    ccGGCAGCAACTGGAAGCCTCCGTGCACGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCg                                                                                                                                             <  2:82399/149‑1 (MQ=255)
      ggCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGAt                                                                                                                                           <  8:291487/149‑1 (MQ=255)
       gcagcaACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGAt                                                                                                                                           >  4:230679/1‑148 (MQ=255)
        cagcaACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATca                                                                                                                                         <  4:146240/149‑1 (MQ=255)
           caACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATccg                                                                                                                                        >  5:302599/1‑145 (MQ=255)
           caACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATcag                                                                                                                                        >  7:66583/1‑147 (MQ=255)
           caACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATcag                                                                                                                                        >  1:264934/1‑147 (MQ=255)
            aaCTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATccgcag                                                                                                                                     <  6:302599/149‑1 (MQ=255)
                 gAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGaca                                                                                                                                 <  1:52556/148‑1 (MQ=255)
                 gAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGaca                                                                                                                                 <  2:264934/148‑1 (MQ=255)
                 gAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGaca                                                                                                                                 >  2:52556/1‑148 (MQ=255)
                  aaGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGTCCTGCTGATCGACCTCGATCCGCAGGgcaac                                                                                                                               <  1:285606/149‑3 (MQ=255)
                   aGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACacga                                                                                                                              >  6:260839/1‑149 (MQ=255)
                     ccTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACacgacg                                                                                                                            <  7:239860/149‑1 (MQ=255)
                              tCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCTAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGa                                                                                                                                            <  5:266404/124‑1 (MQ=255)
                              tCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGa                                                                                                                                            >  6:266404/1‑124 (MQ=255)
                                          aTAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAggc                                                                                                       >  7:255857/1‑149 (MQ=255)
                                                tGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAggcggc                                                                                                    >  7:8219/1‑146 (MQ=255)
                                                    ccAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGAtt                                                                                              >  2:177905/1‑148 (MQ=255)
                                                    ccAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGAtt                                                                                              >  3:319526/1‑148 (MQ=255)
                                                      aaGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGTATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGAt                                                                                           <  6:230541/149‑1 (MQ=255)
                                                           gCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAg                                                                                       <  8:8219/148‑1 (MQ=255)
                                                             aCTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAggcggc                                                                                                    <  3:323360/133‑1 (MQ=255)
                                                             aCTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAggcggc                                                                                                    >  4:323360/1‑133 (MQ=255)
                                                              cTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtgg                                                                                   <  3:230679/149‑1 (MQ=255)
                                                                ccccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGtt                                                                                 <  6:209534/149‑1 (MQ=255)
                                                                 cccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGtt                                                                                 <  7:70680/148‑1 (MQ=255)
                                                                   cGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACTGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggt                                                                                     <  3:275622/142‑1 (MQ=255)
                                                                   cGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACTGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggt                                                                                     >  4:275622/1‑142 (MQ=255)
                                                                        aGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGc                                                                          <  8:66583/148‑1 (MQ=255)
                                                                             gTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCacac                                                                    >  7:63991/1‑149 (MQ=255)
                                                                                  aGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACaccacc                                                                <  8:63991/148‑1 (MQ=255)
                                                                                    aCCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCtt                                                              >  2:153840/1‑148 (MQ=255)
                                                                                     ccGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTc                                                            <  8:255857/149‑1 (MQ=255)
                                                                                      cGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCt                                                           >  1:356438/1‑149 (MQ=255)
                                                                                      cGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCt                                                           >  7:266811/1‑149 (MQ=255)
                                                                                       gAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGAGCGAGGGCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTg                                                          <  5:51824/149‑1 (MQ=255)
                                                                                        aaTGCTCGAGCTCGTGCTTGTCCACACCGACTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCg                                                                                                                                             <  1:59610/66‑1 (MQ=255)
                                                                                        aaTGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCg                                                                                                                                             >  2:59610/1‑65 (MQ=255)
                                                                                             tCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTcgc                                                    >  6:276062/1‑149 (MQ=255)
                                                                                                 gCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCg                                                                                                                    <  8:176280/81‑1 (MQ=255)
                                                                                                 gCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCg                                                                                                                    >  7:176280/1‑81 (MQ=255)
                                                                                                     gTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGc                                              <  5:260839/147‑1 (MQ=255)
                                                                                                              cacaCCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTa                                    <  3:24281/148‑1 (MQ=255)
                                                                                                              cacaCCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAg                                   >  7:17092/1‑149 (MQ=255)
                                                                                                              cacaCCGCTGCCCATGGTGGCGTTGCCATGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTTTg                                                          >  6:38661/1‑126 (MQ=255)
                                                                                                               acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGc                                  >  7:44791/1‑149 (MQ=255)
                                                                                                                 acCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCa                                >  7:189803/1‑149 (MQ=255)
                                                                                                                   cGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAtt                              >  1:98526/1‑149 (MQ=255)
                                                                                                                   cGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAtt                              >  1:94012/1‑149 (MQ=255)
                                                                                                                         gcagcggtggCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGc                         >  4:311847/6‑149 (MQ=255)
                                                                                                                         cccATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAttgttg                          <  4:319526/148‑1 (MQ=255)
                                                                                                                          ccATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGc                         >  3:202557/1‑148 (MQ=255)
                                                                                                                               gtggCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAgtgt                   <  1:153840/149‑1 (MQ=255)
                                                                                                                                  gCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTcc                <  5:276062/149‑1 (MQ=255)
                                                                                                                                      tgtccTGCTGATCGACCTCGATCCGCAGGGCACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGt            >  7:281465/4‑149 (MQ=255)
                                                                                                                                          ctgctgATCGACCTCGATCCGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATATAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATg        >  1:225422/1‑149 (MQ=255)
                                                                                                                                           tGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGc       >  6:316378/1‑149 (MQ=255)
                                                                                                                                                aTCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTgc  >  4:51296/1‑149 (MQ=255)
                                                                                                                                                aTCGAGGTCGATCAGCAGGACACGACGCTTGGTCGAAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTgc  >  1:233940/1‑149 (MQ=255)
                                                                                                                                                 tCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCa              <  6:10654/136‑1 (MQ=255)
                                                                                                                                                 tCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCa              >  5:10654/1‑136 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTGC  >  NC_002947/1558‑1849

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: