Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,774,222 A→G 11.7% intergenic (+28/‑202) 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,2220AG11.7% 155.2 / 9.5 67intergenic (+28/‑202)mtlA/mtlDmannitol‑specific PTS enzyme IICBA component/mannitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (40/18);  new base G (3/5);  total (44/23)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

TCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGG  >  NC_000913/3774083‑3774364
                                                                                                                                            |                                                                                                                                               
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:184802/143‑1 (MQ=255)
                                                      tCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAAt                                                                                        >  1:787436/1‑143 (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)
                                                                                                                                          gaAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGttt    >  1:684183/1‑143 (MQ=255)
                                                                                                                                            aGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGATCGGTTACCAATGTGATGAATGTACACTTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTgg  <  2:406677/143‑1 (MQ=25)
                                                                                                                                            |                                                                                                                                               
TCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGG  >  NC_000913/3774083‑3774364

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