Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,205,327 C→T 6.0% D386D (GAC→GAT
pseudogene (927/1137 nt)
ymfN →
aaaE →
e14 prophage; chimeric replication protein/phage terminase YmfN
phage terminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,205,3270CT 6.0% 262.0 / 4.2 83D386D (GAC→GAT
pseudogene (927/1137 nt)
ymfN
aaaE
e14 prophage; chimeric replication protein/phage terminase YmfN
phage terminase
Reads supporting (aligned to +/- strand):  ref base C (38/40);  new base T (3/2);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

CCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGC  >  NC_000913/1205178‑1205475
                                                                                                                                                     |                                                                                                                                                    
cccATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGAc                                                                                                                                                      <  1:41205/150‑1 (MQ=255)
cccATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGAc                                                                                                                                                      >  1:1076264/1‑150 (MQ=255)
cccATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGAc                                                                                                                                                      >  1:1336792/1‑150 (MQ=255)
  cATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACgg                                                                                                                                                    <  2:281898/150‑1 (MQ=255)
      gACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAAt                                                                                                                                                >  2:650348/1‑150 (MQ=255)
          ccTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCca                                                                                                                                            >  1:134310/1‑150 (MQ=255)
                   cGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCGCCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCt                                                                                                                                   <  2:280379/150‑1 (MQ=255)
                      ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGATGGGATTGCCGTCATGatgaaag                                                                                                                                <  1:1375215/150‑7 (MQ=255)
                      ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgt                                                                                                                                >  1:1074873/1‑150 (MQ=255)
                      ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgt                                                                                                                                <  1:1439135/150‑1 (MQ=255)
                      ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgt                                                                                                                                <  2:1090818/150‑1 (MQ=255)
                       ggtgCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgtt                                                                                                                               <  2:363606/147‑2 (MQ=255)
                       cGGGCTGTCACATGACCTTGCTGATGACGACCGGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgtg                                                                                                                               <  1:1121455/150‑1 (MQ=255)
                       cGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGgtg                                                                                                                               <  1:650348/150‑1 (MQ=255)
                             gTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCgg                                                                                                                         <  2:121940/150‑1 (MQ=255)
                              tCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGc                                                                                                                        <  2:525518/150‑1 (MQ=255)
                               cacaTGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCAATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCa                                                                                                                       >  2:1131504/1‑150 (MQ=255)
                                acaTGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCaa                                                                                                                      >  2:1184373/1‑150 (MQ=255)
                                  aTGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACg                                                                                                                    >  1:1563600/1‑150 (MQ=255)
                                   tGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGt                                                                                                                   <  1:989884/150‑1 (MQ=255)
                                      ccTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGt                                                                                                                >  2:953245/1‑150 (MQ=255)
                                      ccTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGt                                                                                                                >  2:82020/1‑150 (MQ=255)
                                         tGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCgg                                                                                                             >  2:552272/1‑150 (MQ=255)
                                           cTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCa                                                                                                           <  2:1074873/150‑1 (MQ=255)
                                              atgaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCaaaa                                                                                                        >  1:911887/1‑150 (MQ=255)
                                                gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                      <  1:179769/150‑1 (MQ=255)
                                                gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                      <  2:1552180/150‑1 (MQ=255)
                                                gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                      >  2:936740/1‑150 (MQ=255)
                                                gaAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAAcc                                                                                                      <  1:1146670/150‑1 (MQ=255)
                                                     ccTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTcc                                                                                                 <  1:1009933/150‑1 (MQ=255)
                                                      cTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCg                                                                                                >  2:188785/1‑150 (MQ=255)
                                                       tGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCgg                                                                                               >  1:662721/1‑150 (MQ=255)
                                                         aaCCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGt                                                                                             >  2:653669/1‑150 (MQ=255)
                                                           ccccGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTaa                                                                                           >  1:11326/1‑150 (MQ=255)
                                                             ccGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGAGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACg                                                                                         >  2:216310/1‑150 (MQ=255)
                                                              cGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACga                                                                                        <  2:371240/150‑1 (MQ=255)
                                                                 cACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACgatga                                                                                     >  1:617183/1‑150 (MQ=255)
                                                                   ccATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGatgatg                                                                                   >  2:572658/1‑150 (MQ=255)
                                                                   ccATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGatgatg                                                                                   >  2:1379944/1‑150 (MQ=255)
                                                                          tCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGaa                                                                            >  1:1378686/1‑150 (MQ=255)
                                                                           cAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAg                                                                           >  2:749051/1‑150 (MQ=255)
                                                                           cAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAg                                                                           >  2:659546/1‑150 (MQ=255)
                                                                                 tACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTc                                                                     <  2:1076264/150‑1 (MQ=255)
                                                                                  acacCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCa                                                                    >  2:1355528/1‑150 (MQ=255)
                                                                                   cgccAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAACACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCaa                                                                   >  2:763939/3‑150 (MQ=255)
                                                                                     ccAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                 >  1:1445788/1‑150 (MQ=255)
                                                                                     ccAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                 <  1:925269/150‑1 (MQ=255)
                                                                                     ccAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                 >  2:863121/1‑150 (MQ=255)
                                                                                     ccAACATGTCCGATCCGATGAAAGACCTGGAAGCGGCGATTGAATCGGAGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACAATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAgg                                                                 >  2:821862/1‑150 (MQ=255)
                                                                                        aCATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGc                                                              >  2:1481131/1‑150 (MQ=255)
                                                                                         cATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCa                                                             <  1:572658/150‑1 (MQ=255)
                                                                                           tGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAgg                                                           <  1:653669/150‑1 (MQ=255)
                                                                                           tGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAgg                                                           >  2:261457/1‑150 (MQ=255)
                                                                                              ccgatccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCgg                                                        <  1:936740/150‑1 (MQ=255)
                                                                                               cgatccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGa                                                       >  1:275878/1‑150 (MQ=255)
                                                                                               cgatccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGa                                                       <  2:322765/150‑1 (MQ=255)
                                                                                                  tccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGaaaa                                                    <  2:143988/150‑1 (MQ=255)
                                                                                                   ccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAAc                                                   <  2:911887/150‑1 (MQ=255)
                                                                                                   ccgatGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAAc                                                   <  1:650681/150‑1 (MQ=255)
                                                                                                       tGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACaaaa                                               <  1:1131504/150‑1 (MQ=255)
                                                                                                        gAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAAt                                              >  1:73245/1‑150 (MQ=255)
                                                                                                          aaGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCg                                            >  2:202618/1‑150 (MQ=255)
                                                                                                          aaGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCg                                            <  1:228232/150‑1 (MQ=255)
                                                                                                            gagCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCAGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGAt                                          <  2:270229/150‑1 (MQ=255)
                                                                                                             cgCTGGAAGCGGCGATTGAAGCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGCTGTATCGGCAACGTGGTCGGCACAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATg                                         <  2:517616/149‑1 (MQ=255)
                                                                                                                 ggAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGc                                     <  2:834433/150‑1 (MQ=255)
                                                                                                                  gAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCa                                    >  1:371052/1‑150 (MQ=255)
                                                                                                                   aaGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAg                                   <  2:1336792/150‑1 (MQ=255)
                                                                                                                       ggcgATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTgc                               >  1:249751/1‑150 (MQ=255)
                                                                                                                        gcgATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTgcg                              >  2:169360/1‑150 (MQ=255)
                                                                                                                               aaTCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGaaa                                                     <  1:1064662/120‑1 (MQ=255)
                                                                                                                               aaTCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTAAAGGAGCAGGCGGaaa                                                     >  2:1064662/1‑120 (MQ=255)
                                                                                                                                  cGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATgg                    <  2:89617/150‑1 (MQ=255)
                                                                                                                                   ggggCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGc                   >  1:977962/1‑150 (MQ=255)
                                                                                                                                   ggggCGCTTTCATCATGACGGCAATCCCATCATGACCTGGGGTATCGGCAACGTGGGCGGCAAAACCATTCCGGGTAACGATGAGGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGc                   <  2:1053463/150‑1 (MQ=255)
                                                                                                                                     ggCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCgg                 <  1:501600/150‑1 (MQ=255)
                                                                                                                                        gcTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTg              >  1:445927/1‑150 (MQ=255)
                                                                                                                                          ccaggtcatGATGGGATTGCCGTCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGACCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGc            >  2:1121955/6‑150 (MQ=255)
                                                                                                                                            aggtcatGATGGGATTGCCGTCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCag          >  2:965231/4‑150 (MQ=255)
                                                                                                                                              gtcatGATGGGATTGCCGTCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCagag        >  2:308260/2‑150 (MQ=255)
                                                                                                                                               tcatGATGGGATTGCCGTCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGc       <  1:308260/150‑1 (MQ=255)
                                                                                                                                                    aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc  <  2:1563600/150‑1 (MQ=255)
                                                                                                                                                    aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc  <  1:82020/150‑1 (MQ=255)
                                                                                                                                                    aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc  <  1:261457/150‑1 (MQ=255)
                                                                                                                                                    aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc  <  1:1355528/150‑1 (MQ=255)
                                                                                                                                                    aCGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGc  <  1:1184373/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCCGGGTAACGATGATGTGGTGAAGCCCGTCAAGGAGCAGGCGGAAAACAAAATCGATGGTGCAGTTGCGCTGATTATGGCGGTTGGCAGAGCCATGC  >  NC_000913/1205178‑1205475

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