Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,987,174 C→G 5.5% D237H (GAT→CAT)  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,1740CG 5.5% 354.4 / 8.0 109D237H (GAT→CAT) hemYputative 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 = 8.79e-01

Only 100 of 111 total aligned reads displayed.

GATAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCG  >  NC_000913/3987025‑3987323
                                                                                                                                                     |                                                                                                                                                     
gatAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATc                                                                                                                                                       >  2:614105/1‑150 (MQ=255)
   aaTTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCAt                                                                                                                                                    <  2:1082261/150‑1 (MQ=255)
    aTTTGCGGGGCAGTATCATGATCGTCACATTCAATAACATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATc                                                                                                                                                   <  2:788559/150‑1 (MQ=255)
     tttGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCa                                                                                                                                                  >  1:206201/1‑150 (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)
              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)
                             cacaTTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCTCGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCa                                                                                                                          >  2:1117373/1‑150 (MQ=255)
                                  tCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCAtt                                                                                                                     >  2:900295/1‑150 (MQ=255)
                                   cAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                    >  2:602489/1‑150 (MQ=255)
                                   cAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                    <  2:1260485/150‑1 (MQ=255)
                                   cAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                    >  1:594362/1‑150 (MQ=255)
                                   cAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCGTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCACCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                    <  2:1417922/150‑1 (MQ=255)
                                   aaaTAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                    <  2:1602464/149‑1 (MQ=255)
                                    aaTAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGc                                                                                                                   >  1:959351/1‑150 (MQ=255)
                                    aaTAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGc                                                                                                                   >  2:188174/1‑150 (MQ=255)
                                     aTAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCa                                                                                                                  <  2:655094/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)
                                               ccGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGttctt                                                                                                        >  2:1671914/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)
                                                              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)
                                                                    cTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCg                                                                                   <  1:978730/150‑1 (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)
                                                                               cgcgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATggg                                                                        <  1:607909/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)
                                                                                   ttttCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGataa                                                                    <  1:1289584/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)
                                                                                            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:619219/150‑1 (MQ=255)
                                                                                                            gTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCAcc                                           >  1:1253242/1‑150 (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)
                                                                                                                           gCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGc                            <  2:749329/150‑1 (MQ=255)
                                                                                                                            cTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGcc                           <  1:836744/150‑1 (MQ=255)
                                                                                                                              aCCGTTATCGGCACGCGCCTGATCCCTCAGGCCAATCCATGCCTGTTGGTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTg                         <  2:856112/150‑1 (MQ=255)
                                                                                                                              aCCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTg                         <  1:1558533/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:1613723/1‑150 (MQ=255)
                                                                                                                                 gTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTc                      <  1:1671914/150‑1 (MQ=255)
                                                                                                                                 gTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTc                      >  1:840054/1‑150 (MQ=255)
                                                                                                                                 gTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGAACCTGTGCGGATATACGCCTGTTc                      >  1:887021/1‑150 (MQ=255)
                                                                                                                                  ttATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTcc                     <  2:1472384/150‑1 (MQ=255)
                                                                                                                                        ggattggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGa               >  1:901411/7‑150 (MQ=255)
                                                                                                                                         cACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGAc              <  2:930863/150‑1 (MQ=255)
                                                                                                                                             cgcCTGATCCATCAGGCCAGTCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTaa          >  1:1405960/1‑150 (MQ=255)
                                                                                                                                             cgcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTaa          >  1:587498/1‑150 (MQ=255)
                                                                                                                                              gcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAAt         <  2:701254/150‑1 (MQ=255)
                                                                                                                                              gcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAAt         >  1:1237063/1‑150 (MQ=255)
                                                                                                                                              gcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAAAATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGATCCGCCAGACGTaaa         >  2:1340752/1‑149 (MQ=255)
                                                                                                                                                     ccATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCg  <  1:1117373/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTCCG  >  NC_000913/3987025‑3987323

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: