Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,083,705 G→A 5.5% T110T (ACG→ACA efeB → deferrrochelatase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,083,7050GA 5.5% 233.0 / 2.3 74T110T (ACG→ACAefeBdeferrrochelatase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (41/28);  new base A (2/2);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

AGGCAGGGATCCTGACGCCACAACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCCACAGATGCCAAAAAAGCTGCAGAAGA  >  NC_000913/1083557‑1083853
                                                                                                                                                    |                                                                                                                                                    
aggcaggGATCCTGACGCCACAACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAACCg                                                                                                                                                      >  1:376022/1‑149 (MQ=255)
  gcaggGATCCTGACGCCACAACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGcc                                                                                                                                                    >  1:1065131/1‑149 (MQ=255)
  gcaggGATCCTGACGCCACAACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGcc                                                                                                                                                    >  2:357892/1‑149 (MQ=255)
       gATCCTGACGCCACAACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATc                                                                                                                                               <  1:411201/149‑1 (MQ=255)
          ccTGACGCCACAACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACACAGGGCGGAGCAGCACCAGAAACGCCAAATCcgc                                                                                                                                            <  2:1431649/149‑1 (MQ=255)
          ccTCACGCCACAACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGAGCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCcgc                                                                                                                                            <  1:38078/149‑1 (MQ=255)
             gACGCCACAACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGcc                                                                                                                                         <  1:1221784/149‑1 (MQ=255)
                gCCACAACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGc                                                                                                                                      >  1:479806/1‑149 (MQ=255)
                  cacaACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGcca                                                                                                                                    <  1:34087/149‑1 (MQ=255)
                    caacaGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGGTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGccacc                                                                                                                                  >  2:1080149/1‑149 (MQ=255)
                     aacaGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGccacca                                                                                                                                 >  1:1105289/1‑149 (MQ=255)
                        aGGCCGCAATGATGCGGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCtctc                                                                                                                              <  1:218474/149‑1 (MQ=255)
                             gCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTc                                                                                                                         <  2:1453246/149‑1 (MQ=255)
                                       tggtggCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCtt                                                                                                               >  2:385306/1‑149 (MQ=255)
                                        ggtggCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTg                                                                                                              <  1:385306/149‑1 (MQ=255)
                                              gTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCt                                                                                                        >  1:761584/1‑149 (MQ=255)
                                               tttGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTa                                                                                                       >  1:527526/1‑149 (MQ=255)
                                                    tgtgCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTg                                                                                                  <  1:1523120/149‑1 (MQ=255)
                                                        cTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCCCCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGcc                                                                                              >  2:1546575/1‑149 (MQ=255)
                                                         ttGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGccc                                                                                             >  1:349275/1‑149 (MQ=255)
                                                         ttGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGccc                                                                                             >  2:391365/1‑149 (MQ=255)
                                                          tGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCg                                                                                            <  2:704318/149‑1 (MQ=255)
                                                           gCCAGCGATAAAGCCGATCTTGAGCGGTTGTTGCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGa                                                                                           <  1:316263/149‑1 (MQ=255)
                                                             cAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATa                                                                                         >  2:1163941/1‑149 (MQ=255)
                                                              aGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGCGAGGAGCAGCACCAGAAACGCCAAATCCTAGCCTGCCACCACTCGATTCCGGCATACTTGGCGGCTACATTGCGCCCGATaa                                                                                        >  1:55249/1‑149 (MQ=255)
                                                               gCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGAGTCTGGTTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCAAAATCCGCGCCTGCCAGCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAt                                                                                       <  2:264421/149‑1 (MQ=255)
                                                               gCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAt                                                                                       <  1:175758/149‑1 (MQ=255)
                                                                 gATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAtct                                                                                     >  2:1059536/1‑149 (MQ=255)
                                                                   tAAAGCCGATCTTGAGCTGTTGTTTCGCTTGTTGACTCATCGTTTTGCTTTTCTGACTCAGGGCGGTGCAGCACCAGATACGCCTAATCGGCGCCTGCGACCACTCGTTTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCa                                                                                   <  1:1080149/149‑1 (MQ=255)
                                                                         cGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCa                                                                             <  2:479806/149‑1 (MQ=255)
                                                                           aTCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACg                                                                           >  1:652362/1‑149 (MQ=255)
                                                                             cttgctcGGTTGGTTCCCTTGTTGACTCAGCGGTTTGGTTTTCTGTCTCCCGGGGGCCGAGCACCAGTAACCCCAAAACCGCCCCTTGCACCAATCGATTCCGGCATTTTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGtt                                                                         <  2:790304/143‑1 (MQ=255)
                                                                              ttGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCCCTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTa                                                                        >  2:230267/1‑149 (MQ=255)
                                                                              ttGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTa                                                                        >  2:411632/1‑149 (MQ=255)
                                                                                gAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCCCCATCCCGTTATc                                                                      >  1:971153/1‑149 (MQ=255)
                                                                                 aGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACACCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCg                                                                     >  2:724329/1‑149 (MQ=255)
                                                                                 aGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACACCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCg                                                                     >  2:730466/1‑149 (MQ=255)
                                                                                  gCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCgg                                                                    >  1:947690/1‑149 (MQ=255)
                                                                                   cGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTGCATTGCGCCCGATAATCTCACCATCACGTTATCGGt                                                                   <  2:23920/149‑1 (MQ=255)
                                                                                     gttgttTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGgggg                                                                 >  1:879393/1‑149 (MQ=255)
                                                                                     gttgttTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTgg                                                                 >  1:1186858/1‑149 (MQ=255)
                                                                                         tttCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCAACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCa                                                             >  1:1174330/1‑149 (MQ=255)
                                                                                                tgttgACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTg                                                      <  2:1030088/149‑1 (MQ=255)
                                                                                                 gttgACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGt                                                     <  1:1122956/149‑1 (MQ=255)
                                                                                                      cTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCTCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTtgat                                                <  2:1027330/149‑1 (MQ=255)
                                                                                                      cTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGt                                                               >  2:1444465/1‑134 (MQ=255)
                                                                                                      cTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGt                                                               <  1:1444465/134‑1 (MQ=255)
                                                                                                      cTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTtgat                                                <  1:600675/149‑1 (MQ=255)
                                                                                                      cTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTtgat                                                >  1:602500/1‑149 (MQ=255)
                                                                                                        cAGCGTTTTGCTTTTCTGCCTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTtgatga                                              <  1:1160977/149‑1 (MQ=255)
                                                                                                            gTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACAATCACGTTATCGGTGGGTCACTCATTGTTTGATGAgcgc                                          >  2:739019/1‑149 (MQ=255)
                                                                                                                tGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAATCGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTTCATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTg                                      <  1:538043/149‑1 (MQ=255)
                                                                                                                 gCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTgg                                     >  2:865208/1‑149 (MQ=255)
                                                                                                                   ttttCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCCCCACTCGATTCCGGCATTCTTGGCGCCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTGATTGTTTGATTAGCGCTTTGGca                                   >  2:79435/1‑148 (MQ=255)
                                                                                                                   ttttCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCCCGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGcc                                   >  1:875410/1‑149 (MQ=255)
                                                                                                                   ttttCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGcc                                   >  1:881168/1‑149 (MQ=255)
                                                                                                                       cTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCCCCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTgc                               >  1:9540/1‑149 (MQ=255)
                                                                                                                        tGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCNCNTGCCACCACTCGATTCCGGNATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTgcg                              >  1:1066848/1‑149 (MQ=255)
                                                                                                                        tGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTgcg                              >  1:1063112/1‑149 (MQ=255)
                                                                                                                        tGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTgcg                              >  1:1072405/1‑149 (MQ=255)
                                                                                                                         gACTCAGGGCGGAGCAGCACCAGAAACACCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTgcgc                             <  1:724329/149‑1 (MQ=255)
                                                                                                                         gACTCAGGGCGGAGCAGCACCAGAAACACCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTgcgc                             <  1:730466/149‑1 (MQ=255)
                                                                                                                          aCTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGcc                            >  2:1545149/1‑149 (MQ=255)
                                                                                                                          aCTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGcc                            >  1:397706/1‑149 (MQ=255)
                                                                                                                            tCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCcac                          >  1:699597/1‑149 (MQ=255)
                                                                                                                            tCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCcac                          >  1:693641/1‑149 (MQ=255)
                                                                                                                            tCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACGCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCTTTGTTTGATGATCGCTTTGGCCTTGCGCcac                          >  1:783475/1‑149 (MQ=255)
                                                                                                                             cAGGGCGGAGCAGCACCAGCACCGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGGGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCcaca                         <  2:1065131/149‑1 (MQ=255)
                                                                                                                               gggCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCTTTGTTTGATGAGCGCTTTGGCCTTGCGCCACAGa                       >  1:1325100/1‑149 (MQ=255)
                                                                                                                                        agcaCCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCCACAGATGCCaaaaa              <  2:376022/149‑1 (MQ=255)
                                                                                                                                        agcaCCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCCACAGATGCCaaaaa              <  1:1059536/149‑1 (MQ=255)
                                                                                                                                             cAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCCACAGATGCCAAAAAAGCTg         >  1:698757/1‑149 (MQ=255)
                                                                                                                                             cAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCCACAGATGCCAAAAAAGCTg         <  2:349275/149‑1 (MQ=255)
                                                                                                                                                aaaCGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCCACAGATGCCAAAAAAGCTGCag      <  2:15807/149‑1 (MQ=255)
                                                                                                                                                aaaCGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCCACAGATGCCAAAAAAGCTGCag      >  1:193355/1‑149 (MQ=255)
                                                                                                                                                    gCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCCACAGATGCCAAAAAAGCTGCagaaga  >  1:1241977/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGCAGGGATCCTGACGCCACAACAGGCCGCAATGATGCTGGTGGCGTTTGATGTGCTTGCCAGCGATAAAGCCGATCTTGAGCGGTTGTTTCGCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCCTTGCGCCACAGATGCCAAAAAAGCTGCAGAAGA  >  NC_000913/1083557‑1083853

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: