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

GTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGA  >  NC_000913/4612189‑4612479
                                                                                                                                              |                                                                                                                                                     
gtggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATAGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTg                                                                                                                                                <  3:602069/149‑1 (MQ=255)
 tggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTgg                                                                                                                                               <  1:917990/149‑1 (MQ=255)
 tggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATAGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTgg                                                                                                                                               <  4:64686/149‑1 (MQ=255)
  ggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGc                                                                                                                                              <  2:847422/149‑1 (MQ=255)
  ggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGc                                                                                                                                              <  2:422830/149‑1 (MQ=255)
  ggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGc                                                                                                                                              <  2:382640/149‑1 (MQ=255)
     cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCttt                                                                                                                                           >  3:832571/1‑149 (MQ=255)
          cGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCgg                                                                                                                                      <  4:832571/149‑1 (MQ=255)
            gTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTAGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGtt                                                                                                                                    <  2:231673/149‑1 (MQ=255)
               ccGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GCACTGGCTTTTGCGGTTTTa                                                                                                                                 <  2:182813/149‑1 (MQ=255)
                cGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAg                                                                                                                                >  4:672252/1‑149 (MQ=255)
                 gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGc                                                                                                                               >  3:389682/1‑149 (MQ=255)
                 gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGc                                                                                                                               >  1:786236/1‑149 (MQ=255)
                  cTGCAGGCGGAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTACAGC‑GGACTGGCTTTTGCGGTTTTAGCg                                                                                                                              >  4:795574/1‑149 (MQ=255)
                     cAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGc                                                                                                                                              <  4:859425/130‑1 (MQ=255)
                     cAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGc                                                                                                                                              >  3:859425/1‑130 (MQ=255)
                             gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAAtt                                                                                                                   >  3:572190/1‑149 (MQ=255)
                              cTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGa                                                                                                                             >  4:869923/1‑138 (MQ=255)
                              cTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGa                                                                                                                             <  3:869923/138‑1 (MQ=255)
                              cTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTg                                                                                                                  >  3:888846/1‑149 (MQ=255)
                               tAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                                 <  4:305026/149‑1 (MQ=255)
                                 aaaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGt                                                                                                               <  2:665318/149‑1 (MQ=255)
                                   aTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTa                                                                                                             >  1:158254/1‑149 (MQ=255)
                                      gTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACtt                                                                                                          >  1:53395/1‑149 (MQ=255)
                                       tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCttt                                                                                                                                           >  2:6625/1‑115 (MQ=255)
                                       tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACttt                                                                                                         >  4:915339/1‑149 (MQ=255)
                                       tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACttt                                                                                                         >  2:762226/1‑149 (MQ=255)
                                       tCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCC‑GCTGGACTGGCttt                                                                                                                                           <  1:6625/115‑1 (MQ=255)
                                         tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGaataat                                                                                                                        >  3:853270/1‑132 (MQ=255)
                                         tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGaataat                                                                                                                        <  4:853270/132‑1 (MQ=255)
                                         tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                                       >  4:445529/1‑149 (MQ=255)
                                          ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTc                                                                                                      >  1:797678/1‑149 (MQ=255)
                                                tcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTa                                                                                                                                 <  1:192655/116‑1 (MQ=255)
                                                tcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTa                                                                                                                                 >  2:192655/1‑116 (MQ=255)
                                                   ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGttt                                                                                                                                   >  4:724506/1‑111 (MQ=255)
                                                   ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGttt                                                                                                                                   <  3:724506/111‑1 (MQ=255)
                                                      ttattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaa                                                                                          >  2:193784/1‑149 (MQ=255)
                                                       tattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaac                                                                                         >  4:22730/1‑149 (MQ=255)
                                                          tCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTca                                                                                             >  4:682453/1‑142 (MQ=255)
                                                          tCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTca                                                                                             <  3:682453/142‑1 (MQ=255)
                                                          tCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaa                                                                                      >  2:765563/1‑149 (MQ=255)
                                                            tGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                    >  1:733859/1‑149 (MQ=255)
                                                              ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GAACTCGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                                  >  2:402217/1‑149 (MQ=255)
                                                                 cTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGa                                                                               >  1:896059/1‑149 (MQ=255)
                                                                        gTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                                                        <  1:765563/149‑1 (MQ=255)
                                                                        gTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                                                        <  4:888846/149‑1 (MQ=255)
                                                                         tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACGAAAACAGGAAGGCagag                                                                       <  4:286238/149‑1 (MQ=255)
                                                                         tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                       >  1:333717/1‑149 (MQ=255)
                                                                         tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                       <  3:797755/149‑1 (MQ=255)
                                                                             ttCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGAGTTACTTTGTCATTCTTTCACAACAAAAACAGTAAGGCAGAGGTgg                                                                   <  3:795574/149‑1 (MQ=255)
                                                                              tCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTggg                                                                  <  2:68367/149‑1 (MQ=255)
                                                                                       gAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAAt                                                         <  1:313289/149‑1 (MQ=255)
                                                                                        aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                                        >  4:924658/1‑149 (MQ=255)
                                                                                         aaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAc                                                       <  3:672252/149‑1 (MQ=255)
                                                                                          aaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAcc                                                      >  4:582488/1‑149 (MQ=255)
                                                                                           aCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCt                                                     >  4:892400/1‑149 (MQ=255)
                                                                                           aCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCt                                                     <  3:443733/149‑1 (MQ=255)
                                                                                             gACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGt                                                   >  1:119284/1‑149 (MQ=255)
                                                                                                 cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTg                                                                    <  1:410343/128‑1 (MQ=255)
                                                                                                 cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTg                                                                    >  2:410343/1‑128 (MQ=255)
                                                                                                 cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACg                                               >  4:749309/1‑149 (MQ=255)
                                                                                                 cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACg                                               >  4:399369/1‑149 (MQ=255)
                                                                                                     aTTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                           >  3:215573/1‑149 (MQ=255)
                                                                                                       ttCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGt                                         <  4:389682/149‑1 (MQ=255)
                                                                                                       ttCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGt                                         <  2:786236/149‑1 (MQ=255)
                                                                                                           gACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                    >  3:443803/1‑102 (MQ=255)
                                                                                                           gACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                    <  4:443803/102‑1 (MQ=255)
                                                                                                            aCCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                             >  3:477959/1‑108 (MQ=255)
                                                                                                            aCCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                             <  4:477959/108‑1 (MQ=255)
                                                                                                               aTTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgc                                 >  2:167176/1‑149 (MQ=255)
                                                                                                               aTTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgc                                 >  3:43693/1‑149 (MQ=255)
                                                                                                               aTTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgc                                 <  3:399369/149‑1 (MQ=255)
                                                                                                               aTTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTCTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                             >  1:483477/1‑105 (MQ=255)
                                                                                                               aTTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTCTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAg                                                                             <  2:483477/105‑1 (MQ=255)
                                                                                                                ttAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcg                                >  4:90844/1‑149 (MQ=255)
                                                                                                                 tAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcgc                               <  3:445529/149‑1 (MQ=255)
                                                                                                                       aaTTTTATCCAAAGCCAGGCCAGA‑GGACTGGCTTTTGCGGTTTTAGCGAATAAAAAATTGCCATACTTCGTCGTTCTTTCACAACAAAAACAGTAAGACAGAGTTGGAACAGCGAATACATCTAACTCAAGGTAACATTGCGCCGGTga                         >  4:140997/1‑148 (MQ=255)
                                                                                                                          tttATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGc                      >  3:290689/1‑149 (MQ=255)
                                                                                                                              tCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAg                  >  3:300095/1‑149 (MQ=255)
                                                                                                                               ccAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGc                                                             >  2:547533/1‑105 (MQ=255)
                                                                                                                               ccAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGc                                                             <  1:547533/105‑1 (MQ=255)
                                                                                                                               caaaaGCCAGTCC‑GCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAgg                 >  2:923714/3‑149 (MQ=255)
                                                                                                                                 aaaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                           <  3:80758/121‑1 (MQ=255)
                                                                                                                                 aaaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                           >  4:80758/1‑121 (MQ=255)
                                                                                                                                  aaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg              >  1:426821/1‑149 (MQ=255)
                                                                                                                                   aGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGACGTTGTCGCTAAGGCCGt             >  2:263738/1‑149 (MQ=255)
                                                                                                                                   aGCCAGTCCAGC‑GGACTGGCTGTTGCGGTTGTAGCGCATAATAAATTGCGGTACTTCGTCATTCTTTCACAACAAAAACAGGTAGCCAGAGGCGGGGCAGCGGATACCTGTAACGCTCCGTAATATTGCGCCGTTGGCGCTAAAGCCgc             >  2:893442/1‑148 (MQ=255)
                                                                                                                                    gCCAGTGCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGtt            >  2:743924/1‑149 (MQ=255)
                                                                                                                                     ccAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTc           >  2:311062/1‑149 (MQ=255)
                                                                                                                                      cAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTcc          >  3:72503/1‑149 (MQ=255)
                                                                                                                                            aGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTg    <  2:896059/149‑1 (MQ=255)
                                                                                                                                              c‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGa  <  1:244194/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGA  >  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: ‑