Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,361,3370AC41.9% 144.4 / 14.8 87G26G (GGT→GGG) rhmDL‑rhamnonate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (27/23);  new base C (0/36);  total (27/60)
Fisher's exact test for biased strand distribution p-value = 8.31e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCC  >  NC_000913/2361193‑2361484
                                                                                                                                                   |                                                                                                                                                   
gTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCAccgccg                                                                                                                                               >  2:85960/1‑151 (MQ=255)
       gCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAgc                                                                                                                                          <  2:1075058/149‑1 (MQ=255)
        cTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAg                                                                                                                                           <  1:1024493/147‑1 (MQ=255)
            gACTTCAACAACCATCGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGGCGGGATAATCTCCGCCGCCAGCGCCtttt                                                                                                                                   <  1:1201318/151‑1 (MQ=255)
              cTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTc                                                                                                                                 >  4:90810/1‑151 (MQ=255)
                  aacaacCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgc                                                                                                                              <  1:545401/150‑1 (MQ=255)
                   acaacCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgctg                                                                                                                            <  4:270314/151‑1 (MQ=255)
                        caacGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGc                                                                                                                        >  1:1208748/1‑150 (MQ=255)
                        caacGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGc                                                                                                                        >  3:123866/1‑150 (MQ=255)
                        caacGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCa                                                                                                                       <  1:436401/151‑1 (MQ=255)
                          acGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCAcc                                                                                                                     <  3:93259/151‑1 (MQ=255)
                          acGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCAcc                                                                                                                     <  3:19514/151‑1 (MQ=255)
                            gTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGc                                                                                                                   <  1:1085491/151‑1 (MQ=255)
                             tGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGcc                                                                                                                  >  1:694053/1‑151 (MQ=255)
                              gCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGcc                                                                                                                  <  1:664216/150‑1 (MQ=255)
                              gCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGcc                                                                                                                  >  2:253834/1‑150 (MQ=255)
                              gCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCa                                                                                                                 <  2:191509/151‑1 (MQ=255)
                                    aaCGTTAATGCAAAATGACTGGCGTGACTGCTCGTAACCGCGGTATTTACTCACCGGGGGGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTaa                                                                                                             <  1:628284/149‑1 (MQ=255)
                                        ttAATGCCAAATGACTGGCGTGACTGCTCTTAATCGCGGTATTTACTCACCGGGGTGG‑A‑AATATGATCGTCTAC‑CCAGTGATTCCCCCCCTGGTCGTGATAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCa                                                                                                         <  2:828616/149‑1 (MQ=255)
                                        taaaTGCAAAATGACGGGCGTGCCTGCTCGTAATCGCGGTATTTCCTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATCCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAg                                                                                                        <  4:221023/148‑1 (MQ=255)
                                         tAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgc                                                                                                      >  4:349172/1‑151 (MQ=255)
                                          aaTGCAAAATGACTGGCGTGCCTGCTCGAATTCGCGGTTTTTCCCCATCGGGGTGG‑A‑AATATGATCTTCTAT‑CCATTTTTTCGCCCCCTGGCCGTGATAACCCCCGCCGCCAGCGCCTTTTTCTGCTGCCGCACCGCCAGTAAACCAGgc                                                                                                      <  2:92019/150‑1 (MQ=255)
                                           aTGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgcgc                                                                                                    >  2:1007262/1‑151 (MQ=255)
                                            tGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgcgc                                                                                                    >  2:443425/1‑150 (MQ=255)
                                              caaaaTGCCGGGCGTGACTGCTCGTAATCGCGTTTTTTACTCATCGGGGTGG‑A‑AATATGATCGTCTAC‑CCAGTGATCCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGaa                                                                                                 <  1:150527/149‑1 (MQ=255)
                                                aaaTGACTGGCGTGATTGCTATAAATCGCGGTATTTACTCA‑CGGGGTTG‑CAAATAGGA‑CCTTTATCCCGGTGATCCCCCCCCTGTTGGTGAAAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGaa                                                                                                 <  2:1148237/149‑1 (MQ=255)
                                                aaaTGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTTATTCGCCCCCTGGTCGTGATAACCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGaa                                                                                                 <  3:128064/149‑1 (MQ=255)
                                                  aTGACTGGCGTGCCGGCTCTAAATCGCGGTATTTCCTCATCGGGGTGG‑A‑AATATGATCGTCTAC‑CCAGTGTTTCCCCCCCTGATCGTAAAAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTg                                                                                             <  2:220240/151‑1 (MQ=255)
                                                         gCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAtt                                                                                      >  1:350580/1‑151 (MQ=255)
                                                           gTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAtttt                                                                                    >  4:338760/1‑151 (MQ=255)
                                                            tGACTGCTCGTAATCGCGGTTTTTACTCATCGGGGGGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAtttt                                                                                    <  1:797690/150‑1 (MQ=255)
                                                            tGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAttttt                                                                                   >  1:827692/1‑151 (MQ=255)
                                                                 gCTCGTAATCGCGTTTTTTAC‑CCTCCGGGTGGCC‑AATATGTTCGGGTAG‑CCAGGTTTTCCCCCCCTGGTCGTGATAATCCCCGCCGCGAGCCCCTTTTTCTGCTGTCCCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTg                                                                                  <  3:267532/147‑1 (MQ=255)
                                                                   tCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAggg                                                                            >  2:280947/1‑151 (MQ=255)
                                                                    cGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGt                                                                           <  4:97007/151‑1 (MQ=255)
                                                                    cGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGt                                                                           >  4:31448/1‑151 (MQ=255)
                                                                    cGTAACCGCGGTATTTACTCATGGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAggg                                                                            <  1:1258096/150‑1 (MQ=255)
                                                                     gTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTc                                                                          >  3:145187/1‑151 (MQ=255)
                                                                          cgcgGTATTTCCTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatga                                                                      <  3:39170/150‑1 (MQ=255)
                                                                            cgGTTTTTACTCATGGGGGTGG‑A‑AATATGATCTTCTAC‑CCATTAATCCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatga                                                                      <  1:1247812/148‑1 (MQ=255)
                                                                             ggTATTTACTCAGCGGGGTGG‑A‑AATATGATCGTCTAC‑CCAGTGATTCGCCCCCGGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGt                                                                   <  2:859035/150‑1 (MQ=255)
                                                                             ggTATTTACCCATCGGGGTGG‑A‑AATATGATCGTCTAC‑CCAGTGATTCGCCCCCTGGTCGGGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGt                                                                   <  1:1038191/150‑1 (MQ=255)
                                                                               tATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTtct                                                                >  1:699030/1‑151 (MQ=255)
                                                                               tATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTtct                                                                >  1:502026/1‑151 (MQ=255)
                                                                                     cccATCGGGGTGG‑C‑AATATGACCGCCAAT‑CCAGTGATTCCCCCCGTGGGCGGGATAATCCCCGCCGCCGGCGCTTTTTTTTGCTGGCGCGCCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTtct                                                                <  2:785172/143‑1 (MQ=255)
                                                                                      ccaccGGGGTGG‑C‑AATATGTTCGTCTAT‑CCGGTGTTTCCCCCCCTGTCCGTAATAACCCCCGCGGCCAGCGCTTTTTTCGGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCAtt                                                           <  2:291447/145‑1 (MQ=255)
                                                                                       cATGGGGGTGG‑A‑AATATTATCGTCTAC‑CCGTTGATTCGCCCCCTGGTCGTGATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCAt                                                            <  2:136411/147‑1 (MQ=255)
                                                                                       cATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGTTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGt                                                         <  1:522074/150‑1 (MQ=255)
                                                                                       cATCGGGGTGC‑A‑AATATGTTCGTCTAC‑CCAGTGATTCGCCCCCTGGCCGTGATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGtt                                                        <  2:173410/151‑1 (MQ=255)
                                                                                        acgggggTGG‑A‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCAtt                                                           <  1:655284/144‑1 (MQ=255)
                                                                                        aTCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTAAAAATCCCCGCCGCCAGCTCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGtt                                                        <  2:216084/150‑1 (MQ=255)
                                                                                          cGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATg                                                     >  2:719910/1‑151 (MQ=255)
                                                                                          cGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATg                                                     >  1:292927/1‑151 (MQ=255)
                                                                                           ggggggG‑C‑AATATGAGCTTCTAC‑CCAGTGATTCCCCCCCTGGTCGTAAAAATCCCCGCCGCGAGCGTCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTAt                                                      <  1:640387/149‑1 (MQ=255)
                                                                                           ggggTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATg                                                     >  4:373046/1‑150 (MQ=255)
                                                                                           ggggTGC‑A‑AATATAAGCGTCTAC‑GCAGTGATTCGCCCCCTGTTCGGAAAAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTcc                                                              <  2:665307/141‑1 (MQ=255)
                                                                                             ggTGG‑A‑AATTTTACCTTCTAC‑CAAGTGATCCCCCCCCGGGTCGTGATAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAgg                                                  <  2:755422/151‑1 (MQ=255)
                                                                                             ggTGG‑A‑AATATGATCTTCTAC‑CCAGTGATCCGCCCCCTGGTCGTGATAATCCCCGCGGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAgg                                                  <  1:903226/151‑1 (MQ=255)
                                                                                                gg‑C‑AATATGATCGTCTAT‑CCAGGGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCtt                                               <  1:1100713/151‑1 (MQ=255)
                                                                                                 g‑C‑AATATGTTCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGAAAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCtt                                               <  4:384481/150‑1 (MQ=255)
                                                                                                   c‑AATATGACCGTCTAC‑CCAGTAATCCGCCCCCTGTCCGTAATAATCACCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTg                                              <  1:817637/150‑1 (MQ=255)
                                                                                                     aaTATGTTCGTTTAT‑CCAGTGATTCGCCCCCTTGTCGTGATAACCCCCGCCCCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTg                                              <  2:133392/149‑1 (MQ=255)
                                                                                                     aaTATGATCGTCTAT‑CCGGTGATTCGCCCCCTGGTCGTAATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTa                                            <  4:130145/151‑1 (MQ=255)
                                                                                                       taggATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTTTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTaa                                           <  1:880076/147‑1 (MQ=255)
                                                                                                        atGCTCGTCTAT‑CCAGTGTTCCCCCCCCGGGTCGTGATAAACCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGt                                         <  2:871203/151‑1 (MQ=255)
                                                                                                          gttCGTCTTT‑CCTGTGTTTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAg                                          <  3:322115/146‑1 (MQ=255)
                                                                                                          gATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCa                                       <  1:662502/151‑1 (MQ=255)
                                                                                                             cGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCaaa                                     <  4:174124/150‑1 (MQ=255)
                                                                                                                     ccggtgtTTCGCCCCCTGGTCGGGAAAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGa                                 <  1:21723/144‑1 (MQ=255)
                                                                                                                     ccggTGATTCGCCCCCTGGTCGTAATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACt                               <  1:371834/146‑1 (MQ=255)
                                                                                                                     ccAGTGATCCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACttt                             <  4:417978/151‑1 (MQ=255)
                                                                                                                       attgttccccccccTGGTGGGGATAATAACCGCCGCCAGCCTCTTTTTTTGCTGCCGCCCCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAgg                                   <  2:80959/134‑1 (MQ=255)
                                                                                                                        gTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTcc                          >  1:1013534/1‑151 (MQ=255)
                                                                                                                         tgtttCGCCCCTGGGTCGTAATACTCACCGCCGCAAGCGCCTTTTTCTGCTTCCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTcc                          <  4:438208/147‑1 (MQ=255)
                                                                                                                           atgcccccccTGGGCGGGATATCCCCCGGCGCCAGCGCCTTTTTCTGCTTTCTCGCCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTc                           <  1:1026421/142‑1 (MQ=255)
                                                                                                                           aTTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATc                       >  1:109959/1‑151 (MQ=255)
                                                                                                                             tCGCCCCCGGTTCGTGAAAAGCCCCGCCGCGAGCGCCTTTTTCTGGGGGCACGCCGCCAGTAAACCAGGCGCGAACCTTTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATc                       <  2:1089000/149‑1 (MQ=255)
                                                                                                                                cccccTGGTCGTGATAATCACCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATcc                      <  3:196458/147‑1 (MQ=255)
                                                                                                                                cccccTGGTCGTGATAACCACCGCCCCCGCCACCTTTTTCGGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCaa                   <  1:956607/150‑1 (MQ=255)
                                                                                                                                cccccGGGTCGTGATAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGATCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCa                    <  3:248386/149‑1 (MQ=255)
                                                                                                                                 ccccTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAAc                  >  4:391838/1‑150 (MQ=255)
                                                                                                                                     tGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACg             >  1:1103288/1‑151 (MQ=255)
                                                                                                                                       gttgtaaTAACCACCACCCCCGCCCCCTTTTTCTGCTATCGCGCCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCa                    <  2:211730/136‑1 (MQ=255)
                                                                                                                                       gTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTg           >  1:1050275/1‑151 (MQ=255)
                                                                                                                                          gtaaTAATCCCCGCCGCCAGCGCCTTTTTCTGCTGCCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGaa         <  2:1071778/147‑1 (MQ=255)
                                                                                                                                                aTCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATccc  >  4:42708/1‑151 (MQ=255)
                                                                                                                                                aTCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATccc  >  1:650656/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGG‑C‑AATATGATCGTCTAT‑CCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCC  >  NC_000913/2361193‑2361484

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑