Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,612,331 C→T 5.5% intergenic (+42/‑80) ytjA → / → yjjU uncharacterized protein/putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,612,3310CT 5.5% 293.0 / 5.0 91intergenic (+42/‑80)ytjA/yjjUuncharacterized protein/putative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base C (49/36);  new base T (3/2);  total (53/38)
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.54e-01

GTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGA  >  NC_000913/4612189‑4612479
                                                                                                                                              |                                                                                                                                                     
gtggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATAGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTg                                                                                                                                                <  3:602069/149‑1 (MQ=255)
 tggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTgg                                                                                                                                               <  1:917990/149‑1 (MQ=255)
 tggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATAGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTgg                                                                                                                                               <  4:64686/149‑1 (MQ=255)
  ggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCTGGACTGGc                                                                                                                                              <  2:847422/149‑1 (MQ=255)
  ggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGc                                                                                                                                              <  2:422830/149‑1 (MQ=255)
  ggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGc                                                                                                                                              <  2:382640/149‑1 (MQ=255)
     cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCttt                                                                                                                                           >  3:832571/1‑149 (MQ=255)
          cGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCgg                                                                                                                                      <  4:832571/149‑1 (MQ=255)
            gTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTAGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGtt                                                                                                                                    <  2:231673/149‑1 (MQ=255)
               ccGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGCACTGGCTTTTGCGGTTTTa                                                                                                                                 <  2:182813/149‑1 (MQ=255)
                cGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAg                                                                                                                                >  4:672252/1‑149 (MQ=255)
                 gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGc                                                                                                                               >  3:389682/1‑149 (MQ=255)
                 gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGc                                                                                                                               >  1:786236/1‑149 (MQ=255)
                  cTGCAGGCGGAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTACAGCGGACTGGCTTTTGCGGTTTTAGCg                                                                                                                              >  4:795574/1‑149 (MQ=255)
                     cAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGc                                                                                                                                              <  4:859425/130‑1 (MQ=255)
                     cAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGc                                                                                                                                              >  3:859425/1‑130 (MQ=255)
                             gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAAtt                                                                                                                   >  3:572190/1‑149 (MQ=255)
                              cTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGa                                                                                                                             >  4:869923/1‑138 (MQ=255)
                              cTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGa                                                                                                                             <  3:869923/138‑1 (MQ=255)
                              cTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTg                                                                                                                  >  3:888846/1‑149 (MQ=255)
                               tAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                                 <  4:305026/149‑1 (MQ=255)
                                 aaaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGt                                                                                                               <  2:665318/149‑1 (MQ=255)
                                   aTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTa                                                                                                             >  1:158254/1‑149 (MQ=255)
                                      gTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACtt                                                                                                          >  1:53395/1‑149 (MQ=255)
                                       tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCTGGACTGGCttt                                                                                                                                           >  2:6625/1‑115 (MQ=255)
                                       tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACttt                                                                                                         >  4:915339/1‑149 (MQ=255)
                                       tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACttt                                                                                                         >  2:762226/1‑149 (MQ=255)
                                       tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCTGGACTGGCttt                                                                                                                                           <  1:6625/115‑1 (MQ=255)
                                         tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaataat                                                                                                                        >  3:853270/1‑132 (MQ=255)
                                         tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaataat                                                                                                                        <  4:853270/132‑1 (MQ=255)
                                         tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                                       >  4:445529/1‑149 (MQ=255)
                                          ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTc                                                                                                      >  1:797678/1‑149 (MQ=255)
                                                tcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTa                                                                                                                                 <  1:192655/116‑1 (MQ=255)
                                                tcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTa                                                                                                                                 >  2:192655/1‑116 (MQ=255)
                                                   ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGttt                                                                                                                                   >  4:724506/1‑111 (MQ=255)
                                                   ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGttt                                                                                                                                   <  3:724506/111‑1 (MQ=255)
                                                      ttattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaa                                                                                          >  2:193784/1‑149 (MQ=255)
                                                       tattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaac                                                                                         >  4:22730/1‑149 (MQ=255)
                                                          tCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTca                                                                                             >  4:682453/1‑142 (MQ=255)
                                                          tCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTca                                                                                             <  3:682453/142‑1 (MQ=255)
                                                          tCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaa                                                                                      >  2:765563/1‑149 (MQ=255)
                                                            tGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                    >  1:733859/1‑149 (MQ=255)
                                                              ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGAACTCGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                                  >  2:402217/1‑149 (MQ=255)
                                                                 cTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGa                                                                               >  1:896059/1‑149 (MQ=255)
                                                                        gTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                                                        <  1:765563/149‑1 (MQ=255)
                                                                        gTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                                                        <  4:888846/149‑1 (MQ=255)
                                                                         tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACGAAAACAGGAAGGCagag                                                                       <  4:286238/149‑1 (MQ=255)
                                                                         tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                       >  1:333717/1‑149 (MQ=255)
                                                                         tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                       <  3:797755/149‑1 (MQ=255)
                                                                             ttCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGAGTTACTTTGTCATTCTTTCACAACAAAAACAGTAAGGCAGAGGTgg                                                                   <  3:795574/149‑1 (MQ=255)
                                                                              tCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTggg                                                                  <  2:68367/149‑1 (MQ=255)
                                                                                       gAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAAt                                                         <  1:313289/149‑1 (MQ=255)
                                                                                        aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                                        >  4:924658/1‑149 (MQ=255)
                                                                                         aaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAc                                                       <  3:672252/149‑1 (MQ=255)
                                                                                          aaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAcc                                                      >  4:582488/1‑149 (MQ=255)
                                                                                           aCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCt                                                     >  4:892400/1‑149 (MQ=255)
                                                                                           aCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCt                                                     <  3:443733/149‑1 (MQ=255)
                                                                                             gACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGt                                                   >  1:119284/1‑149 (MQ=255)
                                                                                                 cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTg                                                                    <  1:410343/128‑1 (MQ=255)
                                                                                                 cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTg                                                                    >  2:410343/1‑128 (MQ=255)
                                                                                                 cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACg                                               >  4:749309/1‑149 (MQ=255)
                                                                                                 cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACg                                               >  4:399369/1‑149 (MQ=255)
                                                                                                     aTTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                           >  3:215573/1‑149 (MQ=255)
                                                                                                       ttCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGt                                         <  4:389682/149‑1 (MQ=255)
                                                                                                       ttCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGt                                         <  2:786236/149‑1 (MQ=255)
                                                                                                           gACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                    >  3:443803/1‑102 (MQ=255)
                                                                                                           gACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                    <  4:443803/102‑1 (MQ=255)
                                                                                                            aCCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                             >  3:477959/1‑108 (MQ=255)
                                                                                                            aCCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                             <  4:477959/108‑1 (MQ=255)
                                                                                                               aTTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgc                                 >  2:167176/1‑149 (MQ=255)
                                                                                                               aTTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgc                                 >  3:43693/1‑149 (MQ=255)
                                                                                                               aTTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgc                                 <  3:399369/149‑1 (MQ=255)
                                                                                                               aTTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTCTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                             >  1:483477/1‑105 (MQ=255)
                                                                                                               aTTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTCTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                             <  2:483477/105‑1 (MQ=255)
                                                                                                                ttAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcg                                >  4:90844/1‑149 (MQ=255)
                                                                                                                 tAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcgc                               <  3:445529/149‑1 (MQ=255)
                                                                                                                       aaTTTTATCCAAAGCCAGGCCAGAGGACTGGCTTTTGCGGTTTTAGCGAATAAAAAATTGCCATACTTCGTCGTTCTTTCACAACAAAAACAGTAAGACAGAGTTGGAACAGCGAATACATCTAACTCAAGGTAACATTGCGCCGGTga                         >  4:140997/1‑148 (MQ=255)
                                                                                                                          tttATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGc                      >  3:290689/1‑149 (MQ=255)
                                                                                                                              tCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAg                  >  3:300095/1‑149 (MQ=255)
                                                                                                                               ccAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGc                                                             >  2:547533/1‑105 (MQ=255)
                                                                                                                               ccAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGc                                                             <  1:547533/105‑1 (MQ=255)
                                                                                                                               caaaaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAgg                 >  2:923714/3‑149 (MQ=255)
                                                                                                                                 aaaGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                           <  3:80758/121‑1 (MQ=255)
                                                                                                                                 aaaGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                           >  4:80758/1‑121 (MQ=255)
                                                                                                                                  aaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg              >  1:426821/1‑149 (MQ=255)
                                                                                                                                   aGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGACGTTGTCGCTAAGGCCGt             >  2:263738/1‑149 (MQ=255)
                                                                                                                                   aGCCAGTCCAGCGGACTGGCTGTTGCGGTTGTAGCGCATAATAAATTGCGGTACTTCGTCATTCTTTCACAACAAAAACAGGTAGCCAGAGGCGGGGCAGCGGATACCTGTAACGCTCCGTAATATTGCGCCGTTGGCGCTAAAGCCgc             >  2:893442/1‑148 (MQ=255)
                                                                                                                                    gCCAGTGCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGtt            >  2:743924/1‑149 (MQ=255)
                                                                                                                                     ccAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTc           >  2:311062/1‑149 (MQ=255)
                                                                                                                                      cAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTcc          >  3:72503/1‑149 (MQ=255)
                                                                                                                                            aGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTg    <  2:896059/149‑1 (MQ=255)
                                                                                                                                              cGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGa  <  1:244194/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGA  >  NC_000913/4612189‑4612479

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