Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,593,004 C→T 5.3% intergenic (+42/+68) ypfN → / ← ypfH putative membrane protein, UPF0370 family/palmitoyl‑CoA esterase activity, uncertain physiological substrate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,593,0040CT 5.3% 303.0 / 3.7 91intergenic (+42/+68)ypfN/ypfHputative membrane protein, UPF0370 family/palmitoyl‑CoA esterase activity, uncertain physiological substrate
Reads supporting (aligned to +/- strand):  ref base C (43/43);  new base T (3/2);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

AAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATA  >  NC_000913/2592858‑2593152
                                                                                                                                                  |                                                                                                                                                    
aaGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCAt                                                                                                                                                    >  1:880177/1‑149 (MQ=255)
 aGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATc                                                                                                                                                   >  1:876555/1‑149 (MQ=255)
  gAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                                                                  >  2:525478/1‑149 (MQ=255)
  gAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATcc                                                                                                                                                  <  1:760294/149‑1 (MQ=255)
     aTTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCACCCGAAGAAGTAAGCTGAGTCCTAAACCGTAGGCCGAATAAGGCGTTAACGCCCCATCCGGc                                                                                                                                               >  2:798492/1‑149 (MQ=255)
     aTTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                                                               <  1:525478/149‑1 (MQ=255)
         ctcGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATGCGGCGTTAACGCCTTATCCGGcctac                                                                                                                                           <  3:648242/149‑5 (MQ=255)
         ctcGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATg                                                                                                                                           <  4:189385/149‑1 (MQ=255)
         ctcGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATg                                                                                                                                           <  4:153191/149‑1 (MQ=255)
         ctcGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATg                                                                                                                                           <  3:178091/149‑1 (MQ=255)
          tcGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATgg                                                                                                                                          >  4:805927/1‑149 (MQ=255)
          tcGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATgg                                                                                                                                          >  2:146042/1‑149 (MQ=255)
              caacaaGCCGGAGCTTCCCCCGCATCGTGATTTCAAGGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGaa                                                                                                                                      >  4:835745/1‑149 (MQ=255)
              caacaaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGaa                                                                                                                                      <  2:450227/149‑1 (MQ=255)
              caacaaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGaa                                                                                                                                      >  1:207275/1‑149 (MQ=255)
              caacaaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGaa                                                                                                                                      >  1:43071/1‑149 (MQ=255)
                acaaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACg                                                                                                                                    >  3:875463/1‑149 (MQ=255)
                   aGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATg                                                                                                                                 <  1:529768/149‑1 (MQ=255)
                      cGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCt                                                                                                                              <  1:882459/149‑1 (MQ=255)
                            ttCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCg                                                                                                                        >  3:143348/1‑149 (MQ=255)
                                cccGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGc                                                                                                                    >  3:550149/1‑149 (MQ=255)
                                        gTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGt                                                                                                            >  2:835232/1‑149 (MQ=255)
                                               gaaCGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCa                                                                                                     >  2:758888/2‑149 (MQ=255)
                                                aaCGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAg                                                                                                    <  1:801003/149‑1 (MQ=255)
                                                  cGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTTAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGc                                                                                                  <  3:835745/149‑1 (MQ=255)
                                                     tacgTGGGACGATGACGACGACTAGCCGAGAAAGGATGGACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACTATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAg                                                                                               >  2:822160/4‑149 (MQ=37)
                                                     tAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAg                                                                                               <  4:114054/149‑1 (MQ=255)
                                                     tAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAg                                                                                               <  3:519084/149‑1 (MQ=255)
                                                         tGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTaa                                                                                           >  1:231031/1‑149 (MQ=255)
                                                         tGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATAAGACGCGTGAGCGTCGCATCAGGcatcgttc                                                                                           <  1:565475/149‑8 (MQ=37)
                                                          gggACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAAt                                                                                          <  4:268802/149‑1 (MQ=255)
                                                           ggACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATc                                                                                         <  2:164934/149‑1 (MQ=255)
                                                            gACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACTCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCt                                                                                        <  2:208956/149‑1 (MQ=255)
                                                            gACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCt                                                                                        <  2:207275/149‑1 (MQ=255)
                                                              cGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGt                                                                                      <  1:715395/149‑1 (MQ=255)
                                                               gATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTc                                                                                     >  4:758998/1‑149 (MQ=255)
                                                               gATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTc                                                                                     >  4:835758/1‑149 (MQ=255)
                                                               gATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTc                                                                                     <  1:165237/149‑1 (MQ=255)
                                                                  gacgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCt                                                                                    <  1:133232/147‑1 (MQ=255)
                                                                  gacgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCt                                                                                    >  2:133232/1‑147 (MQ=255)
                                                                   acgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCtttt                                                                                 <  4:289834/149‑1 (MQ=255)
                                                                    cgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTa                                                                                <  2:483380/149‑1 (MQ=255)
                                                                      acgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAca                                                                              <  4:633911/149‑1 (MQ=255)
                                                                       cgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTc                                                                                                           >  1:102810/1‑119 (MQ=255)
                                                                       cgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTc                                                                                                           <  2:102810/119‑1 (MQ=255)
                                                                               cGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCaa                                                                     >  3:665600/1‑149 (MQ=255)
                                                                                gAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCtttt                                                                                 >  2:454495/1‑136 (MQ=255)
                                                                                gAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCtttt                                                                                 <  1:454495/136‑1 (MQ=255)
                                                                                gAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAAt                                                                    <  1:28817/149‑1 (MQ=255)
                                                                                gAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAAt                                                                    <  3:768546/149‑1 (MQ=255)
                                                                                gAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAAt                                                                    >  1:757415/1‑149 (MQ=255)
                                                                                gAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAACGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAAt                                                                    <  1:459531/149‑1 (MQ=255)
                                                                                   aaaGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATcac                                                                 <  4:457877/149‑1 (MQ=255)
                                                                                      ggATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCAcatc                                                              <  1:835232/149‑1 (MQ=255)
                                                                                        aTCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCAcatcat                                                            >  2:306845/1‑149 (MQ=255)
                                                                                           aaCCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCAcatcatcat                                                         >  3:175931/1‑149 (MQ=255)
                                                                                             ccGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCg                                                       >  3:683767/1‑149 (MQ=255)
                                                                                             ccGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCg                                                       >  1:233548/1‑149 (MQ=255)
                                                                                               gaagaaGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGcc                                                     >  3:336702/1‑149 (MQ=255)
                                                                                               gaagaaGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGcc                                                     >  4:567290/1‑149 (MQ=255)
                                                                                               gaagaaGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGcc                                                     >  4:342122/1‑149 (MQ=255)
                                                                                                aagaagTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCaa                                                                     >  4:777239/1‑132 (MQ=255)
                                                                                                aagaagTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCaa                                                                     <  3:777239/132‑1 (MQ=255)
                                                                                                aagaagTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCa                                                    >  1:737025/1‑149 (MQ=255)
                                                                                                 agaagTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAg                                                   <  4:808479/149‑1 (MQ=255)
                                                                                                   aagTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGGTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGc                                                 <  4:70890/149‑1 (MQ=255)
                                                                                                      tAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCAcgc                                                                                                              <  2:330373/85‑1 (MQ=255)
                                                                                                      tAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCAcgc                                                                                                              >  1:330373/1‑85 (MQ=255)
                                                                                                                  tCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGcc                                         >  4:257741/1‑142 (MQ=255)
                                                                                                                  tCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGcc                                         <  3:257741/142‑1 (MQ=255)
                                                                                                                  tCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAgc                                  >  3:447461/1‑149 (MQ=255)
                                                                                                                  tCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAgc                                  >  4:126398/1‑149 (MQ=255)
                                                                                                                    aaCCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAgcgc                                <  3:37581/149‑1 (MQ=255)
                                                                                                                       cGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTtc                             >  2:158114/1‑149 (MQ=255)
                                                                                                                        cattgCCGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTtcg                            >  2:332697/5‑149 (MQ=255)
                                                                                                                         attgCCGGATGCGGCGTTAACGGCCTATCCGGCAATGCTGAACGATGCCTAATGCGACGCACCCGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACCTGATTTCAATCAAATCATCATCGCCAGGCTTACTCCCACTTAGCGCTtcga                           >  2:446417/4‑148 (MQ=255)
                                                                                                                          aGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTtcgtc                          <  2:276720/149‑1 (MQ=255)
                                                                                                                            gCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGa                        >  4:112398/1‑149 (MQ=255)
                                                                                                                            gCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGa                        >  3:307756/1‑149 (MQ=255)
                                                                                                                              cGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGaaa                      <  3:758998/149‑1 (MQ=255)
                                                                                                                              cGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGaaa                      <  2:737025/149‑1 (MQ=255)
                                                                                                                                gATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAaata                    >  1:62661/1‑149 (MQ=255)
                                                                                                                                 aTAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAaataa                   <  1:146042/149‑1 (MQ=255)
                                                                                                                                        gTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCg            >  3:33579/1‑149 (MQ=255)
                                                                                                                                         ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCgg           >  1:850219/1‑149 (MQ=255)
                                                                                                                                         ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCgg           >  2:530765/1‑149 (MQ=255)
                                                                                                                                         ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCgg           >  3:644987/1‑149 (MQ=255)
                                                                                                                                          tAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGa          <  1:685938/149‑1 (MQ=255)
                                                                                                                                           aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGaa         <  4:175931/149‑1 (MQ=255)
                                                                                                                                            aCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTCATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAAt        <  3:834491/149‑1 (MQ=255)
                                                                                                                                            aCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAAt        <  2:876555/149‑1 (MQ=255)
                                                                                                                                              gccgcATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCTATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAAtag      >  3:598144/1‑147 (MQ=255)
                                                                                                                                              gccgcATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATgg      >  4:812439/1‑149 (MQ=255)
                                                                                                                                              gccgcATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATgg      >  3:865142/1‑149 (MQ=255)
                                                                                                                                               ccgcATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGt     <  3:706107/149‑1 (MQ=255)
                                                                                                                                                  cATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGtata  <  1:758888/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATA  >  NC_000913/2592858‑2593152

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