Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,612,331 C→T 7.4% intergenic (+42/‑80) ytjA → / → yjjU DUF1328 domain‑containing protein YtjA/putative patatin‑like phospholipase YjjU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,612,3310CT 7.4% 495.8 / 18.5 161intergenic (+42/‑80)ytjA/yjjUDUF1328 domain‑containing protein YtjA/putative patatin‑like phospholipase YjjU
Reads supporting (aligned to +/- strand):  ref base C (75/74);  new base T (8/4);  total (83/78)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 162 total aligned reads displayed.

GGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCT  >  NC_000913/4612182‑4612476
                                                                                                                                                     |                                                                                                                                                 
    ttGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGAc                                                                                                                                               <  2:964889/150‑1 (MQ=255)
       gtgggCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGGCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTgg                                                                                                                                            <  2:1735386/150‑1 (MQ=255)
         ggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCt                                                                                                                                          <  1:1774984/150‑1 (MQ=255)
             tGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCTGGACTGGCtttgg                                                                                                                                      <  2:270811/150‑3 (MQ=255)
              ggCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGc                                                                                                                                     >  2:245857/1‑150 (MQ=255)
              ggCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCAATAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGc                                                                                                                                     >  2:1261270/1‑150 (MQ=255)
               gcgggTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAATTTTTATCCAAAGCCAGTCCGCTGGACTGGCtttggat                                                                                                                                    <  2:1490333/147‑5 (MQ=255)
               gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCg                                                                                                                                    <  2:580542/150‑1 (MQ=255)
                ccGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCgg                                                                                                                                   <  1:2473214/150‑1 (MQ=255)
                ccGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCgg                                                                                                                                   <  2:461112/150‑1 (MQ=255)
                 cGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGt                                                                                                                                  <  1:1461195/150‑1 (MQ=255)
                  ggTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGtt                                                                                                                                 >  1:2271023/1‑150 (MQ=255)
                    tACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGtttt                                                                                                                               <  2:707393/150‑1 (MQ=255)
                      ccGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAg                                                                                                                             <  2:1270484/150‑1 (MQ=255)
                       cGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGc                                                                                                                            >  1:1410945/1‑150 (MQ=255)
                           gCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaat                                                                                                                        >  2:2075472/1‑150 (MQ=255)
                              ggCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGCCCCTAGATTTCCGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGCCTGGCTTTTGCGGTTTTAGCGaattat                                                                                                                     >  2:1158407/1‑150 (MQ=255)
                              ggCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaataat                                                                                                                     <  2:2378490/150‑1 (MQ=255)
                                 gcAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATaaa                                                                                                                  >  1:1180973/1‑150 (MQ=255)
                                      tAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCg                                                                                                             <  2:1417922/150‑1 (MQ=255)
                                       aaaaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGt                                                                                                            >  2:456793/1‑150 (MQ=255)
                                          aTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTAc                                                                                                         >  1:995063/1‑150 (MQ=255)
                                            tGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACtt                                                                                                       <  2:1823801/150‑1 (MQ=255)
                                              tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTg                                                                                                     <  2:1693156/150‑1 (MQ=255)
                                                tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTc                                                                                                   <  2:1041945/150‑1 (MQ=255)
                                                    tcgtcgGGTTTATTATGTTCCTGGTTTTTTTGTTCATGTGCCGAAAACGACCCTAGATTTCAGTCCATTAAGTAAATTTTATCCAAAGCCATTCCAGCGGACTGGCTTTTGAGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTc                                                                                               <  2:157187/150‑1 (MQ=255)
                                                     cgtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCt                                                                                              >  1:1675756/1‑150 (MQ=255)
                                                     cgtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCt                                                                                              <  1:1262006/150‑1 (MQ=255)
                                                       tcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCttt                                                                                            <  1:1094187/150‑1 (MQ=255)
                                                        cgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                                           <  1:2366924/150‑1 (MQ=255)
                                                        cgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                                           >  2:59780/1‑150 (MQ=255)
                                                          ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcac                                                                                         <  2:886316/150‑1 (MQ=255)
                                                          ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcac                                                                                         <  2:1519300/150‑1 (MQ=255)
                                                          ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaa                                                                                         >  2:2283713/1‑149 (MQ=255)
                                                           gATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaca                                                                                        >  2:197701/1‑150 (MQ=255)
                                                            attattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaa                                                                                       >  1:1914637/1‑150 (MQ=255)
                                                              tattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaaca                                                                                     <  2:184042/150‑1 (MQ=255)
                                                               attCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCAcaacaa                                                                                    >  2:875516/1‑150 (MQ=255)
                                                                ttCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaa                                                                                   >  2:1779312/1‑150 (MQ=255)
                                                                   tGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAAc                                                                                <  1:926088/150‑1 (MQ=255)
                                                                      tCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAgg                                                                             <  2:107241/150‑1 (MQ=255)
                                                                         tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGACAGTACAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAATCAGTAAg                                                                          <  2:1360530/150‑1 (MQ=255)
                                                                          ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAgg                                                                         <  1:88656/150‑1 (MQ=255)
                                                                          ggTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTGCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAgg                                                                         <  2:1109128/150‑1 (MQ=255)
                                                                              aGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                                                     <  2:586472/150‑1 (MQ=255)
                                                                              aGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                                                     <  1:1030403/150‑1 (MQ=255)
                                                                               gTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                    >  1:1716771/1‑150 (MQ=255)
                                                                               gTTTGTTAATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCAAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                    >  2:2130912/1‑150 (MQ=255)
                                                                                tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAgg                                                                   <  1:2376083/150‑1 (MQ=255)
                                                                                tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAgg                                                                   >  2:681371/1‑150 (MQ=255)
                                                                                    ttCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTggg                                                               >  2:433125/1‑150 (MQ=255)
                                                                                      cATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGc                                                             <  2:2042196/150‑1 (MQ=255)
                                                                                         gggCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGc                                                          <  2:1370462/150‑1 (MQ=255)
                                                                                          ggCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCg                                                         <  2:821787/150‑1 (MQ=255)
                                                                                          ggCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCg                                                         <  1:2197788/150‑1 (MQ=255)
                                                                                          ggCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCg                                                         >  1:1620355/1‑150 (MQ=255)
                                                                                           gCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGa                                                        >  2:734049/1‑150 (MQ=255)
                                                                                           gCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGa                                                        >  2:402673/1‑150 (MQ=255)
                                                                                              gAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                                     >  1:1319158/1‑150 (MQ=255)
                                                                                              gAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                                     <  1:1265891/150‑1 (MQ=255)
                                                                                               aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAc                                                    >  1:196832/1‑150 (MQ=255)
                                                                                                aaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAcc                                                   <  1:2075472/150‑1 (MQ=255)
                                                                                                aaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAcc                                                   <  1:1222337/150‑1 (MQ=255)
                                                                                                 aaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCt                                                  <  2:1866996/150‑1 (MQ=255)
                                                                                                        cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGc                                           <  1:402673/150‑1 (MQ=255)
                                                                                                        cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGc                                           >  2:1796661/1‑150 (MQ=255)
                                                                                                            aTTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTgt                                       <  1:2120844/150‑2 (MQ=255)
                                                                                                             tttCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGt                                      <  2:1319158/150‑1 (MQ=255)
                                                                                                              ttCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGTGAATACCTGTAACGCTTGGTa                                     <  1:1472688/150‑1 (MQ=255)
                                                                                                                cAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCAGAGTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGTTAat                                   >  2:1922769/1‑150 (MQ=255)
                                                                                                                  gACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAatat                                 >  1:2263075/1‑150 (MQ=255)
                                                                                                                  gACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAatat                                 >  2:2370776/1‑150 (MQ=255)
                                                                                                                   aCCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATAtt                                >  2:632136/1‑150 (MQ=255)
                                                                                                                       ttAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcgc                            <  1:1838817/150‑1 (MQ=255)
                                                                                                                        tAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGcc                           >  1:2412313/1‑150 (MQ=255)
                                                                                                                          aGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGt                         >  2:125869/1‑150 (MQ=255)
                                                                                                                           gAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGtt                        >  2:550202/1‑150 (MQ=255)
                                                                                                                           gAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGtt                        >  1:2408127/1‑150 (MQ=255)
                                                                                                                            aaaaTTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTg                       <  1:61518/150‑1 (MQ=255)
                                                                                                                                ttttATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGc                   <  2:1410945/150‑1 (MQ=255)
                                                                                                                                 tttATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCTGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCt                  >  2:219260/1‑150 (MQ=255)
                                                                                                                                 aaccgcaaaaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCt                  >  1:57067/8‑150 (MQ=255)
                                                                                                                                   ccgcaaaaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTGGGTAATATTGCGCCGTTGTCGCTaa                >  1:243049/6‑150 (MQ=255)
                                                                                                                                    cgcaaaaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAATTTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGGGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAg               >  2:1211593/5‑150 (MQ=255)
                                                                                                                                    cgcaaaaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAg               >  1:2047303/5‑150 (MQ=255)
                                                                                                                                    aTCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAg               <  1:1265188/150‑1 (MQ=255)
                                                                                                                                     tCCAAAGCCAGTCCAGCGGCCTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAgg              >  1:800823/1‑150 (MQ=255)
                                                                                                                                     gcaaaaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAgg              >  1:1783600/4‑150 (MQ=255)
                                                                                                                                      ccAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTGTGTCATTCTTTCACAACAAAAACAGGATGGCAGAGGTGGGGCAGCGAATACCTGTAAGGCTTGGTAATATTGCGCCGTTGTCACTAAGGc             >  2:2013024/1‑150 (MQ=255)
                                                                                                                                       cAAAGCCAGTCCAGCGGACTGGCTTTTGTGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc            >  1:2575/1‑150 (MQ=255)
                                                                                                                                       cAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc            >  1:684838/1‑150 (MQ=255)
                                                                                                                                       cAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc            >  1:480844/1‑150 (MQ=255)
                                                                                                                                       cAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc            >  2:127643/1‑150 (MQ=255)
                                                                                                                                         aaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGt          >  1:232882/1‑150 (MQ=255)
                                                                                                                                         aaGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGt          <  2:196832/150‑1 (MQ=255)
                                                                                                                                          aGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGtt         >  1:929800/1‑150 (MQ=255)
                                                                                                                                            ccAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTATCACAACAAAAACAGGAAGGCAGACGTGGCGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTcc       >  1:800819/1‑150 (MQ=255)
                                                                                                                                             cAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCa      <  1:632136/150‑1 (MQ=255)
                                                                                                                                               gTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGc    >  2:2453779/1‑150 (MQ=255)
                                                                                                                                                tCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGcc   <  2:2206504/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCT  >  NC_000913/4612182‑4612476

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