Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,840,777 T→C 6.0% intergenic (+3/+6) ynjE → / ← ynjF molybdopterin synthase sulfurtransferase/CDP‑alcohol phosphatidyltransferase family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,840,7770TC 6.0% 206.7 / 3.1 67intergenic (+3/+6)ynjE/ynjFmolybdopterin synthase sulfurtransferase/CDP‑alcohol phosphatidyltransferase family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base C (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

GCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCG  >  NC_000913/1840629‑1840925
                                                                                                                                                    |                                                                                                                                                    
gCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACAt                                                                                                                                                      <  6:989/149‑1 (MQ=255)
 cTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACAtt                                                                                                                                                     <  5:166392/149‑1 (MQ=255)
    gcgcgcGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGc                                                                                                                                                   >  8:215951/1‑148 (MQ=255)
       cgcgTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAAt                                                                                                                                                >  5:375071/1‑148 (MQ=255)
        gcgTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTa                                                                                                                                              >  2:175037/1‑149 (MQ=255)
              gAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACg                                                                                                                                        <  8:117738/149‑1 (MQ=255)
                aaCCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCt                                                                                                                                      >  7:232647/1‑149 (MQ=255)
                  ccTTTATGTACGCACGCGCCATGGGTTGGAAGAATTTTTCCTTGTATGACGGCGGCTGGTACGAATGGAGCAGCTATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCAGGACAGTAGCAAATAACATTGCAATTACTGACGCTgg                                                                                                                                    <  3:30692/149‑1 (MQ=255)
                    tttATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAg                                                                                                                                  >  8:355445/1‑149 (MQ=255)
                    tttATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACACTGCAATTACTGACGCTGGAg                                                                                                                                  >  8:87346/1‑149 (MQ=255)
                    tttATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACACTGCAATTACTGACGCTGGAg                                                                                                                                  >  4:93683/1‑149 (MQ=255)
                      tATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGc                                                                                                                                 <  2:210877/148‑1 (MQ=255)
                         gTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACt                                                                                                                             >  1:53821/1‑149 (MQ=255)
                         gTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACt                                                                                                                             >  4:215169/1‑149 (MQ=255)
                            cgcacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCa                                                                                                                          <  5:119980/149‑1 (MQ=255)
                            cgcacgcGCCATGGGGTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATCGAGCATCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGTCCCTGACAGTAGCAAATTACATTGCAATTACTGACGCTGGAGCGACttg                                                                                                                           <  4:66702/148‑2 (MQ=255)
                             gcacgcGCCAGGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAg                                                                                                                         >  3:151634/1‑149 (MQ=255)
                              cacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                                                                                                                        <  1:241465/149‑1 (MQ=255)
                                       tGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAc                                                                                                               <  8:252511/149‑1 (MQ=255)
                                          gTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCa                                                                                                            >  7:14271/1‑149 (MQ=255)
                                              gaagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTcc                                                                                                        <  2:53821/149‑1 (MQ=255)
                                                agaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcac                                                                                                      >  4:208674/1‑149 (MQ=255)
                                                 gaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcaca                                                                                                     <  3:33169/149‑1 (MQ=255)
                                                 gaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcaca                                                                                                     <  7:36627/149‑1 (MQ=255)
                                                    tGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACAccc                                                                                                  >  8:13543/1‑149 (MQ=255)
                                                     gTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcg                                                                                                 >  3:221824/1‑149 (MQ=255)
                                                      tttCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgc                                                                                                <  2:335960/149‑1 (MQ=255)
                                                      tttCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACACTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcg                                                                                                 <  3:93683/148‑1 (MQ=255)
                                                      tttCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACACTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcg                                                                                                 <  7:87346/148‑1 (MQ=255)
                                                       ttCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgcg                                                                                               >  4:124901/1‑149 (MQ=255)
                                                        tCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCAAAAAAATCCGGTAGCAAACGGTGAACGCGGCCCGGACAGTAGCAAATTACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCAACTCCACACCcgcg                                                                                               >  5:208283/1‑148 (MQ=255)
                                                          cGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCTGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTa                                                                                             >  7:291035/1‑148 (MQ=255)
                                                          cGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTa                                                                                             >  2:94783/1‑148 (MQ=255)
                                                           gtgtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTaaa                                                                                           <  8:232647/149‑1 (MQ=255)
                                                           gtgtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTaaa                                                                                           <  7:215951/149‑1 (MQ=255)
                                                            tgtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAg                                                                                          <  6:336330/149‑1 (MQ=255)
                                                               aTGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtc                                                                                       <  6:223279/149‑1 (MQ=255)
                                                                tGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcg                                                                                      <  7:159307/149‑1 (MQ=255)
                                                                tGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtc                                                                                       <  6:228980/148‑1 (MQ=255)
                                                                 gACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTaaa                                                                                           <  1:13222/143‑1 (MQ=255)
                                                                 gACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTaaa                                                                                           >  2:13222/1‑143 (MQ=255)
                                                                 gACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcgt                                                                                     >  6:148180/1‑149 (MQ=255)
                                                                     gcggcTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATAACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATc                                                                                 >  7:194965/1‑149 (MQ=255)
                                                                           ggTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGca                                                                            >  1:178605/1‑148 (MQ=255)
                                                                                gAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGc                                                                      >  5:298887/1‑149 (MQ=255)
                                                                                  aTGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCAc                                                                    <  3:124901/149‑1 (MQ=255)
                                                                                        cagcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtc                                                                                       <  4:40629/124‑1 (MQ=255)
                                                                                        cagcGATCCAAAAAATCCGGTAGCAACAGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtc                                                                                       >  3:40629/1‑124 (MQ=255)
                                                                                          gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGGCCGGACAGTAGCTAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACCCCCGCGTAAAGGTGGTCATCCAGCCCAACGCACCGAAGCTc                                                            >  6:91304/1‑149 (MQ=255)
                                                                                           cGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATcc                                                           <  6:375071/149‑1 (MQ=255)
                                                                                           cGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATcc                                                           <  5:148180/149‑1 (MQ=255)
                                                                                             aTCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccac                                                         >  3:147099/1‑149 (MQ=255)
                                                                                             aTCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccac                                                         <  2:178605/149‑1 (MQ=255)
                                                                                                cAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATcc                                                           >  3:121959/1‑144 (MQ=255)
                                                                                                cAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATcc                                                           <  4:121959/144‑1 (MQ=255)
                                                                                                  aaaaaTCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaa                                                    >  8:208073/1‑149 (MQ=255)
                                                                                                       tCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaaccacg                                               <  3:143773/149‑1 (MQ=255)
                                                                                                                   cGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAg                                   <  7:355445/149‑1 (MQ=255)
                                                                                                                   cGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAg                                   <  3:208674/149‑1 (MQ=255)
                                                                                                                   cGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAg                                   <  4:221824/149‑1 (MQ=255)
                                                                                                                   cGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTAGTCATCCAGCACAACGCACCGAAGATCCACTCAAACCACGCAAACCATGCAg                                   <  8:14271/149‑1 (MQ=255)
                                                                                                                    ggTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAATGTCGTCATCCAGCACAACGCACCGAAGATCCACTCAAACCACGCAAACCATGCAgg                                  <  4:151634/149‑1 (MQ=255)
                                                                                                                         aCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCacgcaaaccacgcaaacc                                         <  6:287621/137‑1 (MQ=255)
                                                                                                                         aCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCacgcaaaccacgcaaacc                                         >  5:287621/1‑137 (MQ=255)
                                                                                                                                 cGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGCCTTCAGCGTCAGCTACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCCCGCAAACCACGCAAACCATGCAGGGAATAAACATcc                      >  6:121835/1‑148 (MQ=255)
                                                                                                                                 cGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATcc                      >  2:380574/1‑148 (MQ=255)
                                                                                                                                        tAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCCCACCCGCGTAAAGGTCGTCATCCCGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGCATAAACATCCCAGcaca               >  5:372372/1‑148 (MQ=255)
                                                                                                                                             aaTAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGt         <  6:298887/149‑1 (MQ=255)
                                                                                                                                                  cATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGt         >  6:374292/1‑144 (MQ=255)
                                                                                                                                                    ttGCAATTACTGACGCTGGAGCGGCTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCg  >  8:390451/1‑149 (MQ=255)
                                                                                                                                                    ttGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCg  >  5:15253/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCG  >  NC_000913/1840629‑1840925

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