Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,987,175 C→G 5.5% M236I (ATG→ATC hemY ← putative protoheme IX synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,1750CG 5.5% 348.6 / 8.9 109M236I (ATG→ATChemYputative protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base C (46/57);  new base G (3/3);  total (49/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 112 total aligned reads displayed.

AATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCGG  >  NC_000913/3987028‑3987324
                                                                                                                                                   |                                                                                                                                                     
aaTTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCAt                                                                                                                                                     <  2:1082261/150‑1 (MQ=255)
 aTTTGCGGGGCAGTATCATGATCGTCACATTCAATAACATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATc                                                                                                                                                    <  2:788559/150‑1 (MQ=255)
  tttGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCa                                                                                                                                                   <  2:206201/150‑1 (MQ=255)
  tttGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCa                                                                                                                                                   >  2:1036031/1‑150 (MQ=255)
  tatGCTTGGCAGTATCCTGCTCGTCAAATTCAATAAGAAGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGGTTTCCGGCTTTGGTGTTTCCACCAGTTACGCATACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCa                                                                                                                                                   <  2:950637/148‑1 (MQ=255)
    tGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAgg                                                                                                                                                 >  2:601768/1‑150 (MQ=255)
    tGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAgg                                                                                                                                                 >  2:607909/1‑150 (MQ=255)
     gCTGGGCAGTATCATGATCGTCACATTCAATTAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGc                                                                                                                                                <  1:642845/150‑1 (MQ=255)
     gCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGc                                                                                                                                                <  2:1493127/150‑1 (MQ=255)
      cTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCCGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCCCGCGCCTGATCCATCAGGcc                                                                                                                                               >  1:1650929/1‑150 (MQ=255)
        gggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcg                                                                                                                                             >  2:1356953/1‑147 (MQ=255)
        gggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcg                                                                                                                                             <  2:1034474/150‑4 (MQ=255)
          gCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcgtg                                                                                                                                           <  1:1356953/150‑6 (MQ=255)
           cAGTTTCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATcc                                                                                                                                          <  2:494050/150‑1 (MQ=255)
             gTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCCCCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACTCAAACCTTCGCTACCGTTATCGGCACGCGCCTGTTCCATCAGGCCAATCCAt                                                                                                                                        <  1:233299/150‑1 (MQ=255)
             gTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCAt                                                                                                                                        <  1:143522/150‑1 (MQ=255)
                tCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGcc                                                                                                                                     <  2:222526/150‑1 (MQ=255)
                 cATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCt                                                                                                                                    >  2:1558533/1‑150 (MQ=255)
                   tGATCGTCACATTCAATAAGATGTTCAGCCATTGCCACCTGCAACGCTACCTGATGCCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTCCGCACACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATCCCtgt                                                                                                                                  >  2:788622/1‑150 (MQ=255)
                      tCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCtgttgt                                                                                                                               <  2:555996/150‑1 (MQ=255)
                        gTCACATTCAATAAGATGTTCCGCCATTGGCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTc                                                                                                                             <  1:1047601/150‑1 (MQ=255)
                         tCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGAGTTTTCCGGCTTTGGTTTTTCCACAAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGACAAACAAAGCATGATGTTcc                                                                                                                            >  1:1384069/1‑150 (MQ=255)
                               tCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCAtt                                                                                                                      >  2:900295/1‑150 (MQ=255)
                                cAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                     >  1:594362/1‑150 (MQ=255)
                                cAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                     <  2:1260485/150‑1 (MQ=255)
                                cAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCGTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCACCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                     <  2:1417922/150‑1 (MQ=255)
                                aaaTAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                     <  2:1602464/149‑1 (MQ=255)
                                 aaTAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGc                                                                                                                    >  2:188174/1‑150 (MQ=255)
                                 aaTAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGc                                                                                                                    >  1:959351/1‑150 (MQ=255)
                                  aTAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCa                                                                                                                   <  2:655094/150‑1 (MQ=255)
                                   tAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCAc                                                                                                                  <  1:620154/150‑1 (MQ=255)
                                     aGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGa                                                                                                                >  1:1472384/1‑150 (MQ=255)
                                        tGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGt                                                                                                             <  2:1624763/150‑1 (MQ=255)
                                        tGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGt                                                                                                             >  1:860368/1‑150 (MQ=255)
                                        tGTTCCGCCATTGCCACCTGAAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTACACCAGTTACGAAAACCTTCGCTACCGATATCGGCACGCGCCTGATCCATCAGGCCAATCGATGCCTGTTGTTACAGCATTGAACGATGt                                                                                                             >  1:187802/1‑150 (MQ=255)
                                             cGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGttcttc                                                                                                        >  1:1124116/1‑150 (MQ=255)
                                             cGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGttcttc                                                                                                        <  1:1571753/150‑1 (MQ=255)
                                              gCCATTGCCACCTGGAACGCTACCTGATGACGCGTTTTCCGGCGTTGGTTTTTCCACCAGATACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCGATGCCTGTTGTTCCAGCATTGCACGATGTTCTtca                                                                                                       <  2:117757/150‑1 (MQ=255)
                                              gCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTtca                                                                                                       <  2:652776/150‑1 (MQ=255)
                                              gCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTtca                                                                                                       <  2:1026895/150‑1 (MQ=255)
                                               ccATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTtcat                                                                                                      <  2:581980/150‑1 (MQ=255)
                                                cATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCTCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTtcatc                                                                                                     >  2:1055572/1‑150 (MQ=255)
                                                    gCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCa                                                                                                 <  2:1017281/150‑1 (MQ=255)
                                                       aCCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACa                                                                                              >  2:711533/1‑150 (MQ=255)
                                                           gCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATggg                                                                                          <  1:1510131/150‑1 (MQ=255)
                                                            cAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGc                                                                                         <  1:369452/150‑1 (MQ=255)
                                                             aaCGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCt                                                                                        >  2:1102024/1‑150 (MQ=255)
                                                              aCGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCtt                                                                                       >  1:369242/1‑150 (MQ=255)
                                                                 cTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCg                                                                                    >  2:608185/1‑150 (MQ=255)
                                                                       gatgaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTg                                                                              <  2:979637/150‑1 (MQ=255)
                                                                         tgaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGAt                                                                            >  1:1234238/1‑150 (MQ=255)
                                                                           aCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATgg                                                                          <  2:1387187/150‑1 (MQ=255)
                                                                            cgcgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATggg                                                                         >  1:930863/1‑150 (MQ=255)
                                                                            cgcgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATggg                                                                         <  1:601768/150‑1 (MQ=255)
                                                                            cgcgTTTTCCGGCTTTGGTTGTTCCTCCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATggg                                                                         <  1:1335598/150‑1 (MQ=255)
                                                                              cgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCAGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGat                                                                       >  2:617279/1‑150 (MQ=255)
                                                                                ttttCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGataa                                                                     <  1:1407970/150‑1 (MQ=255)
                                                                                  ttCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGataata                                                                   >  1:701254/1‑150 (MQ=255)
                                                                                    ccGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATc                                                                 <  1:1055572/150‑1 (MQ=255)
                                                                                      ggCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCca                                                               <  2:1513509/150‑1 (MQ=255)
                                                                                       gCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcag                                                              <  2:1644234/150‑1 (MQ=255)
                                                                                         tttGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcagca                                                            >  1:864833/1‑150 (MQ=255)
                                                                                                tttCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCTCGCGCCTGATCCATCAGGCCAAACCATGCCTGTTGTTGCAGCATTGCACGAAGTGCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACt                                                     <  2:254624/150‑1 (MQ=255)
                                                                                                   ccaccaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAAGTCCa                                                  >  2:166818/1‑150 (MQ=255)
                                                                                                   ccaccaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCATTCGCCATTGATGGGAGAATATCCAGCAGCGAACTCCa                                                  >  2:1215651/1‑150 (MQ=255)
                                                                                                     accaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATg                                                <  1:656264/150‑1 (MQ=255)
                                                                                                         gTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCAcc                                            >  1:1253242/1‑150 (MQ=255)
                                                                                                         gTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCAcc                                            <  1:619219/150‑1 (MQ=255)
                                                                                                            aCGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCCTGCCTGTTGTACCAGCATACCACGATGTTCTACATCACCAACATGGGCTTTCGCCATAGATGGGATAAAATCCAGCAGCGAACTCCATGCACCtgt                                         >  1:203048/1‑150 (MQ=255)
                                                                                                             cGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCtgtg                                        >  1:461/1‑150 (MQ=255)
                                                                                                              gCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGc                                       <  2:959351/150‑1 (MQ=255)
                                                                                                               cAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTAATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCg                                      >  1:1239204/1‑150 (MQ=255)
                                                                                                                 aaCCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACTATTTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGa                                    <  2:159722/150‑1 (MQ=255)
                                                                                                                  aCCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGat                                   >  1:632640/1‑150 (MQ=255)
                                                                                                                    cTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCCTGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGatat                                 >  2:1102861/1‑150 (MQ=255)
                                                                                                                      tCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATAc                               <  1:1652369/150‑1 (MQ=255)
                                                                                                                       cGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACg                              <  1:900295/150‑1 (MQ=255)
                                                                                                                           aCCGTTATCGGCACGCGCCTGATCCCTCAGGCCAATCCATGCCTGTTGGTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTg                          <  2:856112/150‑1 (MQ=255)
                                                                                                                           aCCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTg                          <  1:614105/150‑1 (MQ=255)
                                                                                                                            ccGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGt                         >  1:544950/1‑150 (MQ=255)
                                                                                                                            ccGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGt                         >  1:6120/1‑150 (MQ=255)
                                                                                                                              gTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTc                       >  1:840054/1‑150 (MQ=255)
                                                                                                                              gTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTc                       >  1:1613723/1‑150 (MQ=255)
                                                                                                                              gTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTc                       <  2:1071919/150‑1 (MQ=255)
                                                                                                                              gTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTc                       <  1:1671914/150‑1 (MQ=255)
                                                                                                                              gTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGAACCTGTGCGGATATACGCCTGTTc                       >  1:887021/1‑150 (MQ=255)
                                                                                                                                 aTCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGGCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGc                    >  1:163997/1‑150 (MQ=255)
                                                                                                                                    ggCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAg                 <  2:594362/150‑1 (MQ=255)
                                                                                                                                     ggattggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGa                >  1:901411/7‑150 (MQ=255)
                                                                                                                                      cACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGAc               <  2:930863/150‑1 (MQ=255)
                                                                                                                                      cACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGAc               <  2:1366660/150‑1 (MQ=255)
                                                                                                                                          cgcCTGATCCATCAGGCCAGTCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTaa           >  1:1405960/1‑150 (MQ=255)
                                                                                                                                          cgcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTaa           <  1:391042/150‑1 (MQ=255)
                                                                                                                                          cgcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTaa           >  1:587498/1‑150 (MQ=255)
                                                                                                                                           gcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAAt          >  1:1237063/1‑150 (MQ=255)
                                                                                                                                           gcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAAt          <  2:701254/150‑1 (MQ=255)
                                                                                                                                           gcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAAAATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGATCCGCCAGACGTaaa          >  2:1340752/1‑149 (MQ=255)
                                                                                                                                                  ccATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCg   <  1:1117373/150‑1 (MQ=255)
                                                                                                                                                   cATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCgg  >  2:635765/1‑150 (MQ=255)
                                                                                                                                                   cATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCgg  <  1:1036031/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCGG  >  NC_000913/3987028‑3987324

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