Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,824,478 (C)6→7 intergenic (‑132/+13) pncC ← / ← mltB NMN aminohydrolase/membrane‑bound lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,824,4721.C98.4% 214.5 / ‑3.3 61intergenic (‑126/+19)pncC/mltBNMN aminohydrolase/membrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (25/35);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

TGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAG  >  NC_000913/2824337‑2824602
                                                                                                                                        |                                                                                                                                   
tGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccc                                                                                                                                <  1:498983/140‑1 (MQ=255)
tGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccc                                                                                                                                <  2:299057/140‑1 (MQ=255)
tGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccc                                                                                                                                <  1:929250/140‑1 (MQ=255)
tGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccc                                                                                                                                >  2:40876/1‑140 (MQ=255)
  tcagtcaTCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccccc                                                                                                                              <  2:20737/140‑1 (MQ=255)
  tcagtcaTCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccccc                                                                                                                              <  2:73834/140‑1 (MQ=255)
  tcagtcaTCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccccc                                                                                                                              <  1:366829/140‑1 (MQ=255)
    agtcaTCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGccccccct                                                                                                                            >  1:239518/1‑138 (MQ=255)
     gtcaTCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGccccccctc                                                                                                                           <  2:1072654/140‑4 (MQ=255)
     gtcaTCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGccccccctc                                                                                                                           >  2:508938/1‑137 (MQ=255)
        atcaCAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCATAGGGAAGCTCACGCCCCCCCTCTTg                                                                                                                         <  1:660580/140‑1 (MQ=255)
          cacaACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTa                                                                                                                       >  1:441706/1‑140 (MQ=255)
              aCCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCTCCCTCTTGTAAATa                                                                                                                   <  2:142002/140‑1 (MQ=255)
               ccTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAg                                                                                                                  <  2:784736/140‑1 (MQ=255)
                ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGc                                                                                                                 <  1:912483/140‑1 (MQ=255)
                ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGc                                                                                                                 >  2:81143/1‑140 (MQ=255)
                   tGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTAc                                                                                                              >  2:510371/1‑140 (MQ=255)
                       agtcaAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTa                                                                                                          >  2:926599/1‑140 (MQ=255)
                         tcaAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACt                                                                                                        >  1:153604/1‑140 (MQ=255)
                         tcaAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACt                                                                                                        <  1:913462/140‑1 (MQ=255)
                                   aTGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCgcgc                                                                                                   <  2:156839/135‑1 (MQ=255)
                                   aTGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCgcgc                                                                                                   >  1:156839/1‑135 (MQ=255)
                                     gCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCgcgccagcgcc                                                                                            <  2:604595/140‑1 (MQ=255)
                                         gACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACgg                                                                                        >  2:269430/1‑140 (MQ=255)
                                          aCAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGc                                                                                       <  1:855965/140‑1 (MQ=255)
                                           cAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCt                                                                                      <  1:269430/140‑1 (MQ=255)
                                             aGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTg                                                                                    <  1:40876/140‑1 (MQ=255)
                                             aGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTg                                                                                    <  1:926599/140‑1 (MQ=255)
                                                 aGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCt                                                                                <  2:640830/140‑1 (MQ=255)
                                                  gCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTa                                                                               <  2:153604/140‑1 (MQ=255)
                                                        tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACgg                                                                                        >  1:183882/1‑125 (MQ=255)
                                                        tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACgg                                                                                        <  2:183882/125‑1 (MQ=255)
                                                        tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                                                         >  1:877777/1‑140 (MQ=255)
                                                        tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                                                         >  2:678945/1‑140 (MQ=255)
                                                        tCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                                                         >  2:292505/1‑140 (MQ=255)
                                                              ggAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCg                                                                   <  2:1070897/140‑1 (MQ=255)
                                                              ggAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCg                                                                   <  1:576044/140‑1 (MQ=255)
                                                                  gATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCAt                                                               <  1:81143/140‑1 (MQ=255)
                                                                     ggTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGc                                                            >  1:172385/1‑140 (MQ=255)
                                                                     ggTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGc                                                            <  2:538423/140‑1 (MQ=255)
                                                                       tgtTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGt                                                          <  1:292505/140‑1 (MQ=255)
                                                                       tgtTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGt                                                          <  1:678945/140‑1 (MQ=255)
                                                                         ttGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTaa                                                        >  2:395599/1‑140 (MQ=255)
                                                                         ttGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTaa                                                        >  2:968332/1‑140 (MQ=255)
                                                                                tGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGc                                                 >  1:172134/1‑140 (MQ=255)
                                                                                tGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGc                                                 <  2:172134/140‑1 (MQ=255)
                                                                                 gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTaa                                                        <  1:1049390/132‑1 (MQ=255)
                                                                                 gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTaa                                                        >  2:1049390/1‑132 (MQ=255)
                                                                                 gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCt                                                >  2:336629/1‑140 (MQ=255)
                                                                                 gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCt                                                <  2:441706/140‑1 (MQ=255)
                                                                                          cGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTa                                       <  1:472684/140‑1 (MQ=255)
                                                                                            aggaggATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAAc                                     >  2:107808/1‑140 (MQ=255)
                                                                                            aggaggATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAAc                                     >  2:773040/1‑140 (MQ=255)
                                                                                               aggaTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACggg                                  <  1:107808/140‑1 (MQ=255)
                                                                                               aggaTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACggg                                  <  1:773040/140‑1 (MQ=255)
                                                                                                 gaTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTg                                >  1:344097/1‑140 (MQ=255)
                                                                                                  aTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                               <  2:541813/140‑1 (MQ=255)
                                                                                                  aTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                               <  1:381057/140‑1 (MQ=255)
                                                                                                  aTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                               <  2:115124/140‑1 (MQ=255)
                                                                                                     cAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                               >  2:61925/1‑137 (MQ=255)
                                                                                                     cAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                               <  1:61925/137‑1 (MQ=255)
                                                                                                        aaaaaaGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGgtgt                         >  1:176472/1‑140 (MQ=255)
                                                                                                             aGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAg                    >  2:688434/1‑140 (MQ=255)
                                                                                                                       ccAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc          <  1:688434/140‑1 (MQ=255)
                                                                                                                           agGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTa      <  1:968332/140‑1 (MQ=255)
                                                                                                                           agGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTa      <  2:239518/140‑1 (MQ=255)
                                                                                                                           agGGAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTa      <  2:933109/140‑1 (MQ=255)
                                                                                                                            gggAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTAc     >  1:50583/1‑140 (MQ=255)
                                                                                                                            gggAAGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTAc     <  2:877777/140‑1 (MQ=255)
                                                                                                                               aaGCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAg  <  1:395599/140‑1 (MQ=255)
                                                                                                                                 gCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTa      >  1:412756/1‑134 (MQ=255)
                                                                                                                                 gCTCACGCCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTa      <  2:412756/134‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
TGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTACCAG  >  NC_000913/2824337‑2824602

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