Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,774,215 A→G 11.1% intergenic (+21/‑209) mtlA → / → mtlD mannitol‑specific PTS enzyme IICBA component/mannitol‑1‑phosphate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,774,2150AG11.1% 166.3 / 10.4 68intergenic (+21/‑209)mtlA/mtlDmannitol‑specific PTS enzyme IICBA component/mannitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (41/19);  new base G (4/4);  total (45/23)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

CGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAG  >  NC_000913/3774074‑3774357
                                                                                                                                              |                                                                                                                                               
cGAGCACATTCAGTTTATCTCCAGCCTGACCAATTCACTGGATGATGATTTCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACAt                                                                                                                                                <  1:475888/143‑1 (MQ=255)
 aacaacacTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCCCACACCACCAGCGTGGCTGAAGTGCGGGAACTGCTGGCAGGTCGTAAGCAATCCACTCCCACACTCTCCACGTg                                                                                                                                               >  2:373345/5‑143 (MQ=255)
       aTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGgagagg                                                                                                                                         >  1:559840/1‑141 (MQ=255)
         tCATGTTATCGCCAGGCGGACCCATGCGCTGGATGATGAGTCCGTTTTCGAGCGTTTGGCACACACCACCCGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCATGTGGAGAGGGt                                                                                                                                       <  2:100032/143‑1 (MQ=255)
          cAGGGTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTg                                                                                                                                      <  1:694818/143‑1 (MQ=255)
           aGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTgg                                                                                                                                     <  1:214735/143‑1 (MQ=255)
             gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGCTGCAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAACCCCACCCTCTCCACATGGCGCAGCTgggg                                                                                                                                   >  2:349325/1‑143 (MQ=255)
             gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgggg                                                                                                                                   >  1:184412/1‑143 (MQ=255)
             gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgggg                                                                                                                                   >  1:534555/1‑143 (MQ=255)
                 tCACCTGCCTGACCAATGCACTGTATGAAGAGTCCGTGATCGAGCGTCTGGCACACCCCACCCTCGTGGATGCAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTGGGAttgg                                                                                                                               <  2:912581/143‑3 (MQ=255)
                   cccAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCGTACACCACCAGCGTGGATGAAGCGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAtt                                                                                                                             <  2:125909/142‑1 (MQ=255)
                   aCCAGCCTGACCAATGCACTGGATGATGAGGCGGTCAGCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAtt                                                                                                                             >  2:556839/1‑143 (MQ=255)
                           gACCACTGCACTGGATGATGAGTCCGTCACCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATg                                                                                                                     <  1:193892/143‑1 (MQ=255)
                                      tGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTa                                                                                                          >  1:406677/1‑143 (MQ=255)
                                       ggATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCGGCGTGGCTGAAGTGCGGGAACTGCTGGGGAGCACAGACTGACCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTAt                                                                                                         >  1:915249/1‑143 (MQ=255)
                                        gatgatGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATc                                                                                                        >  1:427383/1‑143 (MQ=255)
                                               aGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTa                                                                                                 >  1:457086/1‑143 (MQ=255)
                                               aGTCCGTCATCGAGCGGCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTa                                                                                                 <  1:302720/143‑1 (MQ=255)
                                                      cATCGAGCGTCGGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATg                                                                                          <  2:865240/143‑1 (MQ=255)
                                                       aTCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGATCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGCAGGTGGGGTTAATTGCCTGATCCGCTACGCTTATAAGGCCTACAGGATGc                                                                                         >  2:206496/1‑143 (MQ=255)
                                                         cGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCAt                                                                                       <  1:457907/143‑1 (MQ=255)
                                                          gAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCCCCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATc                                                                                      >  1:443334/1‑143 (MQ=255)
                                                          gAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGCTGCATc                                                                                      >  1:211429/1‑143 (MQ=255)
                                                           aGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGGGAAGCTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATca                                                                                     >  1:668677/1‑143 (MQ=255)
                                                           aGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATca                                                                                     >  1:582241/1‑143 (MQ=255)
                                                             cGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCCCCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATcaca                                                                                   <  2:591013/143‑1 (MQ=255)
                                                               tCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAAt                                                                                 >  1:787436/1‑143 (MQ=255)
                                                               tCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAAt                                                                                 <  1:184802/143‑1 (MQ=255)
                                                                 tGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTAGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTCCAGGATGCATCACAAttt                                                                               >  2:370202/1‑143 (MQ=255)
                                                                   gCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCCGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTGAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgt                                                                             >  1:311332/1‑143 (MQ=255)
                                                                   gCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgt                                                                             >  1:896854/1‑143 (MQ=255)
                                                                     gccgaccgccAGGGTGGATGCAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGTGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgttg                                                                           <  2:208273/139‑1 (MQ=39)
                                                                        caccacCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAt                                                                        >  1:824799/1‑143 (MQ=255)
                                                                         accaccAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAtt                                                                      >  1:806356/1‑143 (MQ=255)
                                                                           caccaGGGTGGATGAAGTGCTGGAACTGCGGGCAGGGCGTAAGCAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTg                                                                    <  1:263399/143‑1 (MQ=255)
                                                                             ccaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGCTGAGCTACGCTTCCCAGCCCTACAGGATGCCTCACAATTTGTCGAATTCGCa                                                                  >  2:889944/1‑143 (MQ=255)
                                                                             ccaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCa                                                                  >  1:78326/1‑143 (MQ=255)
                                                                                   tGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTCATCAGGCCTACAGGATGCATCACACTCTGTTGAATTTGCACGTTCt                                                            >  2:88221/1‑143 (MQ=255)
                                                                                    ggctGAAGTGCTCGAACTGCCGGCAGGTCGTAAGTAATCCAATGCGTCCCTCTCCGCATGCAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCtt                                                           <  1:763565/140‑1 (MQ=38)
                                                                                    ggATGCAGGGGTGGAACTGCTGGCGGGTGGTAATATTCGCATCTTCGGCCTCTCCACATGGAGAAGGTGGGGTTGATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCtt                                                           <  2:645716/143‑1 (MQ=25)
                                                                                       tGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTa                                                        <  2:713114/143‑1 (MQ=255)
                                                                                            tGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCg                                                   <  1:215736/143‑1 (MQ=255)
                                                                                              cTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTCGTTGAATTTGCACGTTCTTGTAGGCCGGa                                                 >  2:706807/1‑143 (MQ=255)
                                                                                                ggAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGACGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAACTCGTTGACTTTGCACGCTCTTGTAGGGCGGaga                                               >  2:563524/1‑141 (MQ=38)
                                                                                                ggAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACACTTTGTTGAATTTGCACGTTCTTGTAGGCCGGCTa                                               >  2:760048/1‑143 (MQ=255)
                                                                                                 gAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATaa                                              >  1:786206/1‑143 (MQ=255)
                                                                                                  aaCTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAg                                             <  2:457086/143‑1 (MQ=255)
                                                                                                   aCTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATGACAAGTTGTTCAATTTGCACGTCCTTGTACGCCGGATGAgg                                            >  1:727527/1‑143 (MQ=38)
                                                                                                   aCTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAgg                                            <  1:282822/143‑1 (MQ=255)
                                                                                                    ctgctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATACGCCCCGCATAACTGGCGTGCAGGATGCAGCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgc                                           <  1:251066/143‑1 (MQ=255)
                                                                                                    ctgctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCAAATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgc                                           >  1:773188/1‑143 (MQ=255)
                                                                                                        tgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCGACGCTTATCAGGCCTACAGGATGCAACAAAGTTCCGGGAAGTGGCACGTTCTTGTAGACCGGATCAGGCGCtg                                       >  2:189874/1‑142 (MQ=37)
                                                                                                              gTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgccg                                 >  1:871390/1‑143 (MQ=255)
                                                                                                                  taagtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACACATTGCTGAAGTTGCGCGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCCTc                             >  2:598478/1‑143 (MQ=37)
                                                                                                                    ggtgatccaatcgCACCCTCTGCGTATGGTGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCg                           <  2:229208/139‑1 (MQ=25)
                                                                                                                        atccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGgcgc                       <  1:556839/143‑1 (MQ=255)
                                                                                                                         tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACTTTCTTGTAGTCCGGATAAGGCCCTTACGCCGCATcctgacct                      >  2:573294/1‑137 (MQ=37)
                                                                                                                         tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCt                      >  2:7658/1‑143 (MQ=255)
                                                                                                                          ttaaccccACCTTCTCCATGTGGAGCGGGAGGAATAAATTGCCTCATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTCGCACGTTCTCGTATGCCGACTAAGGCGCGTACGCCGCATCCGGCGCTg                     >  2:282630/6‑143 (MQ=16)
                                                                                                                           caatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGc                    >  2:73338/1‑143 (MQ=255)
                                                                                                                           caatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCAAAACCGTCGCTGATTATATGACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGc                    >  1:209333/1‑143 (MQ=25)
                                                                                                                             atccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGccc                  >  1:776642/1‑143 (MQ=255)
                                                                                                                              tccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGcccc                 >  1:725052/1‑143 (MQ=255)
                                                                                                                                 ccACCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGTGTATTCACACTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCctc               >  1:910642/1‑143 (MQ=25)
                                                                                                                                 ccACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTCCTGCGCTACGCTTATCAGGCCTACAGGATGCATCACAAGTTGTTGACTTTGCACGTTCTAATCCGCTTTTCTATGCGCTTACGCCGCACCCGCCGCTGccccgc               >  2:190099/1‑141 (MQ=255)
                                                                                                                                  cACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCctct              >  2:475298/1‑143 (MQ=255)
                                                                                                                                   aCCTTCTCCATGTGGAGAGGGTGGGATTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGCTGCATCACAATTTGTTGAATTTGCACGTTCTCGTCCAAAAGCAATGGATCTTTCGCCGCATCCGGCGCTGCCCCcctc             >  1:167835/1‑143 (MQ=16)
                                                                                                                                     ccTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCGTACAGGATGCATTACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTctccg           >  2:618601/1‑142 (MQ=255)
                                                                                                                                              aTGGGGAGGGTGGGTTTAAGTGCGTGATGGGTGACGCTTATCAGGCTTACGGGATGCATCACAATTGGTTGAGTTTGCACGTGCTGGTCGCCTAGCCTGAGCTCGGAGGCCGCATCCGGCGCTGCCCCTCTCCTCACGGagag  <  2:197277/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAG  >  NC_000913/3774074‑3774357

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: